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Abstract 


Normalization  and  Prediction 
of  Geotechnical  Properties 
Using  the  Cone  I'enetrometer  Test  (CPT) 
by 

Richard  Scott  Olsen 

Doctor  of  Philosophy  in  Civil  Engineering 
University  of  California  at  Berkeley 
Professor  James  K.  Mitchell,  Chair 


The  objectives  of  this  research  were  to  develop  techniques  for  (1)  stress 
normalization  of  CPT  measurements  (and  geotechnical  properties)  and  (2)  CPT 
prediction  of  geotechnical  properties  using  cone  and  sleeve  friction  resistance  values. 
Stress  normalization  allows  a  variable  geotechnical  property  to  be  reduced  to  an 
equivalent  vrlue  at  a  standard  confining  stress. 

A  new  concept,  the  Stress  Focus,  was  identified  which  provides  a  basis  for 
understanding  soil  strength  as  a  function  of  confining  stress.  This  study 
demonstrated  that  sand  friction  angles  for  different  initial  relative  densities  converge 
to  a  Stress  Focus  at  high  confining  stress  (approximately  100  atm),  where  the 
strength  behavior  is  similar  to  that  of  a  sedimentary  rock.  Dilation  of  dense  sands 
decreases  with  increased  confining  stress  until  the  Stress  Focus  is  reached,  as 
confirmed  using  historic  high  pressure  triaxial  test  data  as  weli  with  CPT 
measurements  from  laboratory  chamber  tests  and  uniform  soil  layers.  The  paths  of 
convergence  to  the  Stress  Focus  are  exponentially  related  to  confining  stress  and  are 
the  basis  for  development  of  CPT  cone  and  sleeve  friction  resistance  normalization 
techniques.  The  overburden  stress  at  the  Stress  Focus  is  soil  type  dependent. 
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The  stress  exponent  for  SPT  normalization  was  shown  to  be  equal  to  the  CPT 
derived  stress  exponent. 

CPT  correlations  to  geotechnical  properties  were  established  using  both  CPT 
cone  resistance  and  friction  ratio.  These  correlations  were  based  on  a  large  database 
which  was  developed  for  this  research  effort.  Statistical  evaluation  during  the 
development  of  these  correlations  concentrated  on  excluding  biased  data.  CPT 
based  correlations  were  established  for  the  following  geotechnical  properties:  SPT 
blow  count,  unconsolidated  undrained  triaxial  test  strength,  field  vane  shear  test 
strength,  and  shear  wave  velocity. 
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les  K.  Mitchell 
dissertation  Committee  Chairman 
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Chapter  1 


Introduction 


General 

Site  characterization,  in  terms  ot  geotechnical  properties,  can  be  the  single  most 
important  task  for  geotechnical  engineering  investigations.  Characterization  is 
defined  as  "To  uveal  and  separate  into  categories"  (Webster’s  New  Collegiate 
Dictionary,  1975).  Once  a  site  has  been  realistically  characterized  in  terms  of 
geotechnical  properly  distribution,  to  the  needs  of  a  project,  the  foundation  design  or 
foundation  performance  evaluation  (e.g.  liquefaction  evaluation)  can  be  done  with 
greater  economy  and  reliability.  Site  characterization  by  drilling  and  sampling  is 
economically  limited  because  of  the  small  proportion  of  the  subsurface  sampled  and 
the  relatively  high  cost  of  laboratory  tests. 

When  the  Cone  Penetrometer  Test  (CPT)  is  used  for  a  site  investigation,  the 
industry  standard  practice  for  CPT  data  reduction  has  generally  been  to  use  only  the 
CPT  cone  resistance  and  calculated  friction  ratio;  The  measured  dynamic  pore 
pressure  is  become  more  common  as  an  additional  CPT  measurement.  Most 
CPT-based  theories  are  derived  either  for  clay  or  sand  and  therefore  are  not  directly 
applicable  in  a  world  of  soil  mixtures,  dirty  sands,  and  silty  clays.  Use  of 
techniques  for  CPT  prediction  of  geotechnical  properties  that  arc  applicable  to  all 
soil  types  (Olsen  1984,  1988,  Olsen  and  Farr,  1980)  have  been  constrained  by  the 
limited  amount  of  published  verification. 
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The  goal  of  this  research  program  was  to  develop  improved  techniques  for  CP  F 
prediction  of  geotechnical  properties.  This  goal  was  achieved  by  creating  a  large 
database  (of  CPT  and  soil  test  results),  developing  new'  stress  normalization 
techniques,  accounting  for  bias  error  due  to  stratigraphy  changes  (between  CPT 
soundings  and  borings),  and  finally,  by  using  the  CPT  soil  characterization  chart 
(OJsen,  1988)  to  develop  predictive  correlations  for  Standard  Penetration  Test  (SPT) 
blow  count,  clay  undrained  strengths  (laboratory  unconsolidated  undrained  triaxial 
test  and  field  vane  shear  test),  and  shear  wave  velocity. 


The  CPT  Test 

The  CPT  test  is  performed  by  pushing  a  3.57  cm  (1.4  inch)  diameter 
instrumented  probe  into  the  earth  at  2  cm/second  while  simultaneously  measuring  the 
cone  resistance  and  sleeve  friction  resistance,  as  illustrated  in  Figure  1.1.  Cone 
resistance  (qc)  is  the  axial  component  of  the  stress  acting  on  the  tip  of  the  probe  (10 
cm2  horizontal  cross  sectional  area)  and  sleeve  friction  resistance  (fs)  is  the  sliding 
stress  developed  on  a  short  cylindrical  section  of  steel  just  above  the  tip  (surface  area 
of  150  cm2).  The  electric  CPT  has  received  wide  acceptance  throughout  the  world 
and  its  use  continues  to  grow.  The  advantages  of  the  electrical  CPT  are:  fast  rate  of 
penetration  (production  rate  up  to  600  ft  per  day),  liigher  accuracy  and  repeatability 
(compared  to  the  Standard  Penetration  Test  (SPT)),  and  use  of  a  computer  data 
acquisition  system  which  allows  for  computer  based  evaluation.  While  there  are 
numerous  additional  sensors  (Figure  1.1)  that  can  be  added  to  a  CPT  probe 
(Licberman,  ct.al,  1991,  Cooper,  et.al.,  1988,  Robertson,  et.al.,  1986,  FUGRO,  1980), 
approximately  60%  of  all  CPI  work  in  the  U.S.  is  done  using  only  the  two  basic 
measurements  (i.c.  qc  and  fs ).  The  CPT  is  also  becoming  the  geotechnical  tool  of 
choice  for  investigations  where  there  arc  local  CPT  contractors  and  the  site  is 
composed  of  clays,  sands,  or  soil  mixtures  containing  iittlc  or  no  gravel. 


Figure  1 
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CPT  measurements 

.1  The  Cone  Penetrometer  Test  (CPT)  system 


The  CPT  can  provide  a  large  quantity  of  low-rost,  repeatable,  well-distributed 
measurements  distributed  throughout  a  site.  Alternatively,  a  few  relatively 
undisturbed  soil  samples  may  be  tested  in  laboratory,  but  site  variability  generally 
overshadows  the  benefits  gained  from  a  few  specialty  tests.  CPT  is  becoming  the 
choice  in  situ  exploration  tool  over  the  SPT  because  of  the  lower  cost,  repeatability, 
continuous  record,  and  quantitative  nature  of  the  data.  The  CPT  is  also  easier  to 
interpret  because  the  measurements  are  more  fundamental.  However,  unlike  the  SPT, 
it  provides  no  soil  sample. 


The  CPT  Measurements 

The  two  CPT  measurements  are  remarkably  unique.  The  cone  resistance  is 
influenced  by  many  geotechnical  properties-it  varies  exponentially,  in  fact,  with  soil 
frictional  behavior.  For  example,  the  cone  resistance  can  be  200  tons/ft7'  (tsf) 

(20  MPa)  in  a  sand,  and  as  low  as  3  tsf  (300  KPa)  in  clays.  The  cone  resistance  can 
therefore  be  considered  an  index  of  the  sand  skeleton  strength  (Douglas  and  Olsen, 
1981).  The  CPT  sleeve  frictioh  measurement  is  a  high-strain  sliding  measurement 
along  a  steel  cylinder  (after  the  soil  has  navigated  around  the  cone  tip)  and  is  a  good 
indicator  of  loose  or  unstable  soil  structures  (Olsen  and  Farr,  1986).  While  the  cone 
resistance  has  been  the  topic  of  extensive  theoretical  and  experimental  research  (too 
numerous  to  reference),  there  has  only  been  limited  research  on  sleeve  friction 
resistance  (Olsen,  1984,  1988,  Olsen  and  Farr,  1986,  and  Masood,  1990). 


Developing  CPT  Prediction  Relationships  for  Geotechnical  Properties 

There  are  numerous  means  for  establishing  correlations  of  CPT  measurements 
with  geotechnical  properties.  Laboratory  chamber  testing,  while  very  useful,  has 
historically  been  done  only  on  clean  sands  (silty  sands  have  recently  been  tested  in 
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laboratory  chambers  (Rahardjo,  1989)).  Also,  the  results  of  CPT  chamber  tests  are 
sensitive  to  chamber  boundary'  conditions.  Theoretical  formulations  are  typically 
derived  only  for  a  particular  soil  type  and  must  always  be  thoroughly  verified  with 
laboratory  and/or  field  test  data.  Laboratory-based  correlations,  while  attractive  for 
pure  sands  and  pure  clays  may  be  of  limited  usefulness  in  evaluating  soils  composed 
of  mixtures  of  these  materials. 

Historically,  quantitative  CPT  correlations  with  geotechnical  properties  have  been 
based  only  on  the  cone  resistance.  CPT  sleeve  friction  resistance  was  used  only  for 
soil  classification  (Douglas  and  Olsen,  1981,  Robertson  and  Campanella,  1984). 
Researchers  have  typically  compared  CPT  cone  resistance  (qc)  and  laboratory 
strength  test  results  (or  field  strength  measurement  results)  in  an  X-Y  scatter  plot  to 
develop  best  fit  correlations.  CPT  prediction  of  geotechnical  properties  should  be 
based  on  combination  of  the  two  basic  CPT  measurements  (qc  and  fs)  because  these 
measurements  are  always  performed  and  each  reflects  different  aspects  of 
geotechnical  behavior. 


Stratigraphic  Influences  on  CPT  Correlations 

Using  field  data  to  correlate  CPT  results  with  geotechnical  properties  is  an 
attractive  approach,  however  there  arc  limitations.  Stratigraphic  variation  of  soil 
types  and  geotechnical  properties  between  a  CPT  sounding  and  a  nearby  borehole 
can  introduce  major  errors  into  a  correlation.  The  use  of  uniform  sites  only  would 
significantly  reduce  correlation  error,  however,  there  are  only  a  limited  number  of 
well-documented  uniform  sites  in  the  world.  These  uniform  sites  also  do  not 
represent  all  soil  types,  and  most  importantly  do  not  represent  a  wide  range  of 
strengths  (e.g.  relative  density  for  sand)  for  each  soil  type.  It  is  therefore  preferable 
to  represent  many  soil  types  and  relative  strength  levels  by  using  sites  with  less  than 
perfect  soil  uniformity.  Site  stratigraphy  (i.e.  uniformity)  must  still  be  accounted  for 
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when  developing  CPT  correlations.  A  subjective  quality  index  can  be  assigned  to 
each  CPT-to-boring  comparison  and  used  as  a  index  for  quantifying  bias  error 
potential  (e.g.  a  CPT  encounters  a  clay  layer  but  the  nearby  borehole  encounters  a 
sand  at  the  same  elevation). 


Accounting  for  Factors  that  Influence  Geotechnical  Properties 

Many  factors  influence  geotechnical  properties  and  should  be  accounted  for  when 
developing  CPT-based  correlations.  The  three  factors  of  strongest  influence  are: 

1)  soil  type,  2)  soil  state,  and  3)  confining  stress  level.  Soil  type  and  confining  stress 
are  self  explanatory.  The  soil  stale  (i.e.  void  ratio  at  a  standard  confining  stress 
level)  can  be  defined  as  the  normalized  geotechnical  property  level  (e.g.  Standard 
Penetration  Test  (SPT)  normalized  blow  count,  Nj).  Therefore,  a  CPT  technique  for 
predicting  geotechnical  properties  must  directly  or  indirectly  account  for  all  three  of 
these  influencing  factors. 


Using  CPT  based  Soil  Classification 

Soil  type  reflects  the  bulk  soil  composition  in  terms  of  sand  and  silt  grain  types 
and  sizes  (if  present)  together  with  the  clay  type  (if  present).  The  CPT  soil 
characterization  chart.  (Douglas  and  Olsen,  1981,  Olsen,  1984,  Olsen  and  Farr,  1986, 
Olsen,  1988)  can  approximately  predict  soil  type,  primarily  based  on  the  observation 
that  the  cone  resistance  is  exponentially  influenced  by  sand  grain  frictional  behavior 
and  is  therefore  an  index  of  the  soil  structure.  The  Olsen  (1988)  CPT  soil 
characterization  chart,  shown  in  Figure  1.2,  represents  the  state-of-the-art  for  CPT 
prediction  of  soil  type  from  knowledge  of  the  cone  resistance  (log  scale)  and  friction 
ratio  (which  is  defined  as  the  sleeve  friction  resistance  as  a  percent  of  the  cone 
resistance  level). 
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Figure  1.2  The  CPT  soil  characterization  (Olsen,  1988) 


Accounting  for  Stress  Effects  in  CPT  Correlations 


Confining  stress  influences  CPT  measurements  and  geotechnical  properties 
differently.  It  will  be  shown  that  confining  stress  influences  are  dependent  on  soil 
type  and  relative  strength.  The  influence  of  confining  stress  must  be  properly 
accounted  for  when  developing  correlations  of  CPT  to  geotechnical  properties. 

These  confining  stress  influences  are  also  important  dur  ing  prediction  of  geotechnical 
properties  using  field  CPT  data.  If  stress  influences  are  not  properly  considered, 
geo'.eclmical  properties  for  shallow  (e.g.  5  feet)  or  deep  (e.g.  200  feet)  situations  may 
be  over-  or  under-predicted  by  50%  or  more.  CPT  prediction  of  normalized 
geotechnical  properties  accounts  for  confining  stress  influences  on  properties. 

The  process  of  using  stress  normalization  to  predict  geotechnical  properties  is 
illustrated  in  Figure  1.3.  CPT  data  arc  initially  normalized  to  an  equivalent  value  at 
standard  vertical  effective  stress  (i.e.  atmospheric  pressure)  using  stress  exponent 
techniques.  The  normalized  geotechnical  property  is  then  predicted  using  the 
normalized  CPT  values.  Finally,  the  in  situ  property  is  determined  using  stress 
exponent  concepts.  The  stress  normalization  technique  for  CPT  measurements  and 
two  geotechnical  properties  (i.e.  SPT-N  value  and  shear  wave  velocity)  are  fully 
examined  in  this  dissertation.  Normalized  geotechnical  properties  will  be  derived 
based  on  tire  CPT  soil  characterization  chart. 
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CPT  Boring 


Figure  1.3  Using  o'  stress  normalization  for  CPT  prediction  of  geotechnical 
properties 
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Outline  of  Dissertation 


High  pressure  triaxial  testing  of  sands  and  the  effect  of  confining  stress  level  on 
Mohr  enve','  pe  curvature  are  reviewed  in  Chapter  2,  because  for  a  given  relative 
density,  the  friction  angle  decreases  with  increasing  confining  stress.  It  will  be 
shown  that  this  friction  angle  decrease  continues  until  a  specific  confining  stress  is 
achieved,  namely  the  Stress  Focus  (corresponding  to  an  approximate  depth  of  a  few 
thousand  feet).  Specimens  of  all  initial  relative  densities  for  a  particular  sand  type 
will  have  approximately  the  same  strength  (and  density)  at  the  Stress  Focus,  which 
corresponds  to  a  uncemented  sedimentary  rock  strength.  It  will  also  be  shown  that 
the  decrease  of  friction  angle  with  increasing  vertical  effective  stress  can  be 
explained  in  tenns  of  dilation  effects,  grain  crushing,  and  grain-to-grain  rolling 
influences. 


The  bar  i  stress  exponent  concept  required  for  stress  normalization  and 
determinatic  i  of  the  Stress  Focus  is  described  in  Chapter  3.  Stress  normalization  is 
achieved  usi  o  a  stress  exponent  on  the  vertical  effective  stress. 

CPT  com  resistance  normalization  formulation  for  all  soil  types  based  on  limit 
equilibrium  ai  d  cavity  expansion  theories  is  developed  in  Chapter  4,  This 
formulation  al.  o  accounts  for  the  observed  exponential  relationship  between  cone 
resistance  and  /ertical  effective  stress. 

CPT  laboratory  chamber  test  data  are  evaluated  in  Chapter  5.  This  evaluation 
supports  the  contention  that  cone  resistance  can  be  represented  by  stress  exponents 
and  the  Stress  Focus.  Standard  Penetration  Test  (SPT)  chamber  test  data  arc  also 
used  to  show  that  the  Stress  Focus  location  is  dependent  on  the  sand  classification. 
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CPT  data  from  uniform  in  situ  soil  layers  are  used  in  Chapter  6  to  establish  a 
technique  for  CPT  prediction  of  the  stress  exponent  (required  for  normalizing  CPT 
data).  These  data  also  show  that  the  Stress  Focus  location  depends  on  soil 
classification. 


New  stress  normalization  techniques  for  SPT  and  the  shear  wave  velocity  are 
given  in  Chapter  7.  For  the  SPT,  the  constant  drilling  mud  height  used  for 
laboratory  SPT  chamber  tests  produces  stress  exponents  that  are  too  low.  The 
CPT-determined  stress  exponent  (from  Chapter  6)  is  thus  recommended  to  be  used 
for  SPT  normalization.  Shear  wave  velocity  discussion  will  emphasize  the  need  for  a 
soil  type-dependent  stress  exponent  that  can  be  estimated  using  the  cone  resistance 
stress  exponent. 


In  Chapter  8,  correlations  are  given  for  normalized  CPT  parameters  with  the 
following  normalized  geotechnical  properties;  SPT  blow  count,  undrained  cohesive 
strength  based  on  the  unconsolidated  undrained  triaxial  test  (TxUU),  undrained 
cohesive  strength  based  on  the  field  vane  shear  (FVane),  and  shear  wave  velocity. 
These  correlations  are  represented  as  contours  on  the  CPT  soil  characterization  chart. 
Each  contour  was  established  using  an  Academic  Quality  Index  (AQ1)  to  subjectively 
remove  bias  due  to  stratigraphic  soil  type  changes  between  CPT  soundings  and 
boreholes.  Increasing  the  inclusionary  AQI  level  increases  the  data  group  quality  by 
removing  lower  quality  data. 


The  dissertation  conclusions  and  recommendations  for  future  research  arc  given 
in  Chapter  9. 
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Chapter  2 


Mohr  Envelope  Curvature  and 
Development  of  the  Stress  Focus  Concept 


Introduction 

Curvature  of  the  Mohr  failure  envelope  (i.e.  a  gradual  decrease  in  failure 
envelope  slope  with  increasing  effecti  ve  stress)  is  well-known.  It  is  incorporated  in 
many  nonlinear  behavior  theories  but  is  rarely  used  in  geotechnical  engineering 
practice.  Mohr  envelope  curvature  can  be  quantified  analytically,  but  there  are  no 
comprehensive  means  for  predicting  the  actual  curvature  parameters  based  on  soil 
type  and  relative  strength  (e.g.  friction  angle  for  sands).  At  present,  these  Mohr 
envelope  curvature  parameters  can  only  be  determined  with  a  series  of  triaxiai  tests 
conducted  using  a  wide  range  of  confining  stresses  for  a  given  sand  type  and  a  given 
initial  relative  density.  Mohr  envelope  curvature  and  the  Stress  Focus  concepts 
presented  in  this  chapter  are  the  foundation  for  the  CPT  stress  normalization 
techniques  in  this  dissertation. 


Baligh  Formulation  of  Curvature  Envelope 

Baligh  (1976)  imfroved  precedent  cavity  expansion  theory  by  incorporating  the 
effects  of  Mohr  failure  envelope  curvature.  Figure  2.1  summarizes  the  effects  of 
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failure  envelope  curvature  observed  by  Baligh  (1976)  in  terms  of  friction  angle  for 
several  types  of  sanus  and  relative  densities.  Baligh  (1976)  quantified  failure 
envelope  curvature  and  incorporated  its  effects  into  cavity  expansion  theory.  He 
proposed  liquation  (2.1)  to  express  failure  envelope  curvature  for  laboratoiy  triaxial 
tests. 


t  =  o 


tan<b  +  tana 


/  i  „  w 

1  ,  o 

xt-IoSio— 

V2-3 


(2.1) 


where 

4>0  =  Reference  friction  angle  at  0‘  =  1  ton/ft2  (tsf) 
a  -  Balighs’  Mohr  failure  envelope  curvature  parameter 

The  curvature  parameter,  a,  can  be  used  to  modify  the  reference  friction  angle,  <t>0, 
at  the  reference  stress  level,  ae,  as  shown  in  Equation  (2.2). 

tan^  =  tan<|>0  -  tana  LoglQav  (2.2) 

where 

o'v  =  vertical  effective  stress  (tsf) 

<t)b  =  friction  angle  at  0y 

This  curvature  parameter  (tan  a)  is  simply  the  semi-log  slope  for  given  relative 
density  trend  in  Figure  2.1.  Baligh  made  no  other  inference  concerning  the  trends 
in  Figure  2.1  concerning  sand  type  or  relative  density  toward  the  Mohr  envelope 
curvature  parameter,  a. 
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Figure  2.1  Failure  envelope  curvature  effects  for  different  relative  densities  and 
sand  types  (Baligh,  1976) 


Yareshenko  (1964)  proposed  an  expression  for  the  Mohr  failure  envelope 
curvature  of  the  form 

t  =  (ka)lln  <2-3> 

Baligh  (1976)  believed  that  Equation  (2.3)  was  elegant  and  convenient  for  analytical 
computations,  but  that  it  was  inadequate  for  practical  applications.  He  maintained 
that  there  were  many  combinations  of  k  and  n  that  could  fit  experimental  results 
because  they  were  obtained  from  log  x  versus  log  a  plots. 


New  Strength  Normalization  Based  on 
Failure  Envelope  Curvature 


A  new  expression  was  developed  by  the  author  (Equation  (2.4))  to  describe 
failure  envelope  curvature.  This  "stress  normalization"  technique  will  be  fully 
explained  in  Chapter  3.  The  normalized  shear  strength,  Tj,  in  Equation  (2.4)  is  the 
shear  strength  at  a  vertical  effective  stress  ot  1  atm.  The  proposed  exponent -based 
normalization  formulation  in  Equation  (2.4)  is  similar  to  that  of  Equation  (2.3), 
except  that  normalization  concepts  are  used. 


T 


(2.4) 


where 

(  ®  v)  atm 

X 

t 


Vertical  effective  stress  in  units  of  atm  pressure 
(approximately  tuns/fl2) 

Shear  strength  at  o[,  in  units  of  atm  pressure 

Shear  strength  at  o^  =  l  atm  (i.e.  tan(4>,)) 
shear  strength  stress  exponent 
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This  formulation  overcomes  the  limitation  stated  by  Baligh  (1976)  concerning  the 
Yareshenko  (1964)  equation,  namely  that  several  combinations  of  k  and  n  can  fit  the 
same  set  of  data.  For  loose  sands  or  clays  which  have  little  or  no  failure  envelope 
curvature,  the  Baligh  Mohr  envelope  curvature  parameter,  a,  is  approximately  zero 
and  the  stress  exponent  is  approximately  1  (or  slightly  less).  For  the  dense  Ottawa 
sand  in  Figure  2.1,  a  is  approximately  8°  and  the  stress  exponent  (l)  is  0.86.  Each 
of  the  constant  relative  density  trends  in  Figure  2.1  can  be  expressed  using 
Equation  (2.4)  with  a  constant  stress  exponent. 


The  Stress  Focus  Concept 

It  would  appear  from  Figure  2.1  that  the  only  means  to  determine  the  stress 
exponent  (e.g.  curvature  parameter)  is  with  a  thorough  laboratory  strength  testing 
program.  However,  close  examination  indicates  that  all  sands  have  approximately 
the  same  Motion  angle  at  some  high  overburden  stres:  between  70  and  300  atm 
irrespective  of  the  initial  relative  density.  The  overburden  stress  (afltl)  where 
specimens  of  all  initial  relative  densities  have  the  same  approximate  friction  angle 
(<(>f),  is  defined  as  the  "Stress  Focus".  The  Stress  Focus  is  therefore  a  "hinge"  to 
which  trends  associated  with  a  given  soil  type  and  all  relative  densities  converge  with 
increased  overburden  stress  as  illustrated  in  Figure  2.2.  This  Stress  Focus  occurs  at 
a  confining  stress  equivalent  to  a  depth  of  several  thousand  feet  below  groi  nd 
surface. 

Data  for  the  Chatahoochee  and  Sacramento  sands  in  Figure  2.1  were  replotted  as 
shown  in  Figure  2.3  and  Figure  2.4,  respectively,  to  illustrate  the  effects  of  failure 
envelope  curvature  and  the  resulting  Stress  Focus.  These  plots  depict  vertical 
effective  stress  on  the  vertical  axis  as  a  log  scale  with  the  values  increasing  in  the 
downward  direction.  The  dashed  least  square  fit  lines  are  also  shown  for  each  initial 
relative  density  group.  These  30  year  old  data  represent  the  best  quality  high  stress 
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NOTE: 
expanded 
scales  and 
simplified 
representation 


Overburden  Stress  at  the  Stress  Focus  = 
70  to  300  atm  (atmospheric  pressure) 
(e.g.  7  to  30  MPa  or  70  to  300  tsf) 


Figure  2.2  Simplified  representation  ol'  tire  Stress  Focus  using  the  Mohr  envelop 
curvature 
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tests  available  on  sand.  Other  published  data  do  not  encompass  the  range  of 
confining  stresses  and  the  number  of  relative  density  groups  represented  in  this  data 
set. 


Figure  2.3  and  Figure  2.4  indicate  that  at  the  Stress  Focus  the  friction  angle  is 
the  same  irrespective  of  the  initial  relative  density  from  which  the  samples  were 
consolidated.  Dense  sands  therefore  have  decreasing  friction  angle  with  increasing 
confining  stress  until  the  Stress  Focus  is  reached.  At  the  Stress  Focus,  initially 
dense  and  loose  sands  have  the  same  friction  angle.  Loose  sands  have  little  Mohr 
envelope  curvature  and  therefore  very  little  change  of  friction  angle  with  depth. 

When  the  figures  are  combined  (Figure  2.5),  the  result  suggests  that  the  Stress  Focus 
location  is  dependent  on  sand  type.  Stress  Focus  location  dependence  on  soil  type 
will  be  further  demonstrated  in  Chapters  4  through  6. 

Figure  2.3  and  Figure  2.4  are  plotted  in  terms  of  friction  angle  whereas 
Equation  (2.4)  is  expressed  in  terms  of  strength.  If  data  for  Sacramento  sand  (from 

Figure  2.4)  is  replotted  in  terms  of  drained  Mohr  Coulomb  strength  (i.e.  o^-tan(<j))) 

versus  a'v,  the  result  is  shown  in  Figure  2.6.  The  effect  of  Mohr  envelope 

curvature  in  Figure  2.6  is  not  as  obvious  as  in  Figure  2.4,  but  the  stress  exponent  (t) 
(i.e.  slope)  and  tj  can  be  determined  directly  from  this  figure  and  used  with 
Equation  (2.4). 


Rock  Mechanics  and  Friction  Angle 
at  High  Confining  Stresses 

The  Stress  Focus  (for  friction  angle)  appears  to  occur  at  a  overburden  stress 
range  between  70  and  300  atm.  Barton  (1976)  summarized  rock  strengths  at  high 
stresses  (greater  than  80  atm)  for  sedimentary  rocks  (e.g.  sandstone,  shale  and 
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Figure  2.5  Effect  of  sand  type  on  Stress  Focus  location 
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siltstone)  as  shown  in  Table  2.1.  This  table  indicates  that  the  friction  angle  (at  high 
confining  stress  levels)  increases  from  a  range  between  25°  and  30°  for  clay  based 
sedimentary  rock  (e.g.  chalk,  shale,  and  slate)  to  a  range  between  30°  and  34°  for 
quartz-based  sedimentary  rock.  The  friction  angles  at  the  Stress  Focus  (<j)f)  for 
Chatahoochee  and  Sacramento  sands  are  approximately  32°  and  31°,  respectively 
(from  Figure  2.3  and  Figure  2.4).  This  is  within  the  range  of  30  to  34°  for  quartz 
based  sedimentary  rock  at  high  confining  stress.  It  therefore  appears  that  sand  starts 
to  behave  similarly  to  a  sedimentary  rock  at  the  Stress  Focus. 


Table  2.1  High  stress  friction  angle  for  sedimentary  materials  (Barton,  1976) 


Material 

High  pressure  friction 
angle  (not  having 
dilative  behavior) 

Reference 

Sandstone 

31-34 

Coulson,  1972 

1 

River  sand 

30 

Vesic  &  Barksdale, 

(normal  stress  of 

30  to  630  atm) 

1964 

Siltstone 

27-32 

Coulson,  1972 

31 

Hobbs,  1970 

Shale 

27 

Ripley  &  Lee,  1962 

32 

Hobbs,  1970 

Chalk 

30 

Hutchinson,  1972 

Mudstone 

27 

Hobbs,  1970 

Slate 

25-30 

Barton,  1971 

— 

Dilation  effects  in  granular  materials  can  explain  the  friction  angle  versus 
confining  behavior  at  stresses  less  than  that  at  the  Stress  Focus.  Dilation  is  defined 
as  volume  increase  with  shear.  Figure  2.7  shows  that  the  observed  friction  angle 
during  shear  at  typical  geotechnical  confining  stress  is  caused  by  the  accumulative 
effect  of:  volume  increase  effects  during  shear  (i.e.,  dilation),  grain-to-grain 
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rearrangement  and  material  frictional  behavior  (<j>  ).  All  influences  except  the 
material  frictional  angle  (<J)  )  vary  with  confining  stress  level  (Mitchell,  1993).  This 
figure  illustrates  that  dense  sands  have  high  measured  friction  angles  due  to  the 
dilation  and  frictional  effects.  On  the  other  hand,  very  loose  sands,  at  the  critical 
void  ratio,  exhibit  a  friction  angle  due  only  to  grain-to-grain  rearrangement  and 
material  frictional  effects.  For  dense  sands  at  low  confining  stresses  (e.g.  1  atm), 
dilation  behavior  is  a  major  contributing  factor  toward  the  strength  level.  For 
example,  gravel  ballast  for  railroads  has  a  low  overburden  stress  but  very  high 
resistance  to  dilation.  Dilation  effects  therefore  increase  the  apparent  friction  angle 
for  dense  sands. 


Initially  dense  sand  appears  to  have  the  same  strength  as  initially  loose  sand  at 
the  Stress  Focus  ac  shown  in  Figure  2.5.  Very  loose  sands  do  not  dilate  at  any 
confining  stress  and  therefore  exhibit  no  dilation  effect  at  the  Stress  Focus.  Dense 


sands  at  the  Stiess  Focus  must  also  have  no  dilation  effect. 


vitvwio 


nro 
U IV 1  C1U1  V 


decrease  with  increased  depth  to  a  minimum  (or  near  zero)  level  when  the  Stress 
Focus  is  achieved  as  illustrated  in  Figure  2.8. 


The  density  for  different  initial  relative  densities  at  the  Stress  Focus  is  an 
important  issue.  Sand  density  will  inciease  with  increasing  consolidation  pressure. 
Initially  loose  sands  will  consolidate  (i.e.  density  increase  with  increasing  confining 
stress  level)  at  a  higher  rate  than  initially  dense  sands.  If  initially  loose  sands  and 
initially  dense  sands  have  the  same  strength  at  the  Stress  Focus,  then  the  densities 
may  also  be  equal  at  the  Stress  Focus.  Moiw  likely,  the  density  difference  between 
initially  loose  and  initially  dense  sands  at  the  Stress  Focus  are  probably  evident  but 
minor. 


The  Stress  ^ocus  represents  the  overburden  stress  where  soil-like  behavior 
becomes  (sedimentary)  rock-like  behavior.  It  also  provides  a  comprehensive  means 
of  relating  the  initial  relative  density  to  confming  stiess  level.  The  Stress  Focus 
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Porosity,  n  (%) 


Figure  2.7  Components  of  sh^ar  strength  in  granular  soils 
(After  Mitchell,  1 C93 ,  and  Rowe,  1962) 
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Proposed  interaction  of  friction  angle  effects  and  the  Stress  Focus 
concept 
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therefore  appears  to  be  a  fundamental  geotechnical  property  which  can  be  useful  for 
soil  behavior  evaluation  at  all  overburden  stress  conditions. 

For  confining  stresses  greater  than  the  Stress  Focus,  strength  is  achieved  by 
material  friction,  grain  crushing,  and  grain  rearrangement  effects  (which  can  be 
considered  as  a  minor  dilative  behavior  (Bruce,  et.al,  1966).  Shear  occurs  through 
the  grains  or  at  the  grain-to-grain  contacts  because  there  is  little  potential  for 
granular  rearrangement  at  high  confining  stress  levels  (Vesic  and  Clough,  1968). 
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Chapter  3 


Stress  Normalization 


Introduction 

Normalization  as  defined  by  the  Webster's  New  Collegiate  Dictionary  (1975)  is 
"to  reduce  to  a  norm  or  standard".  In  geotechnical  engineering,  stress  normalization 
allows  a  variable  geotechnical  property  to  be  reduced  to  an  equivalent  value  at  a 
standard  confining  stress.  The  easiest  obtained  standard  confining  stress  for 
geotechnical  stress  normalization  is  the  vertical  effective  stress.  It  is  the  only 
confining  stress  parameter  that  can  be  accurately  calculated.  The  standard  confining 
stress  unit  should  also  be  atmospheric  pressure  because  English  or  metric  units,  are 
arbitrary  stress  units. 

Stress  normalization  enables  an  equation  together  with  correlation  curves  to 
cover  a  wide  range  of  confining  stress  conditions  (analogous  to  dimensionless 
analytical  solutions).  Therefore,  normalized  CPT  measurements  can  be  used  to 
directly  predict  normalized  geotechnical  properties  even  though  the  confining  stress 
influences  on  the  CPT  parameters  and  on  the  geotechnical  properties  may  be 
different. 

This  chapter  will  describe  the  basic  stress  normalization  concept  using  stress 
exponents.  While  the  exponent  is  a  basic  mathematical  concept,  a  full  explanation 
using  geotechnical  engineering  terminology  will  simplify  the  discussion  in  this 
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dissertation.  The  Stress  Focus,  introduced  in  Chapter  2,  will  also  be  described  in 
terms  of  the  stress  exponent. 


The  Stress  Exponent  for  Geotechnical  Engineering 


CPT  measurements  (in  fact,  all  geotechnical  properties)  have  some  degree  of 
nonlinear  dependence  on  increasing  vertical  effective  stress.  Curvature  of  the  Mohr 
envelope  was  quantified  using  a  stress  exponent  in  Chapter  2.  When  the  behavior  of 
a  geotechnical  property  is  represented  with  an  exponent  of  the  vertical  effective 
stress,  die  representation  is  referred  to  as  stress  exponent-based  stress  normalization. 
Linear  and  exponential  behaviors  shown  at  the  top  of  Figure  3.1  become  straight 
lines  on  a  log-log  plot  as  shown  at  the  bottom.  Linear  behaviors  have  slopes  of  one 
to  one  and  exponential  behaviors  have  slopes  which  are  not  one  to  one.  The 
exponential  slope  is  equal  to  the  horizontal  stress  exponent  (e)  for  a  given  log  unit  of 
vertical  overburden  stress. 

An  exponentially  curved  line  in  Figure  3.1  becomes  a  straight  line  on  a  log-log 
plot  and  can  be  expressed  numerically  with  the  following  equation; 

H  =  C  V'  (3-1) 


v/here 

H 

= 

horizontal  ax's 

V 

— 

vertical  axis 

e 

= 

V  exponent 

C 

Constant 
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Figure  3.1 
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Representing  geotechnical  properties  in  terms  of  exponential  behavior 
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This  expression  can  be  rewritten  to  represent  the  horizontal  value  when  the  vertical 
axis  is  equal  to  1 : 


H 


1 


Ve 


(3.2) 


where 

H,  —  Horizontal  axis  value  v/hen  V  =  1 

The  exponent  in  Equation  (3.2)  must  be  specified  with  a  reference  such  as  the 
"V  exponent"  because  the  exponent  is  applied  to  V  and  the  slope  is  equal  to  the 
fraction  (or  multiple)  of  the  log10  H  for  one  log10  V  cycle.  This  V  exponent  (e)  is 
the  transformed  slope  on  the  log-log  chart,  namely  the  distance  e  for  a  given  V  as 
shown  in  Figure  3.1.  Lines  which  are  nearly  vertical  have  very  little  V  axis 
dependence  and  therefore  a  V  exponent  of  approximately  zero.  A  linear  relationship 
with  a  1:1  slope  (i.e.  45°  line)  has  a  V  exponent  of  1. 

In  geotechnical  engineering  the  vertical  axis  is  generally  taken  to  represent 
vertical  effective  stress  ( o  'v )  and  the  exponent  is  referred  to  as  the  stress  exponent. 
The  stress  exponent  is  thus  equal  to  a  fraction  or  multiple  of  the  Alog10(  property  ) 
over  one  log,0  o'v  cycle.  The  normalized  property  value  (e.g.  H,  in  Equation  (3.2)) 

is  simply  the  H  data  trend  projected  from  a  uniform  depth  zone  to  the  — 1 
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horizontal  line  as  shown  in  Figure  3.2.  Stress  exponent- based  stress  normalization 
for  geotechnical  engineering  is  therefore  defined  as: 


Hi 


(3.3) 


where 

o[, 

H, 

H 


vertical  effective  stress  (atmospheric  pressure  units) 
The  K  property  at  0^  =  1  atm 

The  H  property  at  the  o  'v  stress  level 
Stress  exponent 
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Figure  3.2  Illustrating  the  parameters  to  describe  the  curvature  of  the  failure 
envelope  for  geotechnical  engineering 
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The  Stress  Focus  and  Stress  Exponent  Interrelationship 


The  Stress  Focus,  introduced  in  Chapter  2,  represents  the  vertical  effective  stress 
where  all  trends  of  relative  density  (i.e.  curved  Mohr  envelopes)  intersect. 

Figure  3.3  is  a  generalization  of  the  Stress  Focus  concept  with  one  stress  exponent 
trend  shown.  Equation  (3.4)  is  a  simple  representation,  in  terms  of  log  scales,  of  the 
stress  exponent  line  in  Figure  3.3  which  can  be  simplified  as  shown  in 
Equation  (3,5). 

Log  H,  *  (log  -  log(  1 ))  e  =  Log  t5-4) 

Hj  -  (0/Hr  0-5) 

Equation  (3.5)  shows  that  the  H  level  at  Stress  Focus  (Hf)  is  equal  to  the 
normalized  value  (i.e.  H1  at  ov=  1  atm)  times  the  confining  stress  at  the  Stress 

Focus  (am)  raised  to  the  e  power.  Only  three  of  the  four  parameters  in 
Equation  (3.5)  are  required  to  define  an  expression.  For  example,  if  Hj  and  e  are 
known  and  am  can  be  reasonably  estimated,  then  the  H  level  at  the  Stress  Focus 
(Hf )  can  be  calculated. 


Stress  Type  and  Stress  Normalization 

Stress  normalization  using  the  vertical  effective  stress  will  be  shown  in  this  section 
not  to  introduce  error  even  when  a  particular  geotechnical  property  is  dependent  on, 
for  example,  lateral  stress.  CPI  measurements  and  geoteclinical  properties  are 
controlled  by  either  vertical,  mean,  octagonal,  or  lateral  stress.  This  stress  type 
merely  reflects  a  multiple  or  fraction  of  vertical  effective  stress.  For  example,  mean 
stress  is  equal  to  the  average  of  body  stresses  on  the  three  principle  axes  as  shown  in 
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LOG  Vertical  effective  stress  (Atm  units) 


1 


10000 


LOG  Measured  Value  (H) 


10  100  1000 


Figure  3.3  Stress  Focus  expressed  using  the  stress  exponent  and  normalized  value 
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Equation  (3.6); 


■  =  °v  *  °h  4  °h  (3.6) 

^mean  ~  2 

All  the  terms  in  Equation  (3.6)  can  be  defined  in  terms  of  vertical  effective  stress 
(oy)  as  shown  by  Equation  (3.7)  and  finally  by  Equation  (3.8). 


o 


mean 


+2  k; 

,  3  , 


(3.7) 


_  p 

°mean  ~  °v  1  mean 


(3.8) 


The  mean  stress  in  Equation  (3.8)  is  a  product  of  av  and  the  stress  factor,  Fmeajl; 


1  +  2K, 


(3.9) 


This  stress  factor  can  be  generalized  as  F.  with  the  *  representing  vertical,  mean, 
octagonal,  or  lateral  stress.  If  K0  does  not  change  appreciably  with  confining  stress 

level  then  F*  will  remain  constant  for  all  stress  levels.  The  o[,F,  expression  can  be 
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substituted  into  the  basic  normalized  stress  equation  (Equation  (3.3))  with  the  result 
shown  below; 


*1 


which  equals; 


H 


1 


where 

(f.r 


(3.10) 


(3.11) 


(3.12) 


Equation  (3.11)  is  now  enual  to  the  basic  stress  normalization  expression  (Equation 
(3.3))  with  Hj  in  terms  of  F*  (Equation  (3.12)).  If  K0  (i.e.  Fmean)  does  not  change 
with  vertical  efiective  stress  level,  then  H,  is  a  constant  (for  a  given  soil  type  and 
relative  strength).  Correlations  by  Kulhaway  and  Mayne  (1990)  have  shown  that  K0 
(e.g.  Ftncan)  changes  very  little  with  increasing  oj,.  Therefore,  die  stress  type 

(mean,  lateral,  etc.,)  or  K0  are  not  required  for  stress  normalization  as  long  as  the 
overconsolidation  level  (or  normally  consolidation  condition)  remains  constant  for  all 

stress  levels.  As  a  result,  <jv  can  be  used  as  the  normalization  stress  type  even  if 

the  geotechnical  property  is  controlled  by  another  confining  stress  factor  (e.g.  mean 
or  lateral  stress). 
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The  Atmospheric  Pressure  Standard 


Standard  atmospheric  pressure  provides  a  convenient  basis  for  expressing 
geotechnical  properties.  There  is  only  one  standard  stress  level  and  that  is  the 
atmospheric  stress— any  other  stress  units  such  a  KPa  or  psi  (pounds/inch2)  are  based 
on  arbitrary  units.  Atmospheric  pressure  (atm)  is  not  a  usual  engineering  pressure 
level  but  is  very  close  to  the  English  ton/foot2  (tsf)  or  one  kilogram  force  per 
cm2  (Kgf/cm2 )  or  bars  units  as  shown  below: 

1  am  =  1.058—  =  1.033  Kg/orce  =1.013  bar  =10132  KPa  (3.14) 

ft1  cm 1 


The  conventional  geotechnical  engineering  method  for  expressing  normalized 
parameters  is  shown  in  Equation  (3.15)  with  Pa  equal  to  the  atmospheric  pressure. 


£1 

,  \ 

Oy 


(3.15) 


Equation  (3.15)  adds  an  additional  factor  of  complexity  to  a  simple  concept. 
Moreover,  if  Equation  (3.15)  expressed  in  terms  of  Hj  the  following  is  produced; 


(3.16) 


This  final  expression  is  more  complicated  than  Equation  (3.3)  for  stress 
normalization  because  stress  units  are  not  atmospheric  pressure  units.  Therefore, 
atmospheric  pressure  (atm)  is  used  for  stress  units  throughout  this  dissertation  to 
keep  the  stress  notation  simple  and  because  normalization  concepts  will  become  more 
complex. 
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Chapter  4 


Developing  the  Cone  Resistance 
Normalized  Formulation 


Introduction 

A  cone  resistance  normalization  formulation  allows  the  equivalent  cone  resistance 
at  q  'v  =  1  atm  to  be  determined.  Cone  resistance  normalization  is  an  integral  part 

of  the  process  toward  predicting  geotechnical  properties,  it  possible,  it  should  have  a 
theoretical  background  rather  than  only  match  the  trend  of  laboratory  data.  The 
formulation  should  be  simple  arid  not  require  elaborate  laboratory  test-based 
properties  as  input  to  the  formulation.  It  should  also  account  for  exponentia' 
behavior  using  the  stress  exponent  concept  from  the  previous  chapters. 

A  CPT  cone  resistance  normalization  formulation,  based  on  generalization  of 
bearing  capacity  formulations,  is  developed  in  this  chapter,  The  bearing  capacity 
formulation  considers  both  limit  equilibrium  techniques  (for  shallow  bearing  failures) 
and  cavity  expansion  theory  (for  deep  bearing  failures).  The  observed  non-linearity 
of  cone  resistance  with  vertical  effective  stress  will  also  be  introduced  and  analyzed 
further  in  Chapters  5  and  6. 


Comparison  of 

Limit  Equilibrium  and  Cavity  Expansion  Theories 


Most  theories  for  the  bearing  capacity  of  both  deep  foundations  (e.g.  piles)  and 
penetrometers  have  severe  shortcomings.  During  the  initial  half  the  20Ih  century, 
bearing  failure  were  evaluated  using  limit  equilibrium  (i.e.  surface  expression) 
theories.  These  classical  limit  equilibrium  theories  assume  soil  movement  upward  to 
the  ground  surface  using  resisting  stresses  on  wedges  (Ve sic,  1972,  Terzaghi,  1943) 
as  shown  at  the  left  of  Figure  4.1.  Limit  equilibrium  theories  are  based  only  on 
Mohr  envelope  strength  parameters  (i.e.  c  and  <(>).  In  the  1 960's  and  1 970’s,  cavity 
expansion  theories  (Vesic,  1963,  1972)  were  shown  to  better  describe  bearing  failure 
at  great  depth.  Cavity  expansion  theory  assumes  a  hypothetical  spherical  or 
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right  of  Figure  4.1.  Cavity  expansion  theories  are  based  on  non-linear  strength 
descriptors  and  volume  change  parameters.  The  cone  resistance  behavior  in 
laboratory  chambers  has  also  been  computer  modeled  using  cavity  expansion  theory 
together  with  hyperbolic  modeling  of  the  elastic  zone  (Tseng,  1989). 


Exponential  Behavior  of  the  Cone  Resistance 

Classical  limit  equilibrium  (surface  expression)  theories  esult  in  bearing  stresses 
which  arc  linearly  proportional  to  vertical  effective  stress  (Terzaghi  and  Peck,  1967, 
Durgunoglu  ana  Mitchell,  1973,  1975).  Cavity  expansion  theories,  on  the  other 
hand,  produce  a  non  -linear  exponent-based  relationship  of  cone  resistance  with 
overburden  stress  (Vcsic,  1972).  Cavity  expansion  theories  also  show  more 
exponential  curvature  when  soil  compressibility  and  Mohr  envelope  curvature  effects 
arc  included  (Baligh,  1975). 
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Limit  Equilibrium 
for  Shallow  depths 

{  O  1  t  y\  «-v>  /v«yi4- 
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toward  ground 
surface) 


The  transition 
between  Limit 
Equilibrium  and 
Cavity  Expansion 
is  called  the 
Critical  Depth 


Cavity  Expansion 
at  depth 


Figure  4.1  Comparison  of  lire  Classical  Limit  Equilibrium  based  bearing  failure  to 
Cavity  Expansion  bear  ing  failure 
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Combining  a  curved  Mohr  envelope  effect  with  the  simplistic  limit  equilibrium 
bearing  formulation,  as  illustrated  in  Figure  4.2,  can  also  create  a  decreasing, 
non  linear  exponential  dependence  of  cone  resistance  on  vertical  effective  stress. 
However,  the  magnitude  of  curvature  is  less  than  observed  in  laboratory  chamber 
tests.  Therefore,  the  effects  of  Mohr  envelope  curvature  provide  only  part  of  the 
answer  toward  defining  cone  resistance  to  vertical  effective  stress  behavior. 


It  will  be  shown  that  an  exponent  of  the  vertical  effective  stress  (introduced  in 
Chapter  3)  can  characterize  the  curved  cone  resistance  behavior  observed  from 
chamber  tests  (in  Chapter  5)  and  from  CPT  tests  in  uniform  soil  layers 
(in  Chapter  6).  This  stress  exponent  is  based  on  observed  data  and  accounts  for 
contributing  factors  such  as  Mohr  envelope  curvature,  cavity  expansion  effects, 
compressibility,  grain  crushing  effects,  and  others.  The  cone  resistance  normalization 
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cl'id  in  Sitii 


data  in  later  chapters. 


Figure  4.2  Approximately  combining  the  classical  limit  equilibrium  bearing 
capacity  formulation  with  non  linear  Mohr  envelope  behavior 


43 


Deveiopment  of  the 

Cone  Resistance  Normalization  Formulation 


Limit  equilibrium  and  cavity  expansion  bearing  theories  are  similar  in  that  they 
both  have  cohesive  and  frictional  bearing  factor  components.  The  classic  limit 
equilibrium  formulation  for  bearing  failure  stress  is  shown  in  Equation  (4.1) 
(Terzaghi,  1943): 


<lo 


oNC  +  ~ 

<?  O 


*ya;c, 


(4.1) 


where 

% 

Nc,  Ny,  Nq 

Sc  ^  ,  S, 


c 

o 

7 

B 


Bearing  capacity 

Bearing  capacity  factors  as  a  function  of 
Bearing  capacity  shape  factors  which  depend  on  the  soil 
angle  of  internal  friction  (during  shear)  and  shape  of  the 
bearing  surface 

Cohesion  (or  shear  strength  intercept) 

Confining  stress 
Soil  unit  weight 
Strip  footing  width 
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The  classic  spherical  cavity  expansion  formulation  is  shown  in  Equation  (4.2) 
(Vesiu,  1972). 


K  =  cf; 


°nFq 


(4.2) 


where 

* 

Pl 

:= 

Cavity  expansion  pressure 

f'c  and  F* 

= 

Dimensionless  cavity  expansion  factors 

= 

Effective  mean  normal  stress 

c 

= 

Soil  cohesion 

F 

_ 

_  ,  ,  4sin<t> 

J  ( 1  +  sin<J)  )  /  -r  \  3(1  ♦  sin<t>) 

3  -  sin4>  rr 

q 

Fc 

(  Fq  -  Dcottj) 

In 

h 

i  x  r  A 

-i-  a.  y 

Ir 

= 

Rigidity  Index  =  /— 1  =  —  — 7 — - - - 

b  J  \SJ  S  2  (1  +  v) 

A 

= 

Average  compressive  volumetric  strain  in  plastic  zone 

G 

= 

Shear  modulus 

E 

= 

Young's  modulus 

S 

= 

Soil  shear  strength  (c  +  a'tamj)) 

The  cavity  expansion  pressure  is  the  pressure  required  to  expand  a  cavity  in  situ 
(i.e.  pressure  required  to  expand  a  balloon  at  some  depth).  Cavity  expansion 
pressure  is  related  to  bearing  pressure  because  a  cavity  must  be  opened  for  the 
penetrometer  to  advance.  While  cavity  expansion  is  simple  in  concept,  converting  it 
to  bearing  pressure  is  not  a  trivial  task.  The  cavity  expansion  pressure  level  is  still 
related  to  the  bearing  pressure  level,  in  general. 
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Both  of  these  formulations  can  be  summarized  for  penetrometers  as  shown  in 
Equation  (4.3)  using  the  Zc  and  Z0  bearing  factors. 

qc  =  zcc  +  zo° 


where 

4c 

_ 

bearing  stress 

Zc 

= 

cohesive  bearing  factor 

Zc 

= 

frictional  bearing  factor 

c 

= 

cohesion  strength 

cr 

= 

overburden  stress 

The  bearing  stress  (q0)  or  cavity  pressure  (p*)  was  replaced  by  qc  in  Equation 

(4.3)  for  the  cone  penetrometer  representation.  The  Z-ByN^  component  in 
Equation  (4.1)  can  be  ignored  because  the  penetrometer  width  (B)  is  very  small 
compared  to  the  depth.  At  any  depth  the  Zc  and  Z0  can  be  determined  (theory 
dependent)  and  the  bearing  stress,  qc ,  may  be  calculated.  The  Z0  and  Zc  bearing 
factors  are  dependent  on  friction  angle  and  cohesive  strength  for  limit  equilibrium  or 
on  several  geotechnical  parameters  for  the  cavity  expansion  theory.  The  principal 
difference  between  the  limit  equilibrium  and  cavity  expansion  theories  lies  in  how 
the  Z  bearing  factors  are  determined.  The  ZCT,  for  example,  is  as  simple  as  the 
equivalent  Nq  factor  (Terzaghi,  1943)  which  is  based  only  on  the  friction  angle,  or  ZG 
can  be  based  on  numerous  geotechnical  properties  which  are  also  influenced  by 
stress  level  (Vesic,  1972,  Baligh,  1976). 

The  following  discussions  separately  derive  a  normalized  cone  resistance 
expression  for  clay  and  sand.  Nonlinear  exponential  elfccts  r  z  then  be  included 
using  the  stress  exponent  concept.  Ihc  final  step  is  a  normalized  cone  resistance 
formulation  for  all  soil  types. 
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Development  of  the  Bearing  Capacity  Formulation  for  Clay 

The  clay  based  normalized  cone  resistance  can  be  derived  from  the  generalized 
bearing  capacity  formulation  in  Equation  (4.3).  Penetration  of  clay  occurs  as 
undrained  behavior.  Therefore,  the  bearing  stress  for  quickly  loaded  saturated  clay  is 
not  influenced  by  confining  stress.  As  a  result,  the  ZCT  frictional  bearing  factor  is 
equal  to  one  (1)  (Terzaghi  and  Peck,  1967)(Vesic,  1972)  which  reduces  Equation 
(4.3)  to: 

a  =  Z  S  +  o  ,  (4.4) 

The  cohesive  strength  (c)  in  Equation  (4.3)  was  replaced  with  the  conventional 
notation  Su  for  undrained  shear  strength.  The  next  step  is  to  define  the  net  cone 
resistance.  If  the  clay  deposit  has  zero  strength  (i.e.  Su=0)  then  Equation  (4.3) 
implies  that  buoyant  total  force  will  act  upward  on  the  cone  bearing  surface: 

4c  =  Ototal  <4-5) 

The  net  cone  resistance  shown  in  Equation  (4.6)  is  that  remaining  after  removal  of 
buoyancy. 

K)m  =  9c  ~  ° total  (4-6) 

Equation  (4.4)  can  now  be  rearranged  as  shown  using  the  net  cone  resistance: 

(«*)*«  =  4C  ~  owuU  =  ZcSu  (4.7) 

The  net  cone  resistance,  (qc)nct ,  is  now  only  dependent  on  the  soil  strength.  The 
next  step  requires  dividing  both  sides  of  Equation  (4.7)  by  the  normalizing  vertical 
effective  stress,  q'v,  to  produce  the  following: 
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4  c  °  total 


=  z_~ 


(4.8) 


The  normalized  clay  undrained  strength,  — 7  ,  in  Equation  (4.8)  is  often  represented 

cm 


by  which  can  also  be  expressed  using  normalization  terminology  as  Sul: 


£ 

P 


(4.9) 


Equation  (4.8)  «  an  now  be  expressed  as  the  normalized  expression  for  clay  as: 


4f 


o 


V 


(4.10) 


Development  of  the  Bearing  Capacity  Formulation  for  Sand 

Bearing  stress  in  sands  is  entirely  a  function  of  die  effective  stress,  since  sands 
are  frictional  materials.  Therefore,  cohesive  undrained  strength,  c,  is  zero  and 
Equation  (4.3)  simplifies  to: 

qc  =  o'v  Zo  (4-H) 

If  a  sand  has  no  shear  strength  (i.e.  4>— 0),  then  the  frictional  bearing  factor  (Za  )  is 
equal  to  one  (i.e.  Nq=l)  (Terzaghi  and  Peck,  1967)(Vesic,  1972).  Therefore,  qc 
equals  o(,  according  to  Equation  (4.11)  because  N  =1.  However,  for  a  <j)=0 

condition, 0 cotal=u  (hydrostatic  stress)  which  means  o^=0,  which  translates  to  qc=0 
according  to  Equation  (4.1 1)  which  is  incorrect  (it  siiouia  instead  equal  a  buoyancy 
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stress).  The  qc  in  Equation  (4.1 1)  must  therefore  reflect  a  net  cone  resistance  similar 
to  that  shown  in  Equation  (4.6)  for  bearing  stress  of  clay.  Equation  (4.1 1)  can  now 
be  rewritten  using  the  net  cone  resistance  as  shown  below: 


Qc  ° total  °V  Z, 

with  rearrangement,  the  normalized  expression  for  sand  is: 


Qc  ~  0 


total 


=  [Z.  1 


0  J sand 


(4.12) 


(4.13) 


Modification  to  the  Bearing  Formulation  for 
Non-linear  Exponential  Behavior 


The  vertical  effective  stress,  o^,  in  the  denominator  ol  Equations  (4.1U) 

and  (4.13)  is  the  normalizing  stress.  This  av  requires  a  stress  exponent  to  reflect  the 
observed  non-linear  behavior  of  cone  resistance  with  vertical  effective  stress  (concept 
introduced  in  Chapter  2  and  to  be  illustrated  in  Chapter  5).  The  resulting  normalized 
cone  resistance  shown  in  Equation  (4,14)  is  equal  to  the  left  side  of  either  Equations 
(4.10)  or  (4.13)  with  a  stress  exponent  (c)  included.  The  normalized  cone  resistance 
(qcle)  subscripts  have  specific  representations:  "c"  for  cone  resistance,  "1"  for 
normalization  to  an  equivalent  value  at  a  vertical  effective  stress  of  1  atm,  and  "e"  to 
represent  that  normalization  accounts  for  exponential  curvature. 


(4.14) 


where 

c  =  Cone  resistance  stress  exponent 
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Bearing  Stress  Formulation  for  All  Soil  Types 


The  normalized  cone  resistance  expression  for  all  soil  types  is  the  combined 
effect  of  the  cohesive  bearing  stress  expression  from  Equation  (4.10)  and  the 
frictional  sand  bearing  stress  expression  from  Equation  (4.13)  as  shown  below: 

-  9,,,  -  Su,Zk  -  Z,  (4.15) 

K) 

This  new  stress-normalized  expression  is  remarkably  similar  in  appearance  to  the 
Terzaghi  bearing  formulation  (expressed  for  deep  penetrations)  as  shown  below: 

q  =  c  N*  +  o  N"  (4*16) 

Equation  (4.15)  is  expressed  in  terms  of  the  normalized  cone  resistance,  normalized 
undrained  cohesive  strength,  and  bearing  factors.  It  represents  a  new  technique  for 
cone  resistance  normalization  that  will  be  validated  using  laboratory  chamber  data 
in  Chapter  5  and  using  uniform  soil  layer  data  in  Chapter  6.  It  is  also  the  basis  for  a 
new  technique  for  prediction  of  undrained  strength  of  clays  in  Chapter  8. 

It  will  also  be  shown,  in  later  chapters,  that  the  bearing  factors  and  the  stress 
exponent  in  Equation  (4.15)  are  all  related  to  the  combination  of  void  ratio  anc  soil 

type  at  0'  =1  atm.  The  bearing  factors  (i.e.  Zk  and  Zq)  define  the  mechanical 

(strength)  behavior  at  a  c'  =1  atm.  And  the  stress  exponent  (c)  defines  the 

exponential  stress  curvature  of  cone  resistance  to  vertical  effective  stress.  Therefore, 
the  bearing  factors  define  mechanical  strength  behavior  and  the  shess  exponent 
defines  the  curvature  of  cone  resistance  to  vertical  effective  stress.  Chapters  '  a  id  6 
will  specifically  show  that  increasing  the  initial  relative  density,  of  a  given  sand, 
creates  a  higher  strength  (at  o^.=l  atm)  and  a  lower  stress  exponent  (c). 
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Chapter  5 


Cone  Resistance  Stress  Focus 
Confirmation  Using 
CPT  Chamber  Test  Data 


Introduction 


The  Stiess  Focus  concept,  introduced  in  Chapter  2.  can  also  be,  used  10  describe 
the  relationship  of  CPT  cone  resistance  with  ovcrbuiden  stiess.  A  cone  resistance 
Stress  Focus  concept  can  be  used  to  generalize  immeious  effects  such  as  cavity 
expansion,  Mohr  envelope  curvature,  grain  crushing,  compressibility,  and  others. 

Stress  Focus  will  be  confirmed  based  on  cone  resistance  measurements  obtained  from 
CPT  laboratojy  la’ge  diameter  chamber  tests.  As  additional  confimvilio.,  of  the 
dress  Focus  for  in  situ  penetrometers,  the  Standard  Penetration  Test  (SP’J)  chamber 
t  resuits  will  also  be  examined.  It  will  be  shown  in  this  chapter  that  the 
ionship  of  cone  resistance  to  vertical  effective  stress  observed  in  laboratory 
cluiii::  r  tests  become  straight  lines  when  plotted  logarithmically  indicating  a 
constant  .  >re$s  exponent.  It  will  also  be  shown  that  the  trends  of  all  relative  densities 
on  this  logarithmic  plot  converge  to  the  Stress  l  ocus  at  great  depth  as  illustrated  in 
Figure  5.1 . 
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Figure  5.1  Proposed  relationship  of  eonc  resistance  versus  vertical  effective  stress 
trending  for  a  given  sand  (or  soil)  type 
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It  is  proposed  to  show  in  this  chapter  that: 


1)  Cone  resistance  may  be  expressed  as  a  function  of  vertical  effective  stress  for 
a  given  initial  relative  density  using  a  constant  stress  exponent  (indicating  a 
straight  line  on  a  logarithmic  plot) 

2)  The  stress  exponent  decreases  with  increased  i  nitial  relative  density  n  sands 

3)  The  relationships  between  cone  resistance  and  vertical  effective  stress 
(regardless  of  initial  relative  density)  converge  at  great  depth  (which  is  defined  as 
the  cone  resistance  Stress  Focus) 

4)  Overconsolidation  causes  a  lower  stress  exponent 


Interpretation  of  Large  Diameter  Test  Chambt:  Data 

The  best  means  of  introducing  and  defining  the  components  of  the  cone 
resistance  Stress  Focus  is  with  a  good  example.  Figure  5.2  (Baldi,  et  al.,  1 98 1 ) 
summarizes  relative  density  contours  for  Ticino  sand  derived  from  a  large  number  of 
chamber  tests  at  various  initial  relative  densities  and  confining  stress  levels.  These 
contours  were  carefully  established  by  the  original  researchers  for  each  relative 
density  range.  The  same  contours  arc  replotted  in  terms  of  logU)  net  cone  resistance 
(e.g.  (qc)ncI )  versus  log]0  in  Figure  5.3  (net  cone  resistance  is  the  measured  cone 

resistance  minus  the  vertical  total  stress  as  described  in  Chapter  4).  The  relationships 
shown  for  different  initial  relative  densities  in  Figu'f*  5.3  all  converge  to  a  common 
Stress  Focus  for  a  given  sand  type.  Relative  density  contours  in  Figure  5.3  are 
remarkably  straight  in  spite  of  the  fact  that  the  (qc)ncl  versus  o  'v  relationship 

(Figure  5.2)  were  produced  by  the  original  authors  from  scattered  test  data  without 
foreknowledge  of  the  Stress  Focus  concept.  These  relationships  have  also  been 
verified  by  rcplotting  the  original  data  as  shown  in  a  later  section  of  this  chapter. 

The  lines  foi  all  relative  densities  converge  to  the  Stress  l  ocus  because  cavity 
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expansion.  Mohr  envelope  curvature,  grain  crushing,  and  compressibility  influences 
are  changing.  A  sand  of  initially  high  relative  density  (dense)  will  thus  have  the 
same  cone  resistance  (and  density)  as  a  initially  low  relative  density  (loose)  sand  at 
the  depth  of  the  Stress  Focus. 


Explanation  of  the  Cone  Resistance  Stress  Focus 


The  cone  resistance  at  the  Stress  Focus  (qcf)  in  Figure  5.3  is  at  2200  atm,  and  the 
vertical  effective  stress  at  the  Stress  Focus  (ofc)  is  140  atm.  The  normalized  cone 
resistance  et  the  Stress  Focus  (qcfl)  is  16  as  shown  in  Equation  (5.1)  and  illustrated 
in  Figure  5.4. 


Qcf 1  Ng 


2200 
_14  0 


16 


(5.1) 


The  qcfl  represents  the  bearing  factor  Z0  (i.e.  Nq)  for  confining  stresses  greater  than 
the  Stress  Focus.  Tins  Z0  represents  an  approximate  friction  angle  between  28  and 
31°  using  the  simplistic  bearing  factors  of  Janbu  and  Senneset  (1974)  (sec  Figure  5.5) 
or  using  cavity  expansion  techniques  (using  a  low  rigidity  index)  (Vesic,  1972). 
litis  friction  angle  is  also  within  Tie  range  (i.e.  30°  to  34°)  of  most  granular 
sedimentary  rocks  tested  at  high  stresses  as  was  reported  in  Chapter  2. 


Figure  5.4  is  an  annotated  version  of  Figure  5.3  intended  to  illustrate  all  aspects 
of  the  stress  exponent  and  Stress  Focus  concepts.  As  previously  stated,  a  sand  of  a 
given  relative  density  will  follow,  with  increased  vertical  effective  stress,  a  straight 
path  on  a  logarithmic  plot  to  the  Stress  Focus.  Dense  sands  have  a  low  stress 
exponent,  whereas  loose  sands  have  the  highest  stress  exponent;  the  exponent  may 
approach  1  (one).  Any  point  along  a  given  path  can  be  defined  by  tire  combination 
of  stress  exponent  (c)  and  the  normalized  corse  resistance  (qclJ.  Alternatively,  ans 
point  can  be  defined  by  the  stress  exponent  (c)  and  the  Stress  Focus  point  (that  is 
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Figure  5.2  Trends  of  cone  resistance  versus  veitical  effective  siiess  lo:  T" 
sand  (Baldi,  ct.al  .  198]  i 
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Figure  5.3  Keplctting  of  Baldi,  ct.cl .  1981  data  curves  from  Figure  5.2  in  terrm 
of  Log i „  net  cone  n\ sis ta ace  ve;'<uj-  jLx>g,(i  vertical  effective  .stress 


defined  by  qcf  and  0).  )  as  described  in  Chapter  3. 

Shallow  and  Deep  Bearing  Stress  in  Terms  of  the  Stress  Focus 

The  bearing  stress  conditions  for  shallow  and  deep  conditions  may  now  be  better 
explained  as  a  result  of  the  Stress  Focus  concept  as  illustrated  in  Figure  5.4. 

Shallow  conditions  are  controlled  by  vertical  expansion  based  bearing  capacity 
failure  (defined  by  limit  equilibrium  theories)  and  deep  conditions  are  controlled  by- 
lateral  expansion  based  cavity  expansion  effects.  Shallow  bearing  stress  is  associated 
with  a  linear  relationship  between  cone  resistance  and  vertical  effective  stress  as 
reflected  by  limit  equilibrium  theories  (Terzaghi,  1943,  Durgunoglu  and  Mitchell, 
1975)  and  shown  by  the  transformed  linear  behavior  lines  at  the  top  of  Figure  5.4. 
Deep  bearing  stress  conditions  are  controlled  by  cavity  expansion  (Vcsic,  1972, 
Baligh,  1976)  and  can  be  expressed  exponentially  (i.e.  transform  slope  of  o)  as 
shown  at  the  bottom  of  this  figure. 

The  transition  between  vertical  (limit  equilibrium)  and  horizontal  (cavity 
expansion)  bearing  expansion  is  the  "Critical  depth"  line.  Historically,  this  transition 
■was  expressed  as  the  D/B  ratio  (penetrometer  diameter  to  depth  ratio)  and  is  typically 
10  to  20.  Increasing  relative  density  causes  a  greater  ciuical  depth  transition  as 
shown  in  Figure  5.4  (Kctisel,  1964,  Olsen,  1992)  for  normal  penetrometer  diameters. 
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Figu.e  5.4  Annotated  description  of  the  cone  resistance  Stress  Focus  together  with 
bear  ing  stress  di  kle  l  into  vertical  expression  failure  and  cavny 
expansion  taiiurc 
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Figure  5.5  Correlations  between  cone  resistance  and  friction  angle  based  on  Inn  t 
equilibrium  theories  and  laboratory  chamber  tests 
(Campanclla  and  Robertson,  1983) 
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Using  Chamber  Data  to  Establish  Trends  of 
Stress  Exponent  and  the  Stress  Focus 


Large  diameter  chamber  tests  provide  the  best  experimental  means  for  evaluating 
stress  effects  on  the  cone  resistance,  because  vertical  and  horizontal  stresses  can  be 
controlled  independently.  Numerous  researchers  over  the  last  30  years  have 
performed  CPT  chamber  tests  using  a  variety  of  sands  and  confining  stress 
conditions.  This  section  will  critically  examine  chamber  data  toward  proving  that  the 
cone  resistance  behavior  indeed  has  a  Stress  Focus.  Standard  Penetration  Test  (SPT) 
chamber  tests  results  are  also  examined  in  a  subsequent  section  because  they  include 
data  for  a  large  number  of  well  documented  sand  types.  SPT  chamber  rest  results 
verify  that  the  Stress  Focus  location  is  sand  type  dependent. 

A  computerized  database  of  chamber  tests  was  created  based  on  a  data  listing 
provided  to  the  author  by  Professor  M.  Jamiolkowski  in  1989.  This  listing  represents 
410  chuir.tcv  ’.eels  performed  by  ENEL,  ISMES,  and  NGI  over  a  period  of  12  years 
and  represents  two  sand  types,  namely  Hokksund  sand  and  several  variations  of 
Ticino  sand. 


Evaluation  of  the  Chamber  Data 

The  laboratory  chamber  data  were  initially  divided  according  to  sand  types  and 
further  divided  into  groups  representative  of  normally  consolidated  and 
overconsolidated  conditions.  Each  sand  type  (or  batch  of  a  sand)  was  then  divided 
into  groups  based  on  relative  densif  each  having  an  adequate  number  of  points 
(e.g.  6  minimum)  for  statistical  e  i.uation.  For  each  relative  density  group  (of  a 


6U 


given  sand  type),  the  individual  chamber  test  data  points  were  plotted  as  log]0  (qc)nct 


versus  log,,,  0)  as  shown  in  Figure  5.6  through  Figure  5.9. 


Using  statistical  least  square  fit  correlations  for  each  relative  density  group,  for 
the  purpose  of  determining  the  Stress  Focus,  was  not  successful  because  one  or  two 
poor  data  points  will  shift  the  best  fit  line;  statistical  correlations  are  easily  biased  by 
poor  data  at  or  beyond  the  limits  of  data  ranges  (Taylor,  1990).  Therefore  a  new 
means  was  required  to  determine  Stress  Focus  location.  The  Stress  Focus  location 
can  be  iteratively  varied  (in  the  (qc)nct  and  directions)  until  the  optimum  location 

is  found.  For  each  Stress  Focus  location,  lines  are  projected  from  the  Stress  Focus 
through  each  relative  density  data  group.  The  optimum  Stress  Focus  location  has  the 
best  data  fit  for  all  of  the  relative  density  lines. 

The  Stress  Focus  locations  from  the  preceding  figures  are  summarized  in 
Figure  5.10.  The  main  observation  from  these  plots  is  that  a  cone  resistance  Stress 
Focus  can  be  established  using  CPT  chamber  test  data. 


CFT  Stress  Exponent  Evaluation 

As  shown  in  the  preceding  figures,  cone  resistance  for  each  relative  density  range 
is  associated  with  a  constant  stress  exponent  relationship  (e.g.  straight  line  on  a 
logarithmic  plot)  to  the  Stress  Focus.  The  relationship  of  cone  resistance  stress 
exponent  (c)  to  relative  density  for  the  various  sands  (i.e.  Figure  5.6  to  Figure  5.9)  is 
summarized  in  Figure  5.1 1.  There  is  an  apparent  trend  of  decreasing  stress  exponent 
with  increasing  initial  relative  density. 

Also  shown  in  Figure  5.1 1  is  the  cone  resistance  stress  exponent  of  e=0.61 
suggested  by  the  late  II.  Bolton  Seed  (Seed,  ct  al.,  1981)  for  relative  densities 
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Figure  5.6  Net  Cone  resistance  versus  vertical  effective  stress  for  Hokksand  sand 


Vertical  Effective  Sires?  (atm  units) 


Net  Cone  Resistance  (atm  units) 

1  10  100  1000  10000 


b,  ■  CO' 


74? 

I),  -  917 


0X4 
'O  XX  4 

X  to 

\ 

n  x  4 
□  x 

10  0  — j 


100  0 


1000  0 


Relative  Density  Ranges  (i)r) 
Ticino  sand  (H.itch  number  1) 
(data  from  Jamiolkowski.  190?) 

□  I)r  •  50  lo  02 n 

X  P,  -  7i  to  7L  7 

A  P,  -  90  lo  9“  " 


.Stress  I  oc  u  > 
for  Ticino  sand 
(batch  number  1) 


Figure  5.7  Net  cone  resistance  versus  vertical  effective  stress  for  Ticino  sand 
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Figure  5.8  Net  cone  resistance  versus  vertical  effective  stress  for  Ticino  sand 
(batch  #2) 
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Figure  5.9  Net  cone  resistance  versus  vertical  effective  stress  for  Ticino  sand 
(batch  #4) 
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betwe  .Wo  and  80%.  The  stress  exponent  is  calculated  from  Cp  using:. 


cp  = 


[o'v)c 


(5.2) 


The  Cp-based  stress  exponent  was  established  by  Seed,  et  al.,  (1981)  based  on 
chamber  tests  performed  at  the  University  of  Florida.  These  chamber  tests  were,  at 
the  time,  interpreted  to  have  a  constant  stress  exponent  for  all  relative  densities 
(Schmertmann,  1978a).  The  Cp  based  stress  exponent  of  c=0.61  in  Figure  5.1  i 
envelops  al)  of  the  chamber  based  relative  densities  from  40%  to  80%  and 
specifically  represents  a  relative  density  of  approximately  60%.  The  overfding 
conclusion  from  this  figure  is  that  the  cone  resistance  stress  exponent  is  inversely 
proportional  to  relative  density. 


Overconsolidation  Effects  on  the  Stress  Exponent 

There  is  a  consistent  observation  from  chamber  test  results  (and  field  uniform 
layer  trends  discussed  in  Chapter  6)  that  overconsolidated  soils  have  low  stress 
exponent  values.  Figure  5.12  is  ar  example  of  overconsolidated  chamber  data  having 
a  wide  range  of  relative  density  levels  and  numerous  overconsolidation  levels.  The 
solid  lines  represent  the  normally  consolidated  trend  (from  Figure  5.9)  while  the  data 
points  and  dashed  lines  represent  increasing  overconsolidation  level.  The  stress 
exponents  (c)  for  overconsolidation  range  from  0.06  to  0.24,  which  also  represents 
the  lower  range  for  normally  consolidated  very  dense  sand.  Overconsolidation  at  any 
initial  relative  density  therefore  produces  low  stress  exponents.  A  simple  means  of 
explaining  overconsolidation  effects  is  to  compare  the  conventional  consolidation 
relationship  (e.g.  odometer  test  results)  to  the  cone  resistance  versus  vertical  effective 
stress  relationship  in  Figure  5.13 — both  have  steep  slopes  associated  with 
overconsolidation.  Interestingly,  dense  sands  and  overconsolidated  conditions  are 
both  located  in  the  upper  right  portion  of  the  CPF  soil  characterization  chart 
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Figure  5.11  CPT  stress  exponent  range  for  several  sands  in  terms  of  the  Relative 
density  and  CPT  tone  resistance  stress  exponent 


(Olsen,  l988)(Figure  1.2). 


Chamber  Boundary  Effects  for  Dense  Sands 


Chamber  size  and  chamber  boundary  conditions  can  affect  CPT  measurements. 
Dense  sands  tested  in  laboratory  (constant  stress  boundary)  chambers  produce 
cone  resistances  that  are  lower  than  measured  in  situ  (Bellotti  et  al.,  1979;  Keaveny, 
1985).  Typical  laboratory  chambers  are  too  small  to  accommodate  the  full  cone 
resistance  stress  field  produced  in  dense  sand  by  a  3.6  cm  diameter  probe.  Shear 
stresses  that  would  be  generated  beyond  the  radius  of  a  typical  chamber  do  not 
develop  because  the  pressurized  confining  water  beyond  the  flexible  boundary  cannot 
sustain  shear  stress  or  generate  elevated  lateral  pressure.  This  results  in  lower 
measured  cone  resistance  in  most  chamber  tests  than  are  measured  in  situ  for  dense 
sand.  Typical  chamber  corrcct:on  factors  based  on  the  sand  type,  ratio  oi  cnaiiiuer 
diameter  to  probe  diameter,  and  sand  relative  density  are  shown  in  Figure  5.14.  This 
figure  implies  that  initially  very  loose,  loose  and  medium  dense  sands  arc  not 
influenced  by  boundary  effects  with  typical  chamber  diameters  (i.e.  diameter  ratios 
greater  than  40). 


The  cone  resistance  Stress  Focus  concept  developed  using  chamber  data  in  this 
chapter  is  valid  for  very  loose  to  medium  dense  sands  but  questionable  for  dense 
sands  because  of  chamber  boundary  effects.  Dense  to  very  dense  sand  must  exhibit 
artificially  low  measured  cone  resistances  because  of  chamber  boundary  effects. 
However,  it  appears  that  cone  resistance  behavior  at  al!  relative  densities  (for 
chamber  tests  on  a  single  sand  type)  converge  to  a  common  Stress  Focus  (Figure  5.7 
to  Figure  5.9).  This  convergence  to  a  Stress  Focus  and  the  fact  that  dense  sands 
must  have  lower-than-expected  cone  resistance  (e  g.  Figure  5.14)  must  be  attributable 
to  chamber-tested  dense  sands  having  elevated  stress  exponents  as  shown  in 
Figure  5.15.  Dense  sand  tested  in  chambers  and  in  situ  may  therefore  have  the  same 
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Figure  5.12  Relationship  of  the  cone  resistance  stress  exponent  for 
overconsolidated  sands  in  chamber  tests 
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Figure  5.13  Overconsolidation  effects  illustrated  by  comparing  the  conventional 
consolidation  representation  to  the  observed  relationship  of  cone 
resistance  versus  vertical  effective  stress 
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Figure  5.14  Boundary  effects  on  the  measured  cone  resistance  level  in  laboratory 
chamber  test  in  sand  (Parkin  and  Lunne,  1982) 
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Stress  Focus  but  different  stress  exponents.  Bulging  at  the  chamber  boundary  limits 
the  size  of  the  bearing  pressure  bulb  in  front  of  the  probe  to  that  developed  i.i  a 
lower  relative  density  sand  (i.e.  having  a  higher  stress  exponent).  A  sand  v'ith  an 
elevated  stress  exponent  is  therefore  equivalent  to  a  lowei  relative  density  sand. 

The  elevated  chamber  stress  exponent  of  dense  sands  in  chamber  tests  can  be 
further  illustrated  by  examining  trends  of  stress  exponent  versus  relative  density. 
Stress  exponent  versus  relative  density  trends  from  Figure  5.1 1  having  at  least  3 
relative  density  groups  were  replotted  in  Figure  5.16.  The  elevated  chamber  stress 
exponents  for  dense  sand  are  shown  by  the  vertical  arrows.  This  approximated 
correlation  of  stress  exponent  to  relative  density  is  more  linear  than  ihe  relationship 
derived  from  chamber  results. 
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Figure  5.15  The  reduced  cone  resistance  for  dense  sands  due  to  "chamber  size 
effects"  can  be  represented  with  higher  stress  exponents  because  all 
relative  densities  appear  to  have  a  common  Stress  Focus. 
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Figure  5.16  Possible  modification  of  the  chamber  based  stress  exponent  for  dense 
and  very  dense  sands 
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SPT  Chamber  Evaluation 


This  section  will  evaluate  SPT  chamber  test  data  on  sand  to  show  that  tbe  Stress 
Focus  location  is  soil  type  dependent.  These  SPT  chamber  data  consisted  of  tests  in 
4  well  documented  sand  types  and  were  used  to  study  sand  type  effects  on  location 
of  the  Stress  Focus.  Data  from  SPT  Chamber  tests  were  therefore  examined  for  two 
purposes;  1)  to  show  that  SPT  data  can  be  expressed  with  stress  exponents  and  the 
Stress  Focus,  and  2)  to  show  that  sand  composition  affects  the  Stress  Focus  location. 
If  the  Stress  Focus  is  a  genuine  geotechnical  property  then  SPT  chamber  tests  results 
should  also  reflect  it. 

The  Standard  Penetration  Test  (SPT)  blow  count  is  defined  as  the  number  of 
blows  required  to  advance  a  split-spoon  SPT  soil  sampler  one  foot  into  the  bottom  of 
a  borehole  using  a  standardized  hammer.  The  SPT  sampler  is  resisted  during 
penetration  by  both  static  and  dynamic  end  bearing  and  side  forces  on  the  sampler 
(Schmertmann,  1979a,  1979b,  Douglas,  Olsen,  and  Martin,  1981,  Olsen  1988). 

A  comprehensive  laboratory  study  of  the  SPT  was  performed  using  large 
diameter  chamber  tests  at  the  U.S.  Army  Waterways  Experiment  Station  (WES)  in 
the  mid  1 970's  by  Bieganousky  and  Marcuson  (1976,  1977).  SPT  chamber  tests 
were  performed  in  a  4-foot  diameter  chamber  very  similar  in  concept  to  the  large 
diameter  chambers  used  for  CPT  testing.  For  an  SPT  chamber  test,  a  sand  specimen 
is  prepared,  confining  stress  applied,  a  borehole  drilled  into  the  sand  using  rotary 
wash  techniques  (and  drilling  mud),  and  finally  the  SP  T  sampler  is  driven  into  the 
sand  beginning  at  the  bottom  cf  the  mud-filled  borehole.  SPT  sampler  blow  counts 
were  obtained  from  different  chamber  depths,  along  the  center  of  the  chamber  as 
well  as  radially  out  from  the  chamber  center.  Only  the  center  test  values  were  used 
for  the  current  evaluation. 


SPT  Chamber  Data  Plotting 


The  procedures  for  evaluating  WES  SPT  chamber  data  are  the  same  as  used  for 
CPT  cone  resistance  data  in  the  last  section.  For  Reid  Bedford  sand,  trends  related 
to  individual  relative  density  groups  are  shown  in  Figure  5.17  to  Figure  5.19;  all 
relative  density  data  and  trends  are  combined  in  Figure  5.20.  For  Standard  Concrete 
sand,  all  relative  density  group  trends  are  shown  in  Figure  5.21,  for  Platte  River 
sand,  all  of  the  relative  density  groups  trends  are  shown  in  Figure  5.22,  and  finally 
for  Ottawa  sand,  only  ore  relative  density  group  is  available  as  shown  in  Figure  5.23 
(all  data  from  Bieganousky  and  Mr  mson,  1976,  1977).  The  correlations  by 
Gibbs  and  Holtz  (1957)  are  shown  i  Figure  5.24  for  comparison. 


SPT  Stress  Focus  and  Sand  Type 

The  SPT  Stress  Focus  information  from  the  preceding  figures  are  summarized  in 
Figure  5.25  together  with  trends  of  mean  grain  size  and  material  composition  in 
terms  of  feldspar  percentage.  Quartz  is  the  primary  material  type  of  these  sands 
with  feldspar  being  the  secondary  material  type — increasing  feldspar  content  may 
reflect  a  lower  overall  sand  strength.  Compared  to  quartz,  feldspar  is  softer,  has 
lower  compressive  strength,  and  has  distinct  cleavage.  In  Figure  5.25,  as  feldspar 
content  increases,  the  depth  for  the  Stress  Focus  inert  ases.  A  deeper  Stress  focus 
(due  to  feldspar  content  increase)  has  a  lower  calculated  normalized  cone  resistance 
at  the  Stress  Focus  (qcf])  (see  page  54  for  definition)  which  means  that  the  friction 
angle  at  high  confining  stress  is  also  lower.  Chapter  6  will  also  show  that  the  cone 
resistance  Stress  Focus  is  soil  type  dependent  and  exhibit  the  same  general  trend  as 
the  sand  mean  grain  size  in  Figure  5.25. 
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Figure  5.17  Cone  resistance  chamber  data  for  a  relative  density  of  33%  for  Reid 
Bedford  sand 
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Figure  5.18  Cone  resistance  chamber  data  tor  a  relative  density  of  45%  for  Reid 
Bedford  sand 
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Figure  5.19  Cone  resistance  chamber  data  for  a  relative  density  of  67%  for  Reid 
Bedford  sand 
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Figure  5.20  Cone  resistance  chamber  data  for  all  relative  densities  for  Reid 
Bedford  sand  together  with  the  Stress  Focus 
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Cone  resistance  chamber  data  for  Platte  River  sand  together  with  the 
corresponding  Stress  Focus 
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Cone  resistance  chamber  data  trends  for  Gibbs  &  Holtz  data 
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Figure  5.25  Compandor  of  SPT  Stress  Focus  locations  for  different  sands  to  the 
resulting  trends  of  Feldspar  percentage  and  Ds„ 


Chamber  Test  Stress  Exponent  Discussion 


The  SPT  stress  exponents  calculated  for  data  at  each  relative  density  (from  the 
proceeding  figures  on  SPT  chamber  test  evaluation)  are  summarized  in  Figure  5.26. 
These  SPT  stress  exponents  confirm  the  general  trends  based  on  normalized  SPT 
stress  exponent  presented  by  Seed  et  al.  (1981)  that  are  also  shown  in  this  figure. 
However,  it  is  shown  in  Chapter  7  that  the  SPT  stress  exponent  should  be  higher 
because  of  mud  pressure  effects  associated  with  SPT  chamber  testing. 

The  stress  exponents  developed  from  SPT  and  CPT  chamber  evaluations  are 
shown  in  Figure  5.27  and  indicate  a  difference  of  approximately  0.18.  This  figure 
will  also  be  referenced  in  Chapter  6  during  the  introduction  of  the  new  SPT 
normalization  technique. 

It  is  also  important  to  note  that  the  distribution  of  SPT  stress  exponent  for  a 
given  relative  density  in  Figure  5.26  may  not  be  due  entirely  to  data  scatter.  The 
stress  exponent  can  be  related  to  sand  type  composition  as  shown  in  Figure  5.28  of 
which  is  plotted  in  terms  of  feldspar  content  at  a  relative  density  of  36%  (likely  a 
contractive  condition).  As  previously  stated  on  page  77  and  shown  in  Figure  5.25, 
increasing  feldspar  content  increases  the  Stress  Focus  vertical  effective  stress. 
Increasing  feldspar  content  also  appears  decrease  the  stress  exponent  as  shown  in 
Figure  5.28.  Therefore,  conceptually,  the  combined  effects  of  increasing  vertical 
effective  stress  and  decreasing  stress  exponent  can  be  illustrated  in  Figure  5.29.  The 
Stress  Focus  location  and  stress  exponents  therefore  appear  to  be  interrelated. 
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Figure  5.26  Relationship  between  relative  density  and  SPT  stress  exponent  for 
different  sands  from  SPT  chamber  tests 
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Figure  5.28  Relationship  between  SPT  stress  exponent  and  Feldspar  percentage  at 
a  relative  density  of  36%  for  chamber  tests 
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interrelation  of  the  Stress  Focus  and  the  stress  exponent  (for  a  very 
loose  sand  condition)  based  on  evaluation  of  CPI'  cone  resistance 
chamber  data  fc  increasing  feldspar  content 


Conclusion 


The  CPT  cone  resistance  can  be  represented  with  the  Stress  Focus  concept.  The 
cone  resistance  versus  confining  stress  relationship  for  different  relative  densities  will 
converge  with  increasing  vertical  effective  stress  to  a  Stress  Focus  tnat  is  dependent 
on  sand  type.  This  sand  type  dependent  Stress  Focus  location  appears  to  be  based  on 
gradation  effects  and  material  composition.  The  Stress  Focus  location  appears  to 
occur  at  a  overburden  stress  ranging  from  70  to  300  Atm.  Dense  sands  have  the 
lowest  stress  exponent  and  loose  sands  have  a  stress  exponent  of  approximately  one. 
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Chapter  6 


CPT  Prediction  of 
CPT  Normalization  Parameters 
Developed  using 
Uniform  Soil  Layer  Data 


Introduction 

Stress  normalization  provides  the  means  to  account  for  confining  stress  influence 
on  in  situ  measurements  and  geotechnical  properties.  However,  to  be  useful  the 
stress  normalization  parameters  must  be  predictable  using  tire  field  CPT 
measurements  without  auxiliary  information.  Data  from  CPT  tests  in  uniform  soil 
layers  can  be  used  to  establish  correlations  for  prediction  of  the  CPT  cone  resistance 
stress  exponent.  Specifically,  field  CPT  daia  from  tests  in  uniform  soil  layers  can  be 
used  for  several  purposes:  to  show  that  the  Stress  Focus  exists,  the  Stress  Focus 
location  (i.e.  equivalent  overburden  stress  and  cone  resistance)  is  dependent  on  soil 
type,  the  cone  resistance  stress  exponent  can  be  estimated  using  field  CPT  data,  and 
to  prove  that  the  CP-"  soil  characterization  chart  concepts  arc  valid. 
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Limitations  and  Merits  of  Using 
In  Situ  Data  for  Stress  Normalization 


Evaluation  of  uniform  soil  layers  is  more  difficult  than  evaluation  of  chamber  test 
data.  The  main  difference  between  in  situ  uniform  layer  CPT  data  and  chamber  test 
data  is  that  in  situ  data  can  represent  all  soil  types,  whereas,  to  date  chamber  data 
can  only  represent  clean  sands  (and  in  some  cases  sand  with  fines)  of  low  to  medium 
relative  density.  On  the  other  hand,  the  initial  relative  density  and  confining  stress 
levels  can  accurately  be  varied  in  chamber  test  programs.  The  principal  merit  of 
in  situ  uniform  soil  layer  data  is  that  it  can  represent  the  full  spectrum  of  soil  types 
and  thus  that  it  can  be  used  to  determine  the  Stress  Foci  for  all  soil  types.  However, 
when  a  very  uniform  soil  layer  is  found,  it  represents  only  one  relative  strength  level; 
e.g.  a  sand  with  a  70%  relative  density  or  a  clay  with  an  undrained  normally 
consolidated  strength  ratio  (c/p)  of  0.32. 

A  single  uniform  soil  layer  has  a  constant  trend  of  log  net  cone  resistance  versus 
log  ov  that  can  be  represented  with  a  constant  exponent  slope.  To  establish  the 

Stress  Focus  for  a  single  soil  type  requires  numerous  uniform  layers  of  varying 
relative  strength  consistency  levels.  A  database  representing  uniform  soil  layers  must 
therefore  be  large  enough  for  the  purpose  of  establishing  the  Stress  Focus  for 
numerous  soil  types. 
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Establishing  In  Situ  Uniform  Layer  Trends 


A  uniform  soil  layer  is  defined  as  a  constant  soil  type  at  a  constant  initial  relative 

c 

strength  (e.g.  friction  angle  for  sand  or  —  for  clay).  A  plot  of  Log  (qc)nct  versus 

P  . 

Log  (ov)  (and  Log  (fs)  versus  Log  (0^) )  will  follow  a  straight  path  over  the  uniform 

soil  layer  interval  as  illustrated  in  Figure  6.1.  These  net  cone  resistance  relationships 
can  be  projected  to  the  normalizing  stress  of  ov=l  atm  to  determine  the  normalized 

cone  resistance,  qclc.  The  log-log  slope  is  the  cone  resistance  stress  exponent  (c). 

Soil  samples  are  not  required  for  this  technique,  although  proper  soil  classification  is 
useful  for  indexing  purposes. 

Logarithmic  plotting  emphasizes  the  lower  stresses,  therefore  special  care  was 
exercised  while  identifying  uniform  soil  layers — For  example,  a  depth  interval  of 

3  to  8  ft  has  approximately  the  same  log  ( 0'v)  differential  (i.e.  0.3)  as  the  depth  range 

from  30  to  60  ft.  Shallow  uniform  soil  depth  zones  are  artificially  emphasized.  It 
would  be  easy  to  use  a  simple  Alog10  (0'v)  criterion  such  as  0.2  to  represent  a 

minimum  uniform  depth  zone.  However,  it  was  important  to  examine  deep  layers 
for  small  ALog]0  ( cjv)»  such  as  0.08,  that  could  represent  a  uniform  12  foot  layer  at 

o  =  2  atm. 


Approximately  600  field  CPT  soundings  were  examined  for  depth  zones  having  a 
constant  relative  strength.  This  was  accomplished  by  plotting  all  CPT  soundings 

using  the  Log10  (net  qc)  versus  Logi0  (0v)  and  then  examining  the  traces  for  any 
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Figure  6.1  Determining  the  cone  resistance  normalization  value  (qelc)  and 
corresponding  stress  exponent  (c)  using  uniform  soil  layers 
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depth  zone  having  constant  slope  (indicating  a  uniform  relative  strength  level). 
Identifying  uniform  soil  layers  using  CPT  data  was  a  developed  skill  that  improved 
as  more  data  was  evaluated.  Only  78  uniform  soil  layers  were  extractable  from  the 
600  CPT  soundings.  An  example  of  a  uniform  soil  layer  is  shown  in  Figure  6.2  with 
the  soil  layer  limits  and  normalized  parameters  shown  by  computer-plotted  lines 
based  on  the  data  from  the  database.  For  each  uniform  soil  layer,  the  following 
information  was  extracted; 

1)  normalized  cone  resistance,  qcle  (e.g.  equivalent  (qc)net  at  0'  =  1  atm) 

2)  qc  log-log  stress  exponent  (c) 

3)  Normalized  sleeve  friction  resistance,  fsle 

(e.g.  the  equivalent  fs  at  0v=  1  atm) 

4)  fs  log-log  stress  exponent  (s) 

5)  Calculated  normalized  friction  ratio,  Rrle,  based  on  fsle  and  qclc 

6t  Soil  tvr>P.  if  thero  is  a  noarbv  Viorino  or  tVip  Hpnocit  tvrv>  if  tnrmm 

/  V  J  *  - —  — - -  J  ~  — ‘O  V  f”  **  “ 

7)  The  Academic  Quality  Index  (AQI)  which  is  based  on  soil  layer  uniformity 

potential  and  quality  of  the  cone  and  sleeve  fiction  resistance  measurements -(see 

Chapter  8  for  description).  The  AQI  is  a  subjective  quality  index  based  on  the 

academic  grading  scale;  75%  is  average,  85%  is  good  and  95%  is  excellent. 

8)  Top  and  bottom  ov  limits  of  the  soil  layer. 


Establishing  Cone  Resistance  Stress  Exponents 
using  the  CPT  Soil  Characterization  Chart 


The  CPT  soil  characterization  chart  was  used  to  develop  predictive  contours  for 
the  cone  resistance  stress  exponents  using  data  from  uniform  soil  layers.  As 
illustrated  in  Figure  6.3,  unifonn  soil  layers  are  identified  (Step  ti\),  normalized  CPI 
parameters  arc  calculated  (Step  til  &  #3),  the  results  arc  assigned  to  a  point  on  the 
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Vertical  Effective  Stress  (atm  units) 


Cone  Resistance  -  Total  Stress 
(atm  units) 


Figure  6.2  Example  of  in  situ  cone  resistance  data  having  a  uniform  soil  layer 
with  the  stress  exponent,  c,  and  normalized  cone  resistance,  q  , 


CPT  soil  characterization  chart  (Step  #4),  and  the  point  is  assigned  die  cone 
resistance  stress  exponent  (c)  (Step  #5).  This  point  represents  one  uniform  soil  layer. 
Figure  6.4  to  Figure  6.13  represent  a  few  interpreted  uniform  soil  layers  for  different 
soil  types  and  relative  strengths.  Also  shown  at  the  bottom  left  comer  of  each  figure 
is  a  representation  of  the  CPT  soil  classification  chart  (logarithmic  3  cycles  by  2 
cycles)  and  the  normalized  point  that  represents  the  uniform  soil  layer.  All  available 
data  from  uniform  soil  layers  are  plotted  in  Figure  6.14  using  the  techniques 
illustrated  in  Figure  6.3.  The  estimated  contours  of  cone  resistance  stress  exponent 
based  on  the  trends  of  plotted  data  points  are  also  shown  in  Figure  6.14. 


Discussion  of  Cone  Resistance  Stress  Exponent  on  the  CPT  Soil 
Characterization  Chart 

Figure  6.15  is  a  replot  of  Figure  6.14  with  the  soil  type  of  each  undorm  soil 
layer  data  point  indicated.  Also  shown  in  this  figure  are  CPT  soil  classification  lines 
from  Olsen  (1988).  lire  soil  classifications  of  the  uniform  soil  layers  approximately 
match  die  CPT  soil  classifications.  Figure  6.16  is  an  annotated  version  of 
Figure  6.14  that  also  describe  various  soil  characterization  zones.  Note  the  lack  of 
uniform  soil  data  for  soil  mixtures  and  silt  in  the  middle  of  Figure  6.14;  Thick 
layers  of  a  soil  mixture  or  silt  are  seldom  found  in  nature  because  the  required 
deposition  dynamics  are  difficult  to  maintain  over  large  soil  depths. 

The  contours  of  cone  resistance  stress  exponent  in  Figure  6.16  reveal  important 
information  about  the  estimated  cone  resistance  contours.  Loose  sands  typically  have 
low  friction  ratios  (Masood,  1990)  because  tire  Nq  bearing  factor  is  low.  Loose 
sands  should  also  exhibit  little  Molir  envelope  curvature  and,  therefore  have  a  stress 
exponent  very  close  to  one  as  described  in  Chapter  2.  The  stress  exponent  contour 
of  one  approaches  the  annotated  loose  sand  zone  on  the  CPT  soil  characterization 
chart,  as  expected.  Also,  normally  consolidated  clays  should  have  stress  exponents 
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Figure  6.4 
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Uniform  sand  layer  at  the  Holmen,  Norway 
(NGI  research  site) 
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Figure  6,8  Uniform  silty  sand  layer  in  the  Aitchafalaya  River  (dredging  study 
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Figure  6.9  Undo; m  silly  clay  hiycr  in  the  Au-l/iTalaya  Rjver  (Dredging  study) 
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Uniform  .sand  iaycr  segment  at  the  McDondalds  Farm,  Sea  Island, 
Canada 
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Vert ica 1 


Normalized  Cone  Resistance  (atm  units) 


Cone  resistance  stress  exponent  (c)  contours 
using  data  from  uniform  soil  layers 
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Figure  6.14 


Establishing  contours  of  cone  resistance  stress  exponent  using  CPT 
data  from  in  situ  uniform  soil  layers  with  constant  relative  strength 
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Normalized  Cone  Resistance  (at.m  units) 


General  soil  type  from  uniform  soil  layers 
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Figure  6.15  Comparing  soil  types  from  uniform  soil  layers  to  the  CPT  soil 
classification  line  (from  Olsen,  1988) 
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Normalized  Cone  Resistance  (aim  units) 


Figure  6.16  Annotated  version  of  the  cone  resistance  stress  exponent  chart  showing 
soil  zones  (and  strength  character)  together  with  CPT  soil 
characterization  line 
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of  one  or  slightly  less;  the  stress  exponent  contour  for  normally  consolidated  clay  is 
between  0.8  and  1,  as  expected.  Dense  and  overconsolidated  sands  have  stress 
exponents  contours  between  0.15  and  0.55,  as  expected.  Overconsolidated  clays  are 
shown  to  have  stress  exponent  contours  between  0.55  and  0.75,  also  as  expected. 
Unstable  silty  clay  (i.e,  silt  structure  with  underconsolidated  clay  matrix)  have  stress 
exponents  of  about  1.1  to  1.5  and  are  located  at  the  bottom  right  of  the  CPT  soil 
characterization  chart,  which  is  a  new  finding. 


The  cone  resistance  stress  exponent  contours  displayed  in  Figure  6.16  are 
approximately  perpendicular  to  the  CPT  soil  characterization  (classification)  lines. 

For  a  given  soil  type,  increasing  the  relative  strength  (e.g.  relative  density  for  sand) 
should  move  the  CPT-based  soil  characterization  (i.e.  qcle  and  FRle  point  on  the 
chart)  along  a  contour  of  constant  soil  type  to  the  upper  right  (Olsen,  1988,  1984).  It 
was  shown  in  Chapter  2  that  as  relative  strength  increases,  the  stress  exponent 
decreases.  Soil  type  and  relative  strength  should  be  independent  of  each  other 
because  they  define  different  aspects  of  a  soil;  i.e.  what  it  is  and  the  slate  that  it's  in. 
For  a  given  soil  type,  a  relative  strength  increase  can  be  represented  on  the  CPT  soil 
characterization  chart  as  a  decreasing  stress  exponent  contour  level  (i.e.  increasing 
relative  strength)  along  a  given  soil  type  contour  as  shown  in  Figure  6.16.  The  CPT 
soil  characterization  chart  is  therefore,  verified  in  part,  therefore,  because  the  soil 
classification  contours  are  approximately  perpendicular  to  the  stress  exponent 


contours. 
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Establishment  of  the 

Sleeve  Friction  Resistance  Stress  Exponent 


The  CPT  sleeve  friction  resistance  is  obtained  through  a  high-strain  strength 
measurement  where  the  confining  stress  on  the  sleeve  is  equal  to  the  cylindrical 
cavity  expansion  pressure  (Masood,  1 990).  Sleeve  friction  resistance  is  also 
influenced  by  many  of  the  same  geotechnical  properties  as  the  cone  resistance. 
However,  the  sleeve  friction  is  a  high  strain  strength  measurement;  whereas  the  cone 
resistance  is  a  bearing  stress  (many  times  higher  than  confining  stress)  that  is 
dependent  on  both  small  and  large  strains.  Therefore,  the  sleeve  friction  stress 
exponent  should  not  be  equal  to  the  cone  resistance  stress  exponent. 


Direct  Correlation  of  Sleeve  and  Cone  Stress  Exponents 

The  first  attempt  toward  developing  a  correlation  between  sleeve  friction  stress 
exponent  and  cone  resistance  stress  exponent  using  uniform  soil  layer  data  was 
simply  a  scatter  plot  of  sleeve  resistance  stress  exponent  (s)  versus  cone  resistance 
stress  exponent  (c)  as  shown  in  Figure  6.17.  This  figure  shows  a  large  data  scatter 
but  suggests  that  both  stress  exponents  are  approximately  equal  without  considering 
soil  type  or  relative  strength  consistency  contributions.  The  large  scatter  does 
indicate  that  other  factors  are  influencing  the  relationship  between  the  two  exponents. 

The  sleeve  friction  is  more  difficult  to  measure  accurately  than  the  cone 
resistance  because  of  mechanical  constraints  such  as  strain-gauge  accuracy,  thermal 
effects,  and  dirt  clogging  the  joints  (Olsen,  1994).  As  a  result,  it  was  not  possible  to 
establish  contours  of  the  sleeve  friction  stress  exponent  (s)  or  the  ratio  of  s/c  on  the 
CPT  soil  characterization  chart.  A  reliable  s/c  ratio  would  allow  the  less  accurate 


115 


Cone  resistance  stress  exponent  (c) 


Correlation  of  the 

sleeve  friction  resistance  stress  exponent  (s) 
to  the 

cone  resistance  stress  exponent  (c) 

AQI  greater  than  88~  (high  quality  data) 
Q  AQI  between  73",  and 
CP  AQ!  less  than  73%  (lower  quality  data) 


"1 

i 

1 


1  -|  - 


□ 


ffi 


■  / 

/ 


tb 


■  i  '  :  ‘  i  i  t 


'"ini't  i  "i  ii  i  i  i  i 


Sleeve  friction  resistance  stress  exponent  (s) 


Figure  6,17  Relationship  of  sleeve  friction  resistance  stress  exponent  (s)  to  the 
cone  resistance  stress  exponent  (c)  using  uniform  soil  layer  data 
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sleeve  resistance  stress  exponent  (s)  to  be  indexed  to  the  more  accurately  determined 
cone  resistance  exponent  (c).  The  only  correlation  that  proved  satisfactory  was  that 
between  the  s/c  ratio  and  normalized  sleeve  friction  resistance,  as  shown  in 
Figure  6.18.  The  s/c  ratio  is  interesting  because  it  is  analogous  to  the  ratio  of 
cylindrical  expansion  pressure  (on  the  sleeve)  to  the  spherical  plus  cylindrical 
expansion  pressure  (for  the  cone)  as  illustrated  by  Equation  (6.1).  The  s/c  appears 
to  correlate  to  the  normalized  sleeve  friction  resistance  quite  well  probably  because 
the  sleeve  measurement  is  primarily  influenced  by  cylindrical  cavity  effects. 

sleeve  stress  exponent  (s)  x  _  cylindrical  pressure _  (6  1) 

cone  resistance  exponent  (c)  cylindrical  pressure  effects  +  spherical  pressure 
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Figure  6.18 


s/c  Relationship 
(AQI  are  quality  indices) 

4*  AQls  from  90%  to  100%  (Excellent  data) 
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Q  AQls  from  0%  to  100%  (AH  data) 


Correlation  of  the  s/e  ratio  (sleeve  to  cone  stress  exponents)  to  the 
normalized  sleeve  friction  resistance  (fsle) 
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Determining  the  Soil  Type  Dependent  Stress  Focus 
using  Uniform  Soil  Layer  Data 


CPT  data  from  tests  in  uniform  soil  layers  was  used  to  show  that  the  Stress 
Focus  location  is  soil  type  dependent.  The  first  step  in  this  process  involved  dividing 
the  in  situ  uniform  soil  layer  database  into  3  soil  classification  groups:  clays,  soil 
mixtures,  and  sands.  This  grouping  of  data  was  accomplished  based  primarily  on 
soil  type  observed  in  nearby  borings  and  secondarily  on  the  information  about  the 
site  geology.  Only  a  small  portion,  of  the  database  was  composed  of  uniform  soil 
mixtures  due  to  the  geologic  scarcity  of  such  deposits.  These  soil  classification 
groups  were  also  further  distinguished  as  either  normally  consolidated  or 
overconsolidated  soils.  Overconsolidated  soils  have  very  low  stress  exponents 
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of  +hc  moderately  overconsolidated  sands  probably  were  classified  as  normally 
consolidated  dense  sands  because  their  behavior  is  similar  to  that  of  a  dense  sand. 


Developing  the  Stress  Focus  Relationship  for  Different  Soil  Types 

All  uniform  clay  layer  data  collected  in  this  research  program  are  shown  in 
Figure  6.19.  The  projections  for  normally  consolidated  uniform  clay  layers  (using 
different  clay  types)  produces  a  Stress  Focus  at  an  overburden  stress  of 
approximately  9  atm. 

A  clay  Stress  Focus  can  also  be  observed  using  Figure  6.20  by  plotting  the  clay 
and  silty  clay  trends  using  the  stress  normalization  format  from  Chapter  3  with  the 
resultant  shown  in  Figure  6.21.  The  clay  trends  in  Figure  6.21  appear  to  converge  to 
a  Stress  Focus  at  an  overburden  stress  of  approximately  5  to  10  atm,  that  is 


119 


ertioa!  Effective  Stress  (atm  unit 


1 


1  000 


Net  Cone  Resistance  (atm  units) 
10  100 


zjunjform  clay  layers 
3  U  0 .  U  — 1 - — 


Figure  6.19  Projections  of  trends  for  uniform  in  situ  day  layers  in  terms  of  log  net 
cone  resistance  versus  log  vertical  effective  stress 


approximately  the  same  as  9  atm  from  Figure  6.19. 


The  data  in  Figure  6.19  indicate  that  normally  consolidated  clays  have  a  stress 
exponent  (c)  of  about  0.8  to  1.0  and  a  normalized  cone  resistance  (qc]e)  of  about 
3.5  atm.  A  cone  resistance  stress  exponent  of  0.8  (also  shown  in  Figure  6.16)  is 
significant  because  prior  geotechnical  property  formulations  (e.g.  prediction  of  Su  for 
clay)  have  historically  implied  a  stress  exponent  of  one.  The  stress  exponent  of  0.8 
to  1  results  because  net  cone  resistance  rather  then  the  raw  cone  resistance  was  used 
to  establish  the  trend — soil  strength  level  is  more  proportional  to  the 
net  cone  resistance  than  the  raw  cone  resistance  (Chapter  4). 

The  Stress  Focus  for  clay  occurs  at  a  vertical  effective  stress  of  approximately 
9  atm.  The  undrained  strength  at  the  clay  Stress  Focus,  assuming  a  c/p  of  0.31,  is 
approximately  5600  psf  (280  KPa).  A  clay  strength  of  4000  psf  (200  KPa)  has  a 
strength  discriptor  of  hard  (Peck,  el  al.,  1974)  and  is  at  the  lower  boundary  for 
classification  as  a  compaction  sha'e.  It  is  therefore  reasonable  to  assume  that  a  shale 
like  behavior  starts  at  the  Stress  Focus  for  clay  deposition. 

Overconsolidated  clays  were  observed  to  have  a  higher  qc)r  and  lower  stress 
exponents  than  for  a  normally  consolidation  condition.  The  overconsolidated  stress 
exponents  typically  ranged  from  0.4  to  0.6,  whereas  normally  consolidated  clay  stress 
exponents  range  from  0.8  to  1.  The  clay  stress  exponent  due  to  overconsolidation  is 
therefore  approximately  50%  to  60%  of  that  associated  with  the  normally 
consolidated  condition. 

All  uniform  sand  layer  projections  are  shown  in  Figure  6.22.  The  uniform  sand 
layer  trends  are  shown  to  project  to  a  general  Stress  Focus  because  different  sand 
types  can  not  be  distinguished.  Overconsolidated  sands  have  low  observed  stress 
exponents,  as  expected. 


121 


Void  ratio 


Effective  Vertical  Stress  (at 


1.0 


10 

100 

1000 


6.2  i  Replot  of  selected  soil  types  from  Figure  6.20  in  using  stress 

normalization  format  to  show  that  clays  have  an  approximate  Str 
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Figure  6.22  Projections  of  trends  for  uniform  sand  layers  in  terms  of  log  net  cone 
resistance  versus  log  vertical  effective  stress 
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All  uniform  soil  mixture  layer  projections  are  shown  in  Figure  6.23.  Soil 
mixtures  arc  more  complex  and  cover  a  wider  range  of  conditions  than  do  clays  or 
sands.  As  a  result,  the  soil  mixture  layers  projections  are  more  complicated  than  for 
clay  or  sand.  For  example,  dense  silty  sand  can  be  stronger  than  a  dense  clean  sand 
and  a  honeycombed  silty  clay  may  have  a  moderately  high  sensitivity  much  like 
salt-leached  clay.  The  locus  of  Stress  Focus  trends  in  Figure  6.23  represents  the 
zone  projection  for  all  soil  mixture  classifications.  Clayey  soil  mixtures  project  to 
the  upper  left  of  the  elliptical  locus  and  stronger  sandy  silts  with  higher  qcle  are 
project  to  the  lower  light  of  the  elliptical  locus.  Unstable  soil  mixtures  are  also 
depicted  and  project  into  the  soil  mixture  Stress  Focus  to  the  upper  left  of  the 
elliptical  locus. 


Soil  Type  Dependent  Trend  of  the  Stress  Focus 

There  is  enough  information  at  this  point  to  support  establishment  of  a 
soil  type-dependent  Stress  Focus  that  will  be  called  the  "Stress  Focus  boundary". 
Shown  in  Figure  6.24  are  the  Stress  Foci  for  clay  (from  Figure  6.19),  son  mixtures 
(from  Figure  6.23),  and  sands  (from  Figure  6.22).  Also  shown  are  the  Stress  focus 
locations  for  Flokksund  and  Ticino  sands  from  laboratory  chamber  test  evaluations 
(from  Chapter  5).  The  soil  type  dependent  Stress  Focus  boundary  is  shown  as  a 
thick  dashed  line  that  passes  through  all  of  the  Stress  Foci.  This  boundary  is  curved 
and  represents  an  increasing  normalized  Stress  Focus  (qun)  (see  Chapter  5  for 
definition)  as  the  soil  classification  changes  from  clay  to  sand.  Conceptually,  the 
Stress  Focus  boundary  represents  the  boundary  between  a  soil-like  and  (sedimentary) 
rock-like  mechanical  behavior. 

'Ihe  soil  type  dependent  Stress  l  ocus  boundary  is  depick  d  again  in  Figure  6.25 
to  show  behavior  at  confining  stresses  Lower  and  higher  than  th'  Stress  Focus.  Fo, 
confining  stress  less  than  the  Stress  Locus,  a  few  example  soil  typos  arc  depicted 
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Figure  6.23  l'i  ejections  of  trends  fur  uniform  soil  mixture  layers  in  terms  of  Log 
cone  resistance  versus  Log  vertical  effective  stress 


together  with  different  relative  strength  trends.  Also  shown  is  the  probable  behavior 
for  all  soil  types  at  confining  stresses  beyond  the  Stress  Focus;  this  behavior  is  likely 
at  a  stress  exponent  of  one  and  in  a  domain  where  strength  behavior  is 
(sedimentary)  rock-like. 


Conclusions 

The  main  requirement  for  normalizing  the  CPT  measurements  is  the  stress 
exponent.  Figure  6.16  can  be  used  to  determine  the  cone  resistance  stress  exponent 
using  field  CPT  data.  The  sleeve  fiiction  resistance  stress  exponent  is  estimated 
using  Figure  6.18  in  conjunction  with  Figure  6.16.  Cone  resistance  stress  exponent 
contours  are  approximately  perpendicular  to  soil  classification  lines  in  the  CPT  soil 
characterization  cnart;  representing  indirect  support  for  the  validity  of  the  chart. 
Finally,  in  situ  uniform  soil  layer  data  and  chamber  test  data  were  used  to  establish 
the  soil  type  dependent  Stress  Focus  boundary. 
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The  soil  type  dependency  irend  for  the  Stress  Focus  based  on  uniform 
soil  layer  and  laboratory  chamber  test  evaluation 
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Chapter  7 


Normalization  of 
Selected  Geotechnical  Properties 


Introduction 

Geotechnical  property  normalization  is  vital  toward  development  of  CPT-based 
prediction  correlations  and  prediction  of  in  situ  geotechnical  properties.  CPT 
prediction  of  geotechnical  properties  requires  accurate  stress  normalization; 
otherwise,  for  very  shallow  or  very  deep  conditions,  the  stress  normalization  will 
itself  induce  errors  into  the  predictive  process.  This  chapter  will  describe  new 
normalization  techniques  for  the  Standard  Penetration  Test  (SPT)  blow  count  and  the 
shear  wave  velocity. 


Standard  Penetration  Test  (SPT)  Normalization 

Normalization  is  required  to  convert  the  measured  SPT  blow  count,  N,  to  a 
representative  value  that  would  be  measured  when  vertical  effective  stress 
equals  1  tsf  (e.g.  approximately  1  atm).  This  normalized  SPT  blow  count  can  then 
be  used  to  predict  a  variety  of  normalized  geotechnical  properties,  such  as 
liquefaction  resistance  and  friction  angle.  A  new  observation  will  be  introduced  that 
results  from  the  fact  that  SPT  laboratory  chamber  tests  arc  subject  to  a  stress  relief 
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associated  with  the  constant  mud  pressure  in  the  chamber  borehole.  During  a  SPT 
laboratory  chamber  test,  the  borehole  mud  height  is  always  the  same  (approximately 
6  feet),  regardless  of  the  chamber  confining  stress  level.  Quantifying  this  pressure 
relief  will  show  that  the  SPT  chamber-derived  stress  exponents  are  too  low.  It  will 
also  be  shown  that  the  correct  SPT-based  stress  exponent  for  field  applications  is 
equal  to  the  CPT-based  stress  exponent. 


Historical  SPT  Blow  Count  Normalization 

The  state-of-the-practice  for  Standard  Penetration  Test  (SPT)  blowcount 
normalization  for  the  last  12  years  was  developed  by  Seed  et  al.  (1983)  and  later 
confirmed  by  Skempton  (1986).  This  SPT  normalization  technique  uses  an  exponent 
of  the  vertical  effective  stress  for  two  relative  density  ranges.  Seed  et  al.  (1983) 
used  data  from  the  Bieganousky  and  Marcuson  (1976,  1977)  study  to  develop  the 
SPT  CN  normalization  parameter  shown  in  Figure  7.1.  The  SPT  blow  count,  N60, 
(i.e.,  measured  using  equipment  that  delivers  60%  of  theoretical  maximum  free-fall 
hammer  impact  energy  or  adjusted  to  simulate  same)  is  converted  to  a 
stress-normalized  value,  (N,)60,  at  an  equivalent  vertical  effective  stress  of  1  ton/ft2 
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using  the  following  equation: 


(■^1^60  ^60 


60 


(7.1) 


where 
(Nj  )6o 

N6o 


Ov 


n 


CN 


Normalized  SPT  blow  count,  equivalent  SPT  at  a  vertical 
effective  stress  of  1  ton/ft2  (approximate  atmospheric  pressure) 
SPT  measured  blow  count  at  60%  of  theoretical  maximum  free- 
fall  energy  which  is  the  US  approximate  average  achieved  in 
practice 

Vertical  effective  stress  in  tons/ft2  (e.g.  approximately 

1  atm,  100  KPa,  0.1  MPa,  etc.,) 

SPT  overburden  stress  exponent 

( n-'-'H  ^  fnr  ^pneitipe  frnm  40  to  AO®/.  onH 

AVI  tviutll  V  UVUkllUVD  A1  Will  V  \J  VV  V  V  /  U  U11M 

n~0.45  for  relative  densities  from  60  to  80%) 

SPT  normalization  factor  (See  Figure  7.1) 


A  range  of  CN  values  back-calculated  from  field  and  laboratory  chamber  test 
results  is  shown  in  Figure  7.2  (Skempton,  1986)  in  terms  of  overburden  stress. 
Overburden  stress  is  always  represented  in  terms  of  tons/ft2  (approximately  1  atm) 
and  CN  is  equal  to  1  when  the  vertical  effective  stress  is  1  tsf  (i.e.  approximately  1 
atm).  For  the  last  12  years,  equivalent  SPT  stress  exponents  of  0.45  for  relative 
densities  from  60  to  80%  and  0.55  for  relative  densities  between  40  to  60%  as  shown 
Figure  7.1  have  been  widely  used.  The  CPT  chamber-based  stress  exponents  are 
typically  0.6  to  0.7  (Schmertmann,  1979a),  with  0.61  commonly  used  for  relative 
densities  between  40  and  80%  (Seed,  et  al.,  1983).  Therefore  different  stress 
exponents  for  CPT  and  SPT  normalization  have  been  inferred  from  chamber  test. 
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There  appear  to  have  been  no  theories  developed  to  explain  why  the  stress  exponents 
for  Equation  (7.1)  decrease  with  increased  relative  density  or  why  there  is  a 
difference  between  CPT  and  SPT  chamber  derived  stress  exponents. 


SPT  Blow  Count  Normalization 

As  background,  SPT  chamber  data  were  evaluated  in  Chapter  5  and  show  that  the 
SPT  Stress  Focus  location  is  sand  type  dependent.  Also  shown  were  correlations  of 
SPT  chamber  derived  stress  exponent  versus  relative  density.  The  SPT  chamber 
testing  procedures  used  by  Waterways  Experiment  Station  (WES)  were  also 
described  in  that  chapter. 


Stresses  at  the  SPT  Sampler 

Distinguishing  the  stress  states  developed  during  penetration  at  and  nearby  the 
SPT  sampler  is  important  for  the  discussions  to  follow.  These  stresses  are  illustrated 
in  .Figure  7.3.  The  SPT  blow  count  reflects  a  complex  combination  of  static  and 
dynamic  forces  acting  at  the  end  and  along  the  side  of  an  SPT  sampler.  The 
combination  of  all  these  forces  on  the  SPT  sampler  determines  the  SPT  blow  count. 
The  SPT  blow  couni  is  also  dependent  on  confining  stress  as  reflected  in  the  Cn 
normalization  technique.  Therefore,  because  the  SPT  blow  count  is  proportional  to 
the  confining  stress,  the  stresses  influencing  the  SPT  sampler  (crsj )  arc  dependent  on 
the  confining  stress  (asc )  adjacent  to  the  SPT  sampler.  If  the  mud  pressure 
influences  the  confining  stress  nezt  to  the  SPT  sampler  (ctsc  )  then  the  stresses  that 
influence  the  SPT  blow  count  (crsi )  will  be  affected. 
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EFFECTIVE  OVERBURDEN  PRESSURE,  k?ps/sq 


Si- 


Figure  7.1 


Cr  versus  vertical  effective  stress  (Kips/ft2)  rclationsliip  developed  by 
Seed  et  al.  1983  (Note:  2  Kips/ft7  =  1  atm) 
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Figure  7,2  Back-calculated  SPT  CN  pararn  Uts  from  field  and  laboratory  chamber 
test  (Skempton,  19116) 
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Mud  Stress  Influence  at  the  SPT  Sampler 


In  the  field,  the  borehole  mud  pressure  is  approximately  equal  to  the  soil  vertical 
effective  stress.  If  the  borehole  mud  pressure  is  low  (which  is  true  in  chamber 
testing)  in  comparison  to  the  soil  vertical  effective  stress,  then  ihe  confining  stress 
(osc)  next  to  the  SPT  sampler  will  be  reduced.  It  is  outside  the  scope  of  this 
dissertation  to  determine  analytically  the  exact  reduced  stress  level  at  the  sampler  due 
to  reduced  mud  pressure  effects;  however,  an  approximation  will  be  examined.  At 
the  bottom  of  the  chamber  borehole,  the  reduced  mud  pressure  decreases  the  effects 
of  the  in  situ  vertical  effective  stress  around  the  SPT  sampler.  The  typical  5-inch 
diameter  borehole  bottom  can  be  thought  of  as  an  equivalent  "reversed"  circular 
footing  exerting  an  upward  stress  equal  to  the  difference  between  the  in  situ  vertical 
effective  stress  and  the  mud  pressure  as  shown  in  Figure  7.3.  At  a  SPT  sampler 
penetration  of  18  inches  (1 14  feet),  the  ratio  of  equivalent  footing  diameter  to  depth 
is  3 14 j  this  results  in  a  "reversed  footing"  stress  transmission  of  approximately  5% 
based  on  the  simplistic  Boussinesq  stress  distribution  theory  (Larobe  and  Whitman, 
1969).  However,  the  SPT  blow  count  is  determined  by  advancing  the  SPT  sampler 
from  6  to  1 8  inches  below  the  bottom  of  the  borehole.  At  a  depth  of  6  inches,  the 
ratio  of  equivalent  footing  diameter  to  the  depth  is  1.1  which  corresponds  to  a 
"reverse  footmg"  stress  transmission  of  appro ximuieiy  30%.  Therefore,  the  mud 
pressure  effects  on  the  SPT  sampler  end  bearing  forces  during  penetration  starts  at  a 
transmission  of  30%  and  reduces  to  5%  at  the  end  of  penetration. 


The  side  frictional  forces  generated  during  18  inches  of  SPT  sampler  penetration 
arc  also  affected  by  tin;  mud  piessure  effects.  The  SP1  sampler  side  friction 
developed  from  0  to  6  inches  below  the  bottom  of  the  borehole  dominates  the 
sampler  side  frictional  force  dining  sampler  pcncUation  from  6  to  18  inches 
(Schmerimann,  1979a).  Moreover,  the  side  friction  at  18  inches  below  tire  bottom  of 
the  borehole  will  only  influence  the  final  blow  count.  Side  friction  contributions 
from  lire  first  6  inches  arc  sevejal  times  more  influential  than  fiom  the  last  6  inches 
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of  penetration.  Sampler  side  friction  influence  therefore  decreases  with  depth  below 
the  bottom  of  the  borehole.  As  a  result,  the  mud  pressure  transmission  from  the 
equivalent  mud  pressure  "reverse  fooling"  on  the  SPT  side  friction  is  probably 
greatei  than  60%  because  the  first  6  inches  of  side  friction  is  so  dominant. 

The  SPT  sampler  is  resisted  during  penetration  by  a  combination  of  end  bearing 
and  side  frictional  forces  (Schmertmann,  1979a,  and  Douglas,  et.al,  1981).  Reduced 
mud  pressure  influences  will  affect  both  of  these  SPT  sampler  forces.  The  stresses 
surrounding  the  sampler  (c$c)  should  therefore  be  reduced  by  an  amount  equal  to  10 
to  50%  of  the  difference  between  the  in  situ  vertical  effective  stress  and  the  mud 
stress  at  the  bottom  of  the  bore  hole  (i.e.  mud  pressure  "reverse  footing"  effect). 

This  10  to  50%  range  is  based  on  the  60%  value  given  for  side  frictional  force 
in,'  rence  and  the  5  to  30%  range  for  end  bearing  forces.  For  the  immediate 
discussion  to  follow,  an  arbitrary  30%  reduction  will  be  assumed.  However,  other 
mud  pressure  reduction  factors  will  also  be  evaluated  at  the  end  of  the  SPT 
normalization  section. 


Stresses  at  the  SPT  Sampler  for  Field  In  Situ  conditions 

For  a  field  condition,  the  mud  pressure  at.  the  bottom  of  a  bore  bole  shown  as 
line  B  in  Figure  7.4.  The  vertical  effective  stress  is  shown  as  line  V.  The  calculated 
confining  stress  at  the  SPT  sampler  for  the  field  condition,  (as,.)f,  is  shown  as  line  F 
using  the  30%  mud  pressure  influence.  The  resulting  confining  stress  on  the  SPT 
sampler  (as,.),  is  linear  and  very  close  to  the  vertical  effective  stress.  Therefore,  mud 
pressure  does  not  significantly  influence  the  SPT  sampler  for  field  conditions. 
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1 3  S 


Vertical  Effective  Stress  (atm  units) 


Figure  7.4  Drilling  mud  pressure  effects  lor  lire  SPT  sampler  in  Held  boreholes 


Stresses  at  the  SPT  Sampler  for  Chamber  Testing  Conditions 

SPT  laboratory  chamber  tests  have  a  constant  borehole  mud  height  for  all 
chamber  overburden  stresses.  This  mud  pressure  is  always  equal  to  a  4  to  6  feet 
column  of  mud  having  an  approximate  pressure  of  0.16  tsf  as  shown  as  Line  M  in 
Figure  7.5.  The  calculated  confuting  stress  at  the  SPT  sampler  for  the  chamber 
condition,  (orsc)c ,  is  shown  as  line  C  using  the  30%  mud.  pressure  influence 
assumption  previously  discussed.  The  slope  of  Line  C  on  the  log  log  plot  in  terms 
of  vertical  effective  stress  is  0.70.  Line  C  represents  the  relationship  between 
confining  stress  surrounding  the  sampler  (for  chamber  tests)  (  c6C)  c  and  of 

SPT  Chamber  tests  yield  blow  count  versus  ov  relationstiips  that  are  analogous 
to  the  stress  influencing  the  SPT  sampler  (cA  versus  a'..-  From  historic  SPT 

-  u<  *  -  V 

chamber  tests,  this  relationship  between  (csi)c  and  0'v  has  a  SPT  stress  exponent  of 

approximately  0.50,  The  0.50  stress  exponent  is  t fie  average  of  0.45  to  0.55  by  Seed 
et  al.  (1983)  or  is  die  average  of  the  range  at  a  relative  density  of  60%  from 
Chapter  5.  A  stress  exponent  is  simply  the  log  log  slope  as  shown  as  Line  L  in 
Figure  7.5  for  the  (asi)c  to  of 


All  the  stress  relationships  that  influence  tb  :  SPT  sampler  have  been  described 
and  are  in  terms  of  q'v  .  Line  C  (mud  pressure  influence)  describes  the  reduced 

confining  stress  next  to  the  SPT  sampler  (os..)  c,  and  line  L  (from  interpretation  of 
chamber  tests)  describes  the  stresses  on  the  SPT  sampler  that  influence  the  blow 
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Vertical  Effective  Stress  (atm  units) 


N, 


N 


(o;) 


0.43  to  0.53  (using  0.5  lor  illustration) 


Figure  7.5  Stress  effects  on  the  SPT  sampler  for  historical  laboratory  chamber 

tests  due  to  a  constant  drilling  mud  stress  at  differing  confining  stress 
levels 
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count,  (osi)  c.  Both  are  in  terms  of  vertical  effective  stress.  Line  C  can  be 
represented  as; 


=  (0^ 
(o'/70 


(7.2) 


where 

(^sc)c 

~ 

SPT  sampler  confining  stress  level 

for  the  chamber  condition 

c, 

= 

normalized  parameter  (equivalent  value  at  1  Atm) 

Likewise,  Line  L  can  be  represented  as: 

L  =  (7.3) 

(o'/50 

where 

(asi)c  =  £PT  sampler  influencing  stresses 

for  the  chamber  condition 

Lj  =  normalized  parameter  (equivalent  value  at  1  Atm) 
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Equations  (7.2)  and  (7,3)  can  be  combined  to  determine  the  ratio  of  stresses 
influencing  SPT  sampler  (i.e.  blow  count)  to  the  confining  stress  next  to  the  SPT 
sampler: 


^1  = 


(°sX 


1 

0.50 

ll  L>  1 

/ 

with  further  reduction  and  condensing: 


*1  = 


(°Jc 
(  \G-71 

(°Jc 


(7.4) 


(7.5) 


where 

X,  =  combined  normalized  effect 


This  formulation  relates  the  stresses  influencing  the  SPT  blow  count,  iosi)  c,  to 
the  confining  stress  next  the  SPT  sampler  in  the  chamber,  (a8C)  c  with  a  stress 
exponent  of  0.71.  Equation  (7.5)  therefore  relates  the  SPT  blow  count  (as  reflected 
by  (osi  )c)  to  the  actual  confining  stresses  at  the  SPT  sampler  ( (asc)c ).  For  the 
field  SPT  condition,  the  confining  stress  next  to  the  device  is  approximately  equal  to 
Ov-  In  both  cases,  SPT  blow  count  is  related  to  the  confining  stresses  surrounding 
the  SPT  sampler  as  shown  in  Figure  7.6.  For  the  field  SPT  condition,  the  confining 
stress  is  equal  to  the  vertical  effective  stress  (o'v)  However,  for  chamber  SPT  tests, 

the  confining  stress  is  equal  to  the  mud  pressure  reduced  confining  pressure  ( (asc)c). 
The  ultimate  goal  for  field  SPT  normalization  is  a  relationship  of  SPT  blow  count 
(N)  to  cv  Equation  (7.5)  reflect  this  N  to  Ov  relationship  because  the  chamber 

reduced  confining  stress  (osc)c  influences  the  SPT  sampler  just  like  the  field  o  'v 
influen  es  the  S^T  sampler  Therefore,  the  stress  exponent  of  0.71  in  Equation  (7.5) 
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(rather  than  0.5  for  Line  L  in  Figure  7.5)  is  the  correct  stress  exponent  for  field  SPT 
situations.  The  next  few  paragraphs  will  show  that  the  SPT  field  condition  (i.e.  SPT 
stress  exponent)  is  equal  to  the  CPT  chamber  conditions  (i.e.  CPT  stress  exponent). 


A  simple  means  of  determining  the  stress  exponent  for  (osc)  c  based  on  the  mud 
pressure  reduction  value  and  the  chamber  test  stress  exponent  can  be  developed 
based  on  logarithmic  construction  with  the  result  shown  below: 


for 


1‘TbJ 

ri  (°sil 


(7.6) 


(7.7) 


where 

b 

= 

Stress  exponent  for  confining  stress  adjacent  to  SPT  sampler 

n 

Stress  exponent  determined  from  SPT  chamber  tests 

(from  historic  SPT  chamber  test.,) 

K 

Percent  reduction  of  confining  stress  on  the  SPT  sampler  due  to 

the  mud  pressure 

Using  the  initial  example  of  n=0.5  and  1^=30%,  the  result  is  b=0.7 1  using 
Equation  (7.6).  For  a  loose  sand,  the  mud  pressure  reduction  (K^,)  at  the  SPT 
sampler  could  be  20%  (low  end  of  the  range  of  5  to  50%  on  page  137),  the  chamber 
stress  exponent  could  be  n=0.6  (see  Chapter  5),  and  the  result  is  calculated  to  be 
b=0.75.  For  a  dense  sand,  the  reduction  factor  could  be  11^=45%  (high  end  of  tire 
range  of  5  to  50%),  the  chamber  stress  exponent  could  be  n=0.22  (see  Chapter  5), 
and  the  result  is  calculated  to  be  b-0.40.  The  average  stress  exponent  (n)  for  the 


Equvalent  Field  Field 
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column 
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confining  stress  at 
the  SPT  sampler 
due  to  mud 
pressure  effects 


Figure  7.6  Comparison  of  confining  stress  surrounding  the  SPT  sampler  (or  CPT 
probe)  for  chamber  and  field  conditions 
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SPT  chamber  stresses  influencing  the  sampler  ((osi)c)  to  for  these  three 

examples  is  approximately  0.44;  which  is  at  the  lower  range  of  the 

Seed,  et  al.  (1983)  range  of  0.45  to  0.55.  The  stress  exponent  (b)  for  the  SPT  feld 

stresses  influencing  the  sampler  ((osi)  f)  to  o|,  for  these  examples  is  approximately 

0.62;  which  is  within  the  range  of  0.6  to  0.7  from  CPT  cone  resistance  chambei 
tests.  For  these  examples,  the  difference  between  stress  exponents  (n)  and  (b)  are 
specifically  0.21,  0.15,  and  0.18,  with  an  average  of  0.18. 

The  difference  between  the  CPT  chamber  stress  exponents  (c)  and  SPT  chamber 
derived  stress  exponent  (n)  are  shown  in  Figure  7.7  (taken  from  Chapter  6).  This 
stress  exponent  difference  (i.e.  c-n)  in  Figure  7.7  is  about  0.18,  which  is  also  the 
approximate  difference  between  the  SPT  stress  exponent  (n)  and  (b)  from  the 
previous  paragraph.  If  c-n  approximately  equals  b-n  then  c=b.  Therefore,  the  cone 
resistance  stress  exponent  (c)  equals  the  field  based  SPT  stress  exponent  (b).  This  is 
additional  evidence  that  the  CPT  cone  resistance  stress  exponent  (c)  is  the  correct 
SPT  stress  exponent  for  field  SPT  data  normalization. 

Conclusion 

The  Seed  SPT  stress  exponents  of  0.45  to  0.55  are  in  error  because  of  mud 
pressure  influences  for  chamber  tests  at  high  confining  stresses.  Drilling  mud 
pressures  at  the  SPT  sample  in  chamber  tests  are  too  low  compared  to  those  that 
exist  in  the  field.  The  best  stress  exponent  for  SPT  normalization  is  the  CPT  cone 
resistance  stress  exponent.  The  CPT  cone  resistance  stress  exponent  also  falls  within 
the  range  of  back-calculated  CN  in  Figure  7.2.  Chapter  6  introduced  a  technique  and 
chart  for  determining  the  CPT  cone  stress  exponents  using  field  CPT  data  from 
which  the  SPT  normalization  stress  exponents  can  be  taken  as  weii.  To  determine 
the  best  stress  exponer'  for  SPT  stress  normalization  therefore  requires  a  nearby  CPT 
sounding  where  CPT  stress  exponent  can  be  estimated. 
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Normalization  of  the  Shear  Wave  Velocity 


Introduction 

CPT  prediction  of  shear  wave  velocity  requires  a  good  technique  for 
normalization  of  the  shear  wave  velocity — and  the  shear  wave  velocity  normalization 
requires  a  technique  for  estimating  the  shea;-  wave  velocity  stress  exponent. 
Correlations  for  CPT  prediction  of  shear  wave  velocity  will  be  developed  in 
Chapter  8.  A  shear  wave  velocity  normalization  formulation  must  be  simple 
(following  the  stress  normalization  concepts  from  Chapter  3)  and  be  relatable  to 
previously  proposed  formulations.  A  variable  shear  wave  velocity  stress  exponent 
will  be  introduced  and  shown  to  be  dependent  on  soil  type,  and  predictable  using  the 
CPT  cone  resistance  stress  exponent. 

This  shear  wave  velocity  normalization  fonnulation  should  also  be  based  on  the 
maximum  shear  modulus  formulation  because  the  two  arc  theoretic  ally  related.  The 
fust  step  will  be  to  introduce  the  maximum  shear  modulus  formulation  and  then 
show  that  die  historical  formulations  can  be  related  to  it.  The  next  step  is  converting 
the  normalized  maximum  shear  modulus  formulation  to  the  normalized  shear  wave 
velocity  fonnulation.  The  final  step  is  to  correlate  historical  stress  exponents  for  the 
wave  velocity  to  the  stress  exponent  for  the  CPT  cone  resistance. 
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Definition  of  Normalized  Maximum  Shear  Modulus 


Using  the  stress  normalization  concepts  in  Chapter  3,  the  normalized  maximum 
shear  modulus  can  be  defined  as  (Olsen,  1988): 


^max  (^max)j  ( °v)atm 


(7.8) 


where 

Gm(lX  =  shear  modulus  (in  Atm  units) 

(Gmax)i  ~  Normalized  shear  modulus  at  an  equivalent  vertical  effective 

stress  of  3  atmospheric  pressure 

m  -  Shear  modulus  exponent  value 

This  equation  requires  verification  that  it  does  represent  a  generalized  maximum 
shear  modulus  formulation.  Verification  is  shown  in  Table  7.1;  (Gmax),  are  shown 
that  are  derived  from  historic  maximum  shear  modulus  formulations.  The  next  step  is 
converting  (Gmax),  to  a  normalized  shear  wave  velocity  (VsI)  using  theoretical 
relationships. 


Relating  Maximum  Shear  Modulus  to  Shear  Wave  Velocity 


For  linear  elastic  behavior,  the  maximum  shear  modulus  is  theoretically  related  to 
shear  wave  velocity,  Vs,  (Telford,  et.al.,  1976)  as  shown  below: 


G 


mu 


(7.9) 


where 

v. 

= 

Shear  wave  velocity  (in  units  consistent  with  Gmax  and  p) 

P 

= 

Mass  density 
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This  equation  can  be  combined  with  the  normalized  shear  modulus  relationship 
(Equation  (7.8))  to  produce  the  following; 


(Gm")j  Mann 


(7.11) 


where 

(o')atj]  =  Vertical  effective  stress  in  tei ms  of  Atmospheric  pressure 

(mean  stress  (omean)  is  historically  used  in  thir  place  however 
amean  is  simPb’  equal  to  Ftnean  (o'v)  (See  Chapter  2)  ) 
m  =  Stress  exponent  for  Gmax 


With  simple  rearrangement,  shear  wave  velocity,  Vs ,  is  equal  to: 


(7.12) 


This  is  the  basis  for  the  normalized  shear  wave  velocity  formulation.  The 
normalized  shear  wave  velocity,  Vsl,  from  Equation  (7.12)  is  defined  as: 


(7.13) 


with  the  shear  wave  velocity  stress  exponent  (v)  equal  to  m/2.  The  normalized  shear 
wave  velocity  formulation  is  therefore  defined  as  the  following,  based  on  Equation 
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(7.12)  and  using  Equation  (7.13)  with  v=m/2: 


(7.14) 


where; 

v  - 

Shear  wave  velocity  exponent  =  — 

Nonnalized  shear  wave  velocity  in  teims  of  feet/second 

(Vsl  can  be  defined  in  any  velocity  units) 

Relationships  of  Shear  Wave  Velocity  Stress  Exponent  to  Soil  Type 

Lee  and  Campbell  (1985)  presented  generalized  relationships  of  shear  wave 
velocity  versus  depth  based  on  15  years  of  project  work.  This  summary  of  shear 
wave  data  obtained  for  differing  site  conditions  are  in  terms  of  log  she<?.»  wave 
velocity  versus  log  depth.  The  publication  does  not  present  the  actual  data;  however, 
it  presents  the  general  trends  and  range  based  on  the  data.  An  example  for  firm 
natural  soils  is  shown  in  Figure  7.8.  These  shear  wave  velocity  verses  depth  trend;: 
were  also  represented  as  shown  below'; 


V  =  K  (d  c)n  (7.15) 

where 

Vs  =  Siiear  wave  velocity,  in  units  of  ft/second 
d  -  depth,  units  of  feet 

c  =  deptli,  accounts  for  the  non-linear  intersection  of  shear  wave  velocity 

at  the  ground  surface 
n  =  deptli  exponent 

K  ~  constant  in  terms  of  ips/ft" 
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The  depth  exponent  (n)  can  be  related  vo  the  shear  wave  velocity  stress 
exponent  (v),  because  depth  is  used  to  calculate  the  vertical  effective  stress  from  the 
effective  unit  weight  py')  as  shown  below; 

(«'.)-  *  <f  h^y  -  (y y  (n^y  <7-,6> 

The  vertical  effective  stress  and  depth  parameters  in  Equation  (7.16)  both  have  the 
same  exponent.  Therefore  the  depth  exponent  (n)  in  Equation  (7.15)  is 
approximately  equal  to  the  shear  wave  velocity  stress  exponent  (v),  at  least  for 
normalization  purposes.  Lee  and  Campbell's  (1985)  data  were  summarized  for  soft 
natural  soils,  intermediate  firm  natural  soils,  and  firm  soils  as  shown  in  Table  7.2. 
Tne  depth  exponent  (i.e.  stress  exponent)  from  this  cable  will  be  related  to  soil  type 
and  the  CPT  cone  resistance  stress  exponent  at  the  end  of  this  discussion. 


Table  7.2  Shear'  wave  velocity  parameters  determiner:  from  daur  presented  by 
Lee  and  Campbell  (1985)  for  differing  soil  conditions 


Soil 

V  (ft/sec)  at  a  dcptlt  of 

Depth  exp 

classification 

20  ft 

for  shear 

Soft  natural  soils 

780 

0.46 

Intermediate  frrn 

950 

0.43 

natural  soils  ‘ 

Firm  soils  " 
s  D.-rrdjrr:, 

1300 

0.  1 

**  Probably  dense  sands 
*tCr  Approximately  Vs)  at  o’,  -  1  AUu 
**+*  Approximately  equal  to  the  shear  wave  velocity  slices  exponent  (v) 


Shear  Wave  Velocity  Stress  Exponent  trend  from  a  Sand  Site 

liuldi.  Brtrezi,  Supcrbo,  Baftaglio  and  Jamio.'i'owski  (1988)  stud’s-'1  tlie  measured 
shear  wave  velocities  fcr  Po  Iliver  sand: .  1  here  tando  have  a  fines  content  (in  terms 


DEPTH  !N  FEET 


1,000 
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10 


Figure  7.8 
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- 

FIRM  NATURAL  SOILS  : 


100  1,000  10,000 
Shear  wave  velocity  in  ft/sec 


Sheai  wave  velocity  versus  depth  for  intermediate  and  saturated  firm 
natural  so.  s  (l  ee  and  Canipbell,  1983) 
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% 

of  the  percent  passing  the  #200  sieve)  ranging  from  3  to  9%.  They  related 
overburden  stress  and  maximum  shear  modulus  in  a  form  indexable  to  Equation  (7.8) 
but  also  related  fines  content  to  the  stress  exponent  (m)  for  the  maximum  shear 
modulus  as  shown  below: 


m  -  0.43  +  0.39 


viooj 


(7.17) 


where 

m  =  Stress  exponent  for  the  maximum  shear  modulus 

P200  =  percent  passing  the  #200  sieve  as  a  percent 


This  equation  can  be  expressed  for  the  shear  wave  velocity  stress  exponent  (v)  as 
shown  below,  because  (v)  is  theoretically  equal  to  half  of  (m). 


v  -  0.215  +  0.195 


\ 


liooj 


(7.18) 


For  a  clean  sand  (i.e.  P2oo=0%)  the  shear  wave  velocity  stress  exponent  (v)  is  0.21 
and  for  a  sand  with  10%  fmes  content  (P200=10%)  the  stress  exponent  (v)  is  0.24. 
Therefore  for  tire  transition  from  clean  sand  to  dirty  sand,  the  shear  wave  velocity 
stress  exponent  changes  from  0.21  to  0.24. 


Relating  Stress  Exponents  for  Shear  Wave  Velocity  &  Cone  Resistance 

Table  7.3  presents  additional  data  relating  the  soil  type  to  the  shear  wave  velocity 
stress  exponent  from  published  sources  and  this  research.  The  table  indicates  that  the 
shear  wave  velocity  stress  exponent  for  clay  ranges  between  0.42  to  0.50  and  for 
sands  ranges  between  0.18  and  0.28.  These  ranges  from  Table  7.3  together  with  the 
Baldi  et  al.(1988)  and  Lee  and  Campbell  (1985)  data  arc  summarized  in  Figure  7.9 
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in  terms  of  soil  type.  Also  shown  in  this  figure  are  typical  cone  resistance  stress 
exponents  from  Chapter  6  based  on  soil  type.  The  trends  from  Figure  7.9  indicate 
that  the  shear  wave  velocity  stress  exponent  (v)  is  approximately  45%  of  the 
estimated  cone  resistance  stress  exponent  (c).  While  this  stress  exponent  prediction 
technique  can  only  be  considered  an  estimate,  it  represent  the  best  method  when 
there  are  no  uniform  layers  to  support  establishment  of  a  stress  trend. 


Table  7.3  Published  stress  exponent  for  shear  wave  velocity 


Basic 
soil  type 

Specific 
Soil  type 

Vsl 

(ft/sec) 

Stress 

exponent, 

V 

Reference 

Sand 

1 

Review  of 
published 
dry  sands 

■ 

m 

Ottawa  sand 

ermalinn 

--j  - 

Hardin  Richart  fl963i 

—  *  "  J  '  ' ' ' 

equation 

aiii 

m 

u  m 

coarse 

quartz 

041 

Hamilton  (1971) 

Clays 

(Normally 

Consolidated) 

0.5 

SMI 

mm 

Clays 

(Over 

Consolidated) 

over 

consolidated 
Alaskan  clay 

HKI 

Harding  Lawson  Assoc 
(1978) 

i 
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Stress  exponent 


Normally  consolidated  clay 


>!^e 


Medium 

dense 

sand 


The  shear  wave  velocity  stress  exponent  (v) 
is  approximately  equal  to  45%  of  the 
cone  resistance  stress  exponent  (c). 


Clay  Soil  Mixtures  Sand 

Soil  Classification 


Legend 


Cane  resistance  stress  exponent  trends  from  Chapter  5 
J  Stokoc,  Lee,  and  Knox,  1985,  and  Table  IT 


Lee  and  Campbell  (1985) 
Baldi,  ehal.,  (1988) 


Figure  7.9  Comparing  the  stress  exponent  ranges  for  shear  wave  velocity  and 
cone  resistance 
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Conclusions 


The  shear  wave  velocity  stress  exponent  for  normalization  is  related  to  soil  type 
and  can  be  estimated  based  on  the  CPT  cone  resistance  stress  exponent.  The  shear 
wave  velocity  (Vs)  stress  exponent  (v)  is  approximately  45%  of  estimated  cone 
resistance  stress  exponent  (c)  (from  Chapter  6).  This  procedure  represents  a  new 
approach  for  r^nnalizing  shear  wave  velocity  and  seems  to  give  good  results.  To 
successfully  normalize  the  shear  wave  velocity,  therefore,  requires  a  nearby  CPT 
sounding  where  the  CPT  stress  exponent  can  be  estimated  and  then  converted  to  the 
shear  wave  velocity  stress  exponent  (v). 
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Chapter  8 


Development  of  Correlations  for 
CPT  Prediction  of  Geotechnical  Properties 


Introduction 

One  of  the  ultimate  goals  of  using  the  CPT  as  a  geotechnical  investigation  tool  is 
direct  prediction  of  geotechnical  properties.  Prediction  of  geotechnical  properties 
using  CPT  data  is  illustrated  in  Figure  8.1;  CPT  data  is  normalized,  the  normalized 
geotechnical  property  is  predicted,  and  finally  the  geotechnical  property  is  calculated 
for  the  in  situ  stress  condition. 

Historically,  research  has  concentrated  on  developing  relationships  between  cone 
resistance  and  various  geotechnical  properties.  In  this  chapter,  CPT  correlations 
using  CPT  cone  resistance  and  sleeve  friction  resistance  are  developed  for  prediction 
of  SPT  blow  count,  clay  undrained  strength,  and  shea  wave  velocity.  Several 
important  ingredients  were  required  for  developing  these  correlations,  namely  1 ) 
stress  normalization  techniques  for  the  CPT  (and  geotechnical  properties),  2)  a  large 
database  of  CPT  and  soil  sample  data,  and  finally,  3)  accounting  for  bias  data  caused 
by  soil  profile  differences  between  CPT  soundings  and  boreholes. 
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CPT 


Boring 


Figure  8.1  Using  CPT  measuremcuis  10  predict  geotechnical  properties 
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The  CPT  and  soil  sample  database 


A  large  database  of  CPT  and  soil  sample  information  that  represents  soils  from 
around  the  world  is  described  in  the  Appendix.  This  database  contains  data  from 
approximately  90  projects,  670  CPT  soundings,  580  borings,  and  approximately  8100 
laboratory  and  field  test  values.  The  soils  include  weak  clays  in  San  Francisco, 
Sweden,  Hong  Kong,  (and  sensitive  clays  in  Norway),  liquefiable  soil  mixtures  from 
China,  stiff  clays  from  South  Carolina,  and  sands  from  Po  River  and  desen  alluvial 
sands  from  Nevada  (this  is  only  a  partial  list  of  the  total  database). 


Using  the  CPT  Soil  Characterization  Chart 
to  Predict  Properties 

Two  measurements  are  always  better  than  one  when  establishing  a  correlation. 
Moreover,  the  CPT  provides  two  measurements  which  are  unique,  repeatable,  and 
accurate.  The  CPT  soil  characterization  chart  shown  in  Figure  8.2  (Olsen,  1988)  is 
based  on  normalized  CPT  measurements.  This  chart  provides  the  means  of 
characterizing  soil  behavior  in  terms  of  soil  type  and  relative  strength — the  relative 
strength  is  in  terms  of  increasing  overconsolidation  at  one  extreme  and  increasing 
clay  sensitivity  or  a  metastable  condition  in  sand  at  the  other  extreme.  While  the 
underlying  normalization  technique  for  this  chart  was  substantially  improved  as  a 
result  of  the  current  research  effort,  the  fundamental  aspects  of  this  chart  arc  still 
valid.  This  chart  therefore  provides  the  two  soil  behavior  indexes  based  on  the 
normalized  CPT  measurements,  namely  soil  type  and  relative  strength  consistency. 
The  cone  resistance  stress  exponent  (c)  was  shown  in  Chapters  5  and  6  to  be 
inversely  related  to  the  relative  sUength  (c.g.  relative  density  for  sands).  In 
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Chapter  6,  the  stress  exponent  (c)  contours  on  the  CPT  soil  characterization  chart 
was  demonstrated  to  be  perpendicular  to  soil  classification  contours.  Soil  type  and 
relative  strength  are  independent  of  each  other  and  this  perpendicularity  of  soil  type 
to  stress  exponent  signifies  independence.  Consequently,  any  soil  type  and  relative 
strength  level  combination  will  correspond  to  a  unique  point  on  this  chart. 

Therefore,  contours  of  geotechnical  property  levels  can  be  established  on  the  CPT 
soil  characterization  chart  because  the  chart  characterizes  soil  type  and  relative 
strength. 


Statistics  and  CPT  Correlations 

Statistical  errors  can  be  divided  into  three  categories;  systematic  (bias)  errors, 
random  (variance)  emors,  and  mistakes  (Taylor,  1990).  Random  errors  are 
fluctuation  errors  about  the  mean  and  are  statistically  characterized  by  the  standard 
deviation.  Bias  errors  are  offset  errors  from  the  mean.  For  CPT  data  evaluation, 
bias  errors  are  caused  by  stratigraphic  soil  type  differences  between  CPT  soundings 
and  nearby  boreholes.  This  type  of  bias  error  will  be  shown  to  account  for  most  of 
the  statistical  error  during  CPT  correlations.  However,  there  is  little  written  about 
bias  errors  (Taylor,  1990,  Mosteller  &  Tukey,  1 977,  Huaslin,  Mosteller,  and  Tukey, 
1983).  'Ihercfore  a  subjective  quality  index  and  evaluation  technique  were  developed 
to  account  for  data  having  bias  error.  Reducing  bias  error  effects  represented  the 
primary  statistical  evaluation  effort  used  in  establishing  CPT  correlations  in  this 
chapter. 
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Corrected  Cone  Resistance  q  ,  =  — — r  in  terms  of  tsf 
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Figure  8.2  CPT  Soil  Characterization  chart  in  terms  of  the  normalized  cone 
resistance  and  normalized  friction  ratio  (Olsen  1988  ' 
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Bias  Error  Results  from  Stratigraphy 


There  is  always  some  geological  difference  between  field  CPT  soundings  and 
nearby  boreholes  at  the  same  elevation.  At  a  given  elevation  and  for  short  lateral 
distances  between  a  CPT  sounding  and  a  borehole  (such  as  10  to  25  feet)  the 
following  can  be  generally  assumed;  1)  the  soil  types  may  be  different,  2)  the 
formational  environment  is  generally  the  same,  3)  the  vertical  effective  stress  is 
constant,  and  4)  the  vertical  stress  history'  (i.e.  overconsolidation  level)  is  generally 
the  same.  Soil  type  changes  between  CPT  soundings  and  boreholes  is  therefore 
more  likely  the  cause  of  bias  error  than  differences  in  overburden  stress, 
overconsolidation,  or  lateral  stress.  The  potential  for  soil  type  change  for  a  given 
elevation  interval  is  dependent  on  the  depositional  environment;  e.g.  a  lake  or  tailings 
deposit  generally  has  uniform  layers  and  lenses  while  a  rapidly  flowing,  meandering 
river  produces  the  most  non-uniform  deposit. 

Stratigraphic  soil  type  change  (i.e.  geologic  change)  is  a  bias  error  because  it 
shifts  the  average  and  increases  the  variance — other  types  of  random  statistical  errors 
only  increase  the  statistical  variance  level  without  affecting  the  average.  For 
example,  clay  lenses  within  a  sand  deposit  will  move  the  average  downward  to  the 
clay  trend  and  away  from  a  sand  trend  as  shown  in  Figure  8.3.  Likewise,  silt  or 
sand  seams  within  a  clay  deposit  will  create  a  higher  average  as  also  shown  in 
Figure  8.3.  Bias  error  direction  is  therefore  different  for  clay  and  sand  correlations 
(or  dense/stiff  and  loose/soft  correlations). 


The  Academic  Quality  Index  (AQI) 

The  Academic  Quality  Index  (AQI)  was  developed  during  this  research  program 
as  tool  for  accounting  for  bias  error.  It  also  provides  a  basis  for  weighting  data  from 


164 


Depth 


SANDy 


Clay  Lenses 
influence 


165 


the  different  sites  in  the  development  of  CPT  to  property  correlations.  This  ranking 
scale  is  almost  a  universal  scale  because  it's  based  on  the  student  academic  grading 
scale  understood  by  all  ages  and  professions.  For  example,  data  with  average  quality 
has  an  AQI  of  76%,  good  data  has  an  AQI  of  80  to  85%  and  excellent  data  has  a 
AQI  of  90  to  95%.  This  index,  while  subjective,  does  allow  excellent  quality  data  to 
be  isolated  from  good  or  poor  quality  data. 

An  Academic  Quality  Index  (AQI)  can  be  individually  assigned  for  stratigraphic 
change  potential,  CPT  work  quality,  and  boring/laboratory  quality.  The  overall  CPT 
AQI  is  equal  to  the  stratigraphic  change  potential  together  with  any  decreasing 
influence  of  measurement  quality.  The  overall  CPT  AQI  therefore  reflects  the 
potential  for  lateral  matching  of  the  soil  type  (and  relative  strength)  between  a 
excellent  quality  CPT  sounding  ana  adjacent  quality  borehole.  Some  of  the 
commercial  CPT  data  in  the  database  has  poor  quality,  for  example,  a  high  capacity 
cone  used  in  weak/soft  soil  where  a  low  capacity,  high  accuracy  cone  would  have 
been  more  appropriate.  If  the  CPT  or  laboratory  testing  data  quality  is  not  excellent, 
the  overall  AQI  is  less  than  the  stratigraphic  AQI.  For  example,  if  the  stratigraphic 
continuity  is  excellent  (i.e.  AQI  of  95%)  but  the  CPT  measurements  are  near  the 
equipment  noise  level,  the  overall  AQi  might  only  be  80%.  If  there  is  little 
information  about  site  geology  or  testing  quality,  then  an  overall  AQI  of  76%  is 
assigned.  While  these  are  obviously  arbitrary  distinctions  they  are  easily 
remembered  and  useable. 

The  stratigraphic  AQI  is  based  on  a  simple  arbitrary  system  equal  to  the 
estimated  percent  of  soil  depth  which  are  continuous  between  a  CPT  sounding  and 
nearby  borehole.  For  example,  if  76%  of  the  soil  layers  are  continuous  between  a 
C  vr  sounding  and  borehole  then  the  stratigraphic  AQI  is  76%  as  shown  in 
Figure  8.4.  This  sti atigraphic  AQI  must  also  account  for  the  CPT-to-borehole 
distance  because  AQI  will  increase  us  the  lateral  distance  decreases.  For  example,  if 
the  AQI  is  equal  to  76%  for  a  lateral  distance  of  20  feet,  then  the  AQI  might 
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Definition  of  a  stratigraphic  AQI  of  76%  (i.e.  76%  of  the  soil  layer 
are  continuous  between  the  boring  and  CPT  sounding) 


increase  to  85%  for  a  distance  of  10  feet.  Research  sites  with  thick  uniform 
horizontal  deposits  can  have  a  stratigraphic  AQ1  of  95%  for  CPT-to-boreholc 
spacings  of  20  feet. 

Overall  AQI  values  were  assigned  for  each  geotechnical  project  in  database  based 
on  the  concepts  just  described.  The  overall  AQI  values  are  shown  in  the  fourth 
column  of  table  A-l  in  the  Appendix. 


Developing  CPT  Predictive  Contours 


Consider  the  ncn-lineai  response  surface  based  on  two  dependent  parameters  as 
shown  in  Figure  8.5.  The  two  dependent  parameters  are  the  X  and  Y  axis  with  the 
contoured  (predictive)  response  surface  (or  blanket  or  contours)  in  the  Z  direction. 
This  response  surface  is  always  non-linear  in  geotechnical  engineering  and  can  be 
determined  by  fitting  the  data  points  in  the  3  dimensional  space  using  gridding 
techniques  (Box  and  Draper,  1987)  if  the  data  noise  (i.e.  error  distribution)  is 
random.  Contouring  techniques  do  not  account  for  bias  error.  The  next  section  will 
describe  a  teclmique  using  the  AQI  quality  index  for  establishing  the  best  fit 
correlation  for  each  contour  by  excluding  biased  data.  The  response  surface  (i.e.  set 
of  contours)  is  established  by  separately  generating  each  individual  contour  and  then 
combining  all  the  contours.  Each  contour  was  developed  using  2D  data  scatter 
plotting  software  program  while  also  accounting  for  the  shapes  of  the  other  contour 
levels. 
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Predictive 


Figure  8.5  Description  of  a  response  surface 
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Process  for  developing  CPT  Predictive  Contours 


Developing  contours  of  geotechnical  properties  on  the  CPT  soil  characterization 
diagram  required  numerous  steps  that  are  described  in  this  section.  All  data  for  a 
given  geotechnical  property  are  extracted  from  the  database  and  placed  into  a  single 
sequential  computer  file  list  (composed  of  normalized  CPT  measurements  and 
corresponding  normalized  geotechnical  properties).  Also  included  in  each  list  are  the 
AQIs  for  each  CPT  to  geotechnical  property  comparison.  The  list  is  then  imported 
into  a  spreadsheet  for  the  goal  of  creating  scientific  graphics.  This  spreadsheet  data 
list  is  initially  sorted  by  the  data  column  having  the  geotechnical  property.  The 
spreadsheet  data  list  is  then  divided  into  6  to  7  geotechnical  property  level  groups 
(i.e.  100  to  200,  200  to  400,  400  to  800,  etc.,)  for  the  purpose  of  establishing  CPT 
based  correlations  based  on  data  groups.  Each  of  these  geotechnical  property  level 
groups  is  then  individually  sorted  based  on  the  AQ!  column.  This  final  sort 
differentiates  potentially  biased  data  from  good  data  for  each  geotechnical  property 
level  group.  To  establish  a  trend  for  a  given  geotechnical  property  level  data,  all  or 
part  of  the  data  within  each  data  group  can  be  plotted.  The  best  correlations  (i.e. 
contours)  are  established  by  using  the  highest  quality  data  which  means  using  the 
highest  AQI  level  for  each  data  group. 


Use  of  the  AQI  for  Developing  Predictive  Contours 

A  low  stratigraphic  AQI  implies  a  higher  likelihood  for  soil  type  change  between 
a  CPT  sounding  and  borehole.  When  a  group  of  data  having  a  minimum  AQI  of,  for 
example,  80%  is  used  for  a  correlation  then  this  is  defined  as  the  "minimum 
inclusionary  AQI"  of  80%.  If  all  data  is  included  for  a  correlation  then  the  minimum 
inclusionary  AQI  is  zero.  By  increasing  this  minimum  inclusionary  AQI,  low  quality 
biased  data  is  excluded  and  the  predictive  correlation  will  shift  away  from  the  biased 


170 


data.  If  the  minimum  AQI  is  raised  too  high,  there  will  be  insufficient  quantity  of 
data  to  develop  a  trend.  For  a  particular  data  set,  the  optimum  AQI  might  be  85% 
but  if  this  minimum  AQI  is  increased  to  95%,  there  might  only  be  2  points  which 
may  not  be  enough  data  to  establish  a  trend.  For  a  given  predictive  property  the 
minimum  inclusionary  AQI  must  therefore  be  carefully  increased  until  the  optimum 
trend  is  esm'  fished. 


An  illustrated  example  will  show  how  the  inclusionary  AQI  concept  is  used  to 
determine  a  best  fit  contour.  Consider  a  hypothetical  data  range  for  a  normalized 
shear  wave  velocity  (Vs])  of  300  to  500  feet/second  (fps)  (with  an  average  of  400 
fps)  shown  in  Figure  8.6-a.  The  displayed  data  points  represent  data  having  a  Vs] 
between  300  and  500  fps  and  are  plotted  in  terms  of  qcle  versus  FRle.  A  Vsl  of  400 
fps  can  be  considered  at  the  low  range  typically  considered  in  geotechnical 
engineering  (i.e,  300  to  2000  fps).  Therefore,  biased  data  for  400  fps  data  will  tend 
to  skew  the  average  to  higher  qc!c  and  f£,e  levels.  The  total  data  in  Figure  8.6-a  can 
be  divided  into  two  inclusionary  AQI  ranges  (as  shown  in  Figure  8.6-b)  in  terms  cf 
the  data  range  and  average  for  each  inclusionary  AQI  range.  With  increasing 
inclusionary  AQI  ranges,  the  lower  quality  data  is  excluded  ai.J  the  average  moves 
away  from  the  skewed  group  average.  Figure  8.6-c  is  probability  distribution  along  a 
cut  through  the  data  set  in  Figure  8.6-b  to  show  how  increasing  the  inclusionary 
AQI  moves  the  data  average  by  excluding  biased  data. 

The  next  example  is  taken  from  the  next  section  on  CPT  prediction  of  the 
normalized  SPT  N,.  All  data  from  the  database  for  SPT  N,  ranging  from  1  to  3 
(average  of  2)  are  shown  in  Figure  8.7.  For  sands,  a  SPT  N,  of  2  is  considered  a 
loose  sand  and  therefore  can  be  biased  by  dense  and  medium  dense  sand  layers. 

When  the  minimum  inclusionary  AQI  is  raised  to  89%  (shown  with  solid  circles  in 
Figure  8.7),  only  the  best  data  are  included  and  clearly  show  a  correlation  at  the 
lower  boundary  of  the  total  data  scatter. 
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Normalized  Cone  Resistance 


(a) 


(b) 


Limits  of  the  data 


(c) 

Figure  8.6  An  ihustration  of  the  iterative  technique  for  removing  bias  data  using 
the  AQI  index 
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Figure  8.7  Example  using  the  Academic  Quality  Index  (AQI)  toward  establishing 
the  best  correlation  line  by  excluding  biased  data 
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CPT  Prediction  of  the  SPT  blow  count 


Introduction 

Prediction  of  the  Standard  Penetration  Test  (SPT)  blow  count  using  CPT  data  is  a 
complex  task  because  the  SPT  sampler  resistance  during  penetration  is  dependent  on 
the  combination  of  end  bearing  and  side  friction  forces  which  are  dependent  on  soil 
type  and  soil  strength  (Douglas,  et.al.,  1981).  The  SPT  sampler  is  resisted  by  the 
same  types  of  forces  that  are  measured  by  the  CPT,  therefore  using  both  CPT 
measurements  to  predict  the  SPT  blow  count  (Schmertmann,  1979a)  is  a  better 
approach  than  using  only  the  ratio  of  cone  resistance  to  SPT  blow  count  (e.g.  qc/N). 
The  forces  acting  on  the  SPT  sampler  are  reflected  in  the  final  CPT-to-SPT 
correlation  if  both  CPT  measurements  are  used  for  the  evaluation. 

It  was  shown  in  Chapter  7  that  the  historic  SPT  stress  exponents  (i.e.  0.45  to 
0.55)  for  normalization  are  incorrect  and  the  actual  values  are  equal  to  the  CPT  cone 
resistance  stress  exponents.  This  3lso  implies  that  the  SPT  normalization  is 
dependent  on  soil  type  and  relative  strength  consistency  as  is  the  CPT  normalization 
(Chapter  6). 

Historically,  there  have  been  two  general  techniques  for  CPT  prediction  of  the 
SPT  blow  count;  1)  the  q^N  ratio  to  predict  the  SPT  N  (Rodin,  et.al.,  1974, 
Robertson,  et.al.,  1983,  Seed  and  De  Alba,  1986,  Kuihawy  and  Mayne,  1990),  and 
2)  using  both  CPT  measurements  to  predict  the  normalized  SPT  N,  (Olsen,  1984, 
1986,  1988).  What  direrentiates  these  two  techniques  that  the  first  technique  is 
based  on  an  empirical  correlation  and  the  second  technique  indirectly  accounts  foi 
the  end  bearing  and  side  friction  forces  that  act  on  the  SPT  sampler. 
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SPT  Prediction  using  Both  CPT  Measurements 


Research  indicates  that  the  resisting  stresses  on  the  SPT  sampler  are  analogous  to 
the  CPT  stresses  as  shown  in  Figure  8.8  (Schmertmann,  1979a).  McLean,  Franklin, 
and  Dahlstrand  (1975)  developed  a  computer  model  of  the  dynamic  forces 
influencing  SPT  sampler  based  on  work  by  Schmertmann  (1971).  Douglas,  Olsen 
and  Martin  (1981)  also  developed  a  computer  model  to  predict  SPT  using  the  CPI 
measurements  based  on  work  by  Schmertmann  (1979a,  1979b)  in  terms  of  the  static 
and  dynamic  forces  on  the  SPT  sampler.  This  work  evolved  to  a  technique  of 
predicting  the  SPT  blow  count  which  uses  both  CPT  measurements  (Olsen,  1984) 
with  the  most  recent  published  version  shown  in  Figure  8.9  (Olsen,  1988).  However, 
in  1 987,  there  was  no  realistic  means  of  showing  the  data  that  was  used  to  develop 
the  contours  of  SPT  N,  on  the  CPT  soil  characterization  chart.  Also,  the  database 
for  this  research  program  is  approximately  7  times  larger  than  existed  in  1987  for 
developing  SPT  N,  contours. 
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Figure  8.9  Prediction  ol  the  SPT  blow  count,  N,  using  both  CP  I"  measurements 
(Olsen,  1988) 


CPT  Prediction  of  N}  using  the  CPT  Soil  Characterization  Chart 

Establishing  correlations  of  normalized  CPT  parameters  to  SPT  N,  was  found  to 
be  very  dependent  on  the  AQI  statistical  technique.  AQI  was  used  to  remove  biased 
data  effects  thereby  allowing  the  best  correlation  of  CPT  to  SPT.  All  SPT  data 
from  the  database,  for  all  projects,  were  divided  into  seven  N,  ranges  for  the  purpose 
of  establishing  individual  correlation  contours  on  the  CP  I  soil  characterization  chart. 
The  results  of  the  SPT  N]  correlations  are  shown  in  Figure  8.10  to  Figure  8.16  for 
average  Nj  levels  of  2,  4,  7,  15,  25,  35,  and  50.  The  best  fit  correlation  contour  in 
each  figure  was  based  primarily  on  highest  inclusrnary  AQI  but  also  accounted  for 
the  lower  quality  data  together  with  the  shape  of  the  other  SPT  N,  contour  groups. 
Figure  8.17  shows  all  the  SPT  N,  predictive  contours  together  with  soil  classification 

i : - £t _ ni  _ _ nnoo\ 

lines  Hum  wisen  ^1700). 


Discussion  of  CPT  Prediction  of  the  SPT  Blow  Count 

The  new  contours  for  predicting  normalized  SPT  blow  count  (Figure  8.17)  in 
general  have  more  curvature  compared  to  the  1988  version  (Figuie  8.9);  however,  for 
faction  ratios  less  than  0.5°/e,  the  contours  arc  parallel  to  cone  resistance,  which  was 
not  expected.  The  contour  shapes  and  intervals  for  the  1988  version  are  symmetric 
due  to  the  lack  of  data;  newer  version  contour  shapes  and  intervals  are  not  in  general 
symmetric.  The  shape  of  the  contours  for  the  newer  version  reflect  soil  type  and 
relative  strength  contribution.  Within  the  sands  area  of  the  chart  for  friction  ratios 
less  than  0.5%,  the  contour  shapes  and  intervals  are  symmetric.  Within  the  clay  area 
of  the  chart  (from  normally  consolidated  to  overconsolidated),  the  contour  shapes  and 
intervals  are  also  symmetric.  However,  the  contour  shapes  are  changing  in  the 
middle  of  the  chart  within  the  soil  mixture  area  (transition  between  clays  and  sands). 
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Figure  8.10  Prediction  of  the  N,=2  contour  on  the  CPT  Soil  Characterization  Chart 
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Normalized  Cone  Resistance  (atm) 
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Figure  8.12  Prediction  of  the  N,=7  contour  on  the  CPI  Soil  Characterization  Chart 


Normal 
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Figure  8.15  Prediction  of  the  Nt=35  contour  on  the  CPT  Soil  Characterization 


When  the  CPT  friction  ratio  is  low,  such  as  0.6%,  the  contours  of  predicted  SPT 
Nj  are  almost  parallel  to  the  normalized  cone  resistance  (qcle)  which  indicates  that 
the  SPT  sampler  end  bearing  force  dominates  the  SPT  resistance.  Therefore,  for 
loose  to  medium  dense  sands  (which  produce  low  friction  ratios,  less  than  0.6%) 
(Douglas,  1982)  the  SPT  sampler  is  resisted  primarily  by  end  bearing  forces. 

For  CPT  friction  ratios  greater  than  3%  in  Figure  8.17,  the  predicted  normalized 
SPT  N,  contours  are  approximately  parallel  to  the  normalized  CPT  sleeve  friction 
resistance,  fslc  (contours  of  constant  normalized  sleeve  friction  resistance  are  shown 
in  Figure  8.18).  This  indicates  that  the  SPT  sampler  resistance  is  primarily 
dependent  on  the  side  friction  resistance  for  dense  and  overconsolidated  sands 
together  with  all  clays.  Dense  sands  and  overconsolidated  soils  also  have  high  lateral 
stress  ratios  (KJ.  For  sands,  a  high  lateral  stress  produces  a  higher  CPT  sleeve 
resistance  (Masood,  l9yu)  and  therefore  wouid  also  produce  a  higher  SPT  sampiei 
resistance.  The  relative  contribution  of  SPT  side  friction  forces  to  SPT  end  bearing 
forces  can  be  calculated  from  the  product  of  factors:  the  area  ratio  (i.c.  SPT  sampler 
side  area  to  SPT  end  bearing  area)  and  stress  ratio  (i.e.  the  CPT  friction  ratio).  'In; 
area  ratio  is  the  SPT  side  friction  area  to  the  SPT  end  bearing  area  which  is  equal  to 
28  at  1  Vi  foot  penetration.  For  clays  and  overconsolidated  sand,  the  CPT  friction 
ratio  is  approximately  5%  (i.e.  stress  ratio  of  0.05),  therefore  the  SPT  side  force  to 
end  bearing  ratic  is  calculated  to  be  1.4  (i.e.  28  times  0.05).  The  SPT  side  friction 
force  is  1 .4  times  higher  than  the  end  bearing  force.  However,  if  side  friction  forces 
inside  the  sampler  are  included,  then  the  area  ratio  increases  to  44  and  the  SPT  side 
force  to  end  bearing  ratio  increases  to  2.2. 

An  example  of  CPT  predicted  (using  Figure  8.17)  versus  field  measured  SPT 
blow  counts  is  shown  in  Figure  8.19  together  with  annotations.  Numerous  other 
examples  of  CPT  predicted  versus  field  measured  SPT  blow  counts  are  shown  in 
Figure  8.20  to  Figure  8.27.  Prediction  of  the  SPT  blow  count,  is  difficult  as 
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Figure  8.18  Contours  of  constant  normalized  sleeve  friction  resistance  on  the  CP  I 
soil  characterization  chart 
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illustrated  by  large  variance  for  each  predictive  SPT  blow  count  contour  (for  example 
Figure  8.14).  However,  as  seen  from  these  predicted  versus  measured  examples,  the 
match  of  predicted  versus  measure  does  support  this  teclmique  of  using  both  CPT 
measurements  to  predict  the  SPT  blow  count.  These  figures  also  show  how  much 
detail  is  missed  by  having  only  discontinuous  SPT  measurement  records. 
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Figure  8.19  Example  of  CPT  Prediction  of  SPT  blow  count  with  comparison  data 
(Barkley  dam  seismic  evaluation  study,  data  from  Olsen,  et  al,  1989) 
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Figure  8.23  CPI  predicted  versus  measured  SPT  N,  for  a  stiff/dense  dirty  sand  site 
(project  code  DE-S-R) 
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Figure  8.24  CPT  predicted  versus  measured  SPT  N,  for  a  site  composed  of  soil 
mixtures  in  Tracy  California  (project  code  WC-TSS) 
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Figure  8.25  CPT  predicted  versus  measured  SPT  N,  for  McDonaids  Farm,  Sea 
Island,  Richmond,  British  Columbia,  Canada  (University  of  British 
Columbia  research  site)(data  from  Campanella,  et  al,  1982) 
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Figure  8.27  CPT  predicted  versus  measured  SPT  Nj  at  the  Coyote  North  site  (near 
San  Jose)  California  (data  from  Martin  &  Douglas,  1980) 
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CPT  Prediction  of  Clay  Undrained  Strength,  Su 


Introduction 

The  unchained  strength,  Su,  of  clay  is  one  of  the  classic  geotechnical  properties 
and  is  critical  for  many  design  applications.  In  die  U.S.,  the  standard  practice  for 
obtaining  clay  undrained  strength  requires  retrieval  of  relatively  undisturbed  soil 
samples,  typically  with  a  3  inch  shelby  tube,  for  unconsolidated  undrained  (UU) 
triaxial  testing.  However,  only  a  few  laboratory  strength  tests  are  typically 
performed  for  most  geotechnical  projects  and  it  is  assumed  that  they  represent  the 
character  of  the  total  site. 

In  Europe,  the  standard  practice  for  obtaining  tire  undrained  strength  of  soft  clays 
is  the  field  vane  shear  test.  However,  the  vane  shear  strength  is  always  too  high 
because  of  s'b  v  -utent,  overconsolidation  effects,  tire  high  vane  rotation  rate,  and/or 
vane  geometric,  effects.  A  correction  factor  (p)  is  required  to  reduce  the  field  vane 
shear  strength  to  the  field  strength  level  (typically  the  equivalent  unconsolidated 
undrained  triaxial  test  level).  This  correction  factor  (p)  is  dependent  on  many 
faetors,  but  generally  established  based  on  site  specific  correlations. 


CPT  Prediction  of  Su 

Clay  undrained  strength  was  the  first  geotechnical  property  predic  ted  using  CPT 
data  in  the  1 930's  (reported  by  Broms  and  Flodin,  1988).  For  the  last  six  decades, 
the  literature  has  been  filled  with  theoretical  and  case  studies  on  the  topic  of 
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undrained  strength  prediction.  It  was  difficult  task  then  and  is  still  a  challenge  today. 
Historically,  the  problem  has  been  to  select  the  correct  bearing  Nk  factor  in  order  to 
calculate  the  undrained  strength  from  the  CPT  cone  resistance  measurement.  This  Nk 
factor,  which  ranges  from  8  to  25  (Aas,  Lacas.se,  Lunne,  and  Hoeg,  1986,  Lunne 
and  KJeven,  1981),  depends  on  overconsolidation,  clay  type,  sensitivity,  silt  content, 
reference  testing  device,  and  to  a  lesser  degree  on  overburden  stress.  A  new 
technique  is  developed  in  this  chapter  that  uses  both  CPT  measurements  to  directly 
predict  the  undrained  strength  without  the  need  to  estimate  the  Nk  factor. 


Historical  Means  of  Estimating  Nk  for  CPT  Prediction  of  Su 

The  classical  means  of  predicting  undrained  strength  using  the  CPT  uses  the  Nk 
factor  which  can  be  estimated  using  Table  8.1  (Olsen,  1994).  However,  this 
procedure  requires  prior  knowledge  of  the  in  situ  state  of  the  soil.  The  best 
procedure  for  calculating  Su  is  to  estimate  Nk  assuming  a  medium  stiff  normudy 
consolidated  condition  (i.e.  Nk  =  13)  then  calculate  Su  uri.ug  Equation  (8.1)  and  then 

Su 

calculate  the  — r  ratio. 
ov 


S  =  *£ _ (8.1) 

Equation  (8.1)  was  introduced  in  Chapter  4  during  the  development  of  the 

Su 

normalized  cone  resistance  formulation.  II  the  calculated  — r  is  0.29  to  0.33  and  Su 

av 

is  250  to  500  psf  (medium  stiff  condition)  then  the  soil  is  probably  normally 
consolidated  with  Nk  equal  to  approximately  13  (within  the  range  of  10  to  16).  If 
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u  .  . 

the  calculated  — r  is  high,  than  the  clay  is  probably  overconsolidated  or  a  stiffer 

av 

condition  prevails  which  requires  a  higher  Nk.  Iterations  might  be  required  until  the 
conditions  regarding  Nk  in  Table  8.1  are  matched.  A  few  reference  laboratory 
strength  tests  for  each  project  should  still  be  required.  This  discussion  has  shown  the 
Su 

importance  of  the  — 7  toward  estimating  the  Nk. 


Determining  Su]  using  CPT  Soil  Characterization  Chart 

The  normalized  undrained  strength  Su!  (i.e.  —  or  —7)  can  be  correlated  tc  CPT 

P  ov 

measurements  using  the  CPT  soil  characterization  chart.  This  section  will  introduce 
the  concept  of  predicting  Su]  based  on  bearing  capacity  formulation  and  the  next 
section  wili  describe  the  process  of  establishing  Sul  contours  on  the  CPT  soil 
characterization  chart. 

The  cone  resistance  bearing  formulation  (Equation  (8.2))  was  developed  in 
Chapter  5  and  is  in  terms  of  the  normalized  undrained  strength  (Sul),  cohesive 
bearing  factor  (Nk).  and  the  normalized  cone  resistance  (qcle): 

c  _  «c/e 


This  equation  can  be  rearranged  as  shown  below: 

Qrle  ~  SUI  l^k 


(8.2) 


(8.3) 


Equation  (8.3)  can  also  be  represented  in  a  graphic  form  based  on  the  CPT  soil 
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Table  8.1  Typical  Nk  values  for  estimating  the  unconsolidated  undrained  strength 
of  clays  using  the  CPT  cone  resistance  (Olsen,  1994) 


Clay  condition 

Unconsolidated  Undrained  (UU) 
Triaxial  Test 
as  the 

reference  strength  standard 

Nk 

range 

Nk 

average 

(typical) 

Normally  s  insolidatcd, 
normal  sevitivity, 
soft  to  me :  urn  stiff 

10  to  16 

n 
*  - 

Normally  o  isolidated, 
moderate^  nsitive, 
soft  to  ver  oft 

9  to  13 

ii 

Moderately  ver 

consolidates 

non-fissured 

15  to  20 

17 

Highly  over 
consolidated 
fissured 

17  to  23 

19 
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characterization  chart  as  shown  in  Figure  8.28.  The  qc]e  in  Equation  (8.3)  is  equal  to 
the  combined  effect  of  Sul  contours  and  Nk  contours  in  Figure  8.28.  For  a  given  Sul 
contour,  increasing  qcle  (vertical  axes)  will  correspond  to  an  increasing  contour 
of  Nk.  Increasing  the  overconsolidation  or  siit  content  will  also  increase  the  Su] 
contour  level  as  also  shown  with  arrows  and  annotations.  Therefore,  any  point  on 
the  CPT  Soil  Characterization  Chart  (for  example  point  "A")  has  a  Nk  contour  value, 
a  Su]  contour  value,  and  a  qcle  value  from  the  vertical  axis.  Contours  of  Sul  can  be 
established  directly  on  the  CPT  Soil  Characterization  Chart  using  laboratory 
unconsolidated  undrained  triaxial  strength  data  ((Sul)TxUU)  or  field  vane  shear 
strength  data  ((Sul)FV)  together  with  the  normalized  CPT  data.  There  is  no  need  to 
determine  the  Nk  contour  because  Sul  is  tire  ultimate  goal. 


Establishment  of  Su]  contours  on  the  CPT  soil  characterization  ' chart 

The  procedure  used  for  establishing  contours  of  normalized  undrained  strength 
(Sul)  for  clays  on  the  CPT  soil  characterization  chart  was  the  same  as  was  used  in 
the  last  section  for  prediction  of  the  SPT  blow  count.  Two  sets  of  SLll  contours  were 
established  using  measured  undrained  strength  data  from  the  database  developed  for 
this  research  program  (and  described  in  the  Appendix).  The  first  set  of  Sul  contours 
were  bused  on  data  from  laboratory  unconsolidated  undrained  triaxial  (TxUU)  tests 
and  the  second  set  were  based  on  data  from  field  vane  (FVanc)  shear  tests.  The  Sul 
contours  were  also  established  based  on  the  knowledge  (from  the  last  section)  that 
the  contours  should  increase  in  value  with  increased  qclc  and  f'slc  levels. 

All  measured  strength  data  from  the  database,  for  all  projects,  were  divided  into 
unconsolidated  undrained  triaxial  and  field  vane  shear  tests  and  then  further  divided 
into  6  strength  ranges  for  the  purpose  of  establishing  individual  correlation  contours 
on  the  CPT  soii  characterization  chart.  Best  fit  contours  of  normalized 
unconsolidated  undrained  triaxial  strengths  are  shown  in  Figure  8.29  to  Figure  8.35 
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Log  Normalized  Cone  Resistance 


Log  Friction  Ratio  (%) 


Figure  8.28  Establishing  Sul  (i.e.  c/p)  contours  on  the  CPI  soil  characterization 
chart,  based  on  bearing  theory  for  clay  penetration 
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based  on  the  minimum  inclusionary  AQI  levels  for  average  (Sul)TxUU  levels  of  0.15, 
0.25,  0.31,  0.39,  0.54,  and  2.0.  The  fmal  plot  having  all  the  best  fit  (Sul)TxUU 
contours  is  shown  in  Figure  8.36.  The  best  fit  contours  of  normalized  field  vane 
strength  are  shown  in  Figure  8.37  to  Figure  8.41  for  average  (Sul)rv  levels  of  0.25, 
0.31,  0.54,  0.80,  and  2.0.  The  final  plot  having  all  (Sul)FVane  contours  is  shown  in 
Figure  8.42. 

Several  examples  of  CPT  predicted  versus  measured  triaxial  unconsolidated 
mlramed  (UU)  tost  results  ((Su)Txirj)  are  shown  in  Figure  8.43  to  Figure  8.49. 

These  examples  represent  .sites  composed  of  soft  to  stiff  clay  including  several  sites 
composed  of  sandy  to  silty  clay.  There  are  also  several  examples  of  desiccated 
layers  with  properly  predicted  (Su)7xUU  levels 


Several  examples  of  CPT  predicted  versus  measured  field  vane  shear  test  results 


Vv.°u^I'vanc  )  ^  ^  ah  1  iap  a  tgu iv 


TTU  Avnmnlnr  rnnmrnnt  itarrf 

liivow  i  UJU1U  wi) 


soft  to  medium  stiff  day  sites  and  contain  a  few  examples  of  desiccated  strength. 
There  is  more  field  vane  strength  data  compare:1,  to  laboratory  triaxial  strength  data, 
probably  because  the  field  vane  test  is  easier  to  perform  and  less  expensive. 


The  contour'!  of  (Suj)lxUU  and  (Su;)rVai„  have  the  same  general  contouring 
characteristics,  namely  increasing  Su)  with  increasing  qcJc  and  fsle.  The  differences 
between  these  contours  is  important.  T  he  ratio  of  (Su,)TxUU  to  (Sul)1Vanc  is  the 
historic  Bjenrum  vane  shear  correction  fact',  r,  p.  (Bjerrum,  1972)  as  shown: 


b  =  - 


(s«! 


UU 


V 


TxUU 


TxUU 


FVarte 


(V 


V°v) 


FVane 


(8.4) 
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Figure  8.29  Predition  of  the  (Sul)-IxUU=0.1  contour  on  the  CPT  Soil 
Characterization  Chart 
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Normalized  Cone  Resistance  (aim  units) 
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Figure  8.30  Prediction  of  the  (S„,)7xUU=0.25  contour  on  the  CPT  Soil 
Characterization  Chart 
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Figure  8.31  Prediction  of  tlie  (Sul)lxUU~0.31  contour  on  the  CPT  Soil 
Characterization  Chart 
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Normalized  Cone  Resistance  (atm) 
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Figure  8.33  Prediction  of  the  (Sul)lxUU:=0.54  contour  on  die  CPT  Soil 
Characterization  Chart 
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Figure  8.34  Prediction  of  the  (Sul)TxUU=0.80  contour  on  the  CPT  Soil 
Characterization  Chart 
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Normalized  Cone  Resistance  (atm  units) 
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Figure  S.35  Prediction  of  the  (Sul)TxUu=2  contour  on  the  CPT  Soil  Characterization 
Chart 
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Figure  8.37  Prediction  of  the  (Sul)Fvanc  =  0.25  contour  on  the  CPT  soil 
characterization  chart 
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Figure  8.38  Prediction  of  the  (Su,),.VllK.  =  0.3  i  contour  on  tlie  CPI  soil 
characterization  chart 
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5'igure  8.39  Prediction  oi  the  (Sul),.-Vailc  =  0.54  contour  on  the  CPI  soil 
char;icte)  ization  chart 
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Prediction  of  the  (Sul),  Vanc  =  0.80  contour  on  the  CTT  soil 
characterization  chart 
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Normalized  Cone  Resistance  (Atm) 
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Figure  8.41  Prediction  of  the  (Stll),  Vanc  =  2  contour  on  tlic  CPI  soil 
characterization  chart 
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Figure  8.42  Contours  of  (Stl))iVanc  0-c-  (Su),  Vanc  /a'v  )  for  the  field  vane  shear  on 
the  CPT  soil  characterization  chart 
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figure  8.43  CPT  predicted  versus  measured  (Sui)lxlJU  for  silty  clays  from  a  eraser 
River  site,  Vancouver  International  Airport.  Canada  (data  from 
Konrad,  ct  al.,  1985) 
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Figure  8.44  CPT  predicted  versus  measured  (^u/h-.uu  for  silly  clay  ^iom  Minima 
Mcllosa,  Sweden  (40  km  north  of  Stockholm)  (data  from  Larsson  & 
Mulabdic,  i  94 1 ) 
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Figure  8.45  CPT  predicted  versus  measured  (Sul)lxUU  for  San  Francisco  Bay  mud 
at  tlic  Muni  Metro  turnaround  facility  (the  measured  strength  data  are 
located  ai  1 0  different  elevations)  (data  from  Dames  &  Moore,  1 992) 
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igurc  8.46  CPT  predicted  versus  measured  tSul)lxUU  for  the  Onsoy  site  Norway 
(NGI  research  site)  (the  measured  strength  data  are  at  four  different 
denths)  (data  from  Rad,  ct  al,  1985,  Gillespie,  ct  al.,  1985) 
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Figure  8.47  CPT  predicted  versus  measured  (Su,)rxUlJ  for  San  Francisco  Bay  mud 
at  Hamilton  AFB  (University  of  California  at  Berkeley  research  site) 
(data  from  project  code  ET-SE) 
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Figure  8.48  CPT  predicted  versus  measured  (Sul),xUU  from  Norrkoping,  Sweden 
(with  measurements  at  1 1  different  depths)  (a  Swedish  Geotechnical 
Institute  research  site)  (data  from  Larsson  and  Mu  abdic,  1991) 
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Figure  8.49  CPT  predicted  versus  measured  (Sul)TxUU  from  the  Hibernia  offshore 
investigation  (data  from  Konrad,  et  al.,  1991) 
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Figure  8.50  CPT  predicted  versus  measured  (Sy,),,^,  for  San  Francisco  bay  mud 
at  the  muni  metro  turnaround  facility'  (with  measurements  from  1 7 
different  depths)  (data  from  Dames  &  Moore,  1985) 
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CPT  predicted  versus  measured  (Sul),.vane  for  a  conference  pile  load 
test  at  Northwestern  University  in  1988  (with  measurements  from  12 
different  depths)  (data  from  Stratigraphies,  Inc,  1988) 
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Figure  8.52  CPT  predicted  versus  measured  (Sul), .-vanr  fur  a  foundation  day  ia) 
for  the  Nerlerk  Penn  liquefaction  slide  study 
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CPT  predicted  versus  measured  (Sul)l  vaiic  at  the  Fraser  River  Delta 
Bridge,  25  km  SW  of  Vancouver,  Canada  (with  measurements  from 
different  depths)  (Data  from  Crawford  &  Campanella,  1991) 
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Hus  [i  factor  is  used  to  reduce  the  measured  field  vane  shear  strength  to  an 
equivalent  undrained  triaxiai  test  strength  level  (i.c.  field  strength  value  for  slope 
stability  evaluation).  The  u  is  therefore  the  ratio  of  (S(:1)TxUli  contours  in  Figure  8.36 
to  the  (Sul)FV  contours  in  Figure  8.42  as  illustrated  in  Figure  8.57,  with  the  results 
shown  in  Figure  8.58.  The  p  contours  in  Figure  8.58  (i.e.  0.5  to  1.0)  have  the  same 
general  range  as  from  historic  observations  in  Figure  8.59  (i.e.  0.4  to  1.0).  The 
Bjerrum  p  factor  was  historically  indexed  to  Plasticity  Index  (PI)  as  shown  in 
Figure  8.59  (Bjerrum,  1972).  PI  is  now  used  only  to  index  overconsolidation 
character  as  shown  at  the  top  of  Figure  8.60  (Aas,  Lacasse,  Lurrnc,  and  Roeg,  1985). 
After  the  overconsolidation  character  is  estimated,  the  Bjerrum  ji  factor  is  then 
estimated  using  the  chart  ui  the  bottom  of  Figure  8.60.  This  figure  shows  that  tire 
u  factor  is  approximately  0.6  to  1.0  for  normally  consolidated  clay  and  0.35  to  1.0 
for  overconsolidated  clay.  The  contours  in  Figure  8.58  show  a  calculated  \x  factor  of 
about  0.9  to  1 .0  for  normally  consolidated  clays  decreasing  to  0..5  with  increasing 
overconsolidation  or  increasing  silt  content.  The  lack  of  held  vane  shear-  data  for 
highly  overconsolidated  clays  prevents  establishing  p.  contours  which  are  lower  than 
0.50. 

Conclusions 

T  his  section  introduced  a  new  technique  for  estimating  undrained  strength  of  clay 
which  is  appears  to  be  as  accurate  as  methods  which  require  estimating  the  cohesive 
Nj  bearing  factor.  These  correlations  were  shown  to  give  good  agreement  between 
measured  and  predicted  values  for  clays  around  the  world. 
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Calculated  contours  of  (Su)lxUU/(Su),.v  on  the  CIT  soil  characterization 
chart  by  taking  a  ratio  of  contours  on  Figure  8.36  to  Figure  8,4.: 
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Figure  8.60  Correlation  ot  (Sul)FV  to  PI  to  establish  the  overconsolidation 
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CPT  Prediction  of  Shear  Wave  Velocity  (Vs) 


Introduction 

As  noted  earlier,  accurate  prediction  of  the  shear  wave  velocity,  V, ,  using  the 
CPT  results  may  not  be  possible  because  Vs  is  a  low  strain  measurement  while  the 
CPT  is  a  high  strain  measurement.  Nonetheless,  there  have  been  several  attempts  to 
establish  useful  correlations.  Typically,  Vs  has  been  predicted  using  the  CPT  cone 
resistance  (Baldi,  et.al.,  1988)  or  more  recently  by  using  both  CPT  measurements 
(Olsen,  1988).  The  Olsen  (1988)  publication  developed  a  technique  for  CPT 
prediction  of  the  maximum  shear  modulus  (Gmax)  (Gmax  and  Vs  are  theoretically 
related).  This  section  will  extend  the  Olsen  (1988)  technique  for  CPT  prediction  of 
V£  based  on  improved  stress  normalization  techniques  (Chapters  6  &  7)  and  use  of 
the  new  larger  database  (Appendix  A).  Vs,  Graax,  (Gmax),,'and  Vsl  were  fully 
described  in  Chapter  7  during  the  development  of  the  shear  wave  velocity 
normalization  stress  exponent  (v). 


Minimal  Stratigraphic  Error  when  Vs  Measured  by  CPT  sounding 

There  are  two  general  means  of  measuring  shear  wave  velocity  in  the  field:  with 
borehole(s),  or  with  the  CPT  probe.  Borehole  based  shear  wave  measurements  are 
performed  cither  using  crosshole  (i.e.  boreholc-to-borehole)  or  downhole  (surface 
energy  source  and  geophones  in  the  borehole)  techniques.  The  CPT  based  shear 
wave  velocity  measurements  use  downhole  measurement  with  the  geophones  inside 
the  CPI  probe  (about  1  meter  up  from  the  cone  tip)  (FUGRO,  1980;  Campanella  and 
Robctson.  1984).  During  a  shear  wavi.  test  with  a  seismic  CPT  probe,  the  shear 
wave  tra-  els  down  through  the  soil  directly  adjacent  to  the  probe.  The  soil  tested  by 
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the  CPT  probe  is  the  same  soil  traversed  by  the  shear  waves,  therefore,  there  should 
be  little  if  any  stratigraphic  bias  when  establishing  correlations.  However, 
developing  a  correlation  is  still  difficult  because  shear  wave  velocity  is  a  low  strain 
measurement  and  the  CPT  measurements  are  high  strain  measurements. 


Historical  use  of  both  CPT  Measurements  to  Predict  Vs  or  Gmax 

The  first  technique  for  estimating  Vs  (i.e.  Gmax)  using  the  normalized  cone  and 
sleeve  resistances  is  shown  in  Figure  8.61  (Olsen,  1988).  The  contours  of  (Gmax)]  in 
this  figure  are  parallel  to  normalized  sleeve  resistance.  Gmax  appears  to  be 
proportional  to  the  sleeve  friction  resistance,  more  so  than  to  the  cone  resistance. 

The  proportionality  of  Vsl  (i.e.  (Gmax),  )  to  fs]e  was  believed  in  1987  because  both 
parameters  are  influenced  by  horizontal  stress.  It  -vTil  be  shown  in  a  later  section 
that  this  proportionality  is  more  likely  because  both  paianictoK-.  aie  dependent  on 
\oid  ratio  and  soil  type. 


CPT  Prediction  of  VsI  using  Both  CPT  Measurements 

'rite  procedures  for  establishing  contours  of  normalized  shear  wave  velocity  (V,.,) 
on  the  CPT  soil  characterization  chart  using  the  minimal  inclusionary  AQ1  technique 
is  the  same  as  used  for  establishing  the  Ns  and  S(ll  contours  in  previous  secii-  ns  of 
this  chapter.  All  shear  wave  velocity  data  from  the  database,  for  all  projects,  were 
divided  into  five  Vsl  ranges  for  the  purpose  of  establishing  individual  correlations 
contours  on  the  CPT  soil  characterization  chart,  flic  best-fit  contours  of  normalized 
shear  wave  velocity  using  the  highest  inclusionary  AQ1  are  shown  in  Figure  8.62  to 
Figure  8.6b  for  average  Vsl  levels  of  150,  400,  650,  900,  and  1200  fee u second.  Th 
final  best  fit  Vsl  contours  are  shown  in  Figure  8  67 
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Figure  8.62  Prediction  of  V  j  =  150  ft/sec  contour  on  the  CPT  soil  characterization 
chart 
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Several  examples  of  CPT  predicted  versus  measured  normalized  shear  wave 
velocity  (Vsl)  are  shown  in  Figure  8.68  to  Figure  8.76.  These  plots  represent  all  soil 
type  arid  all  relative  strength  levels.  The  important  observation  from  these  Figures  is 
that  the  CPT  predicted  Vs,  agrees  well  with  both  the  measured  values  and  'heir 
variation  with  depth.  It  is  also  important  to  note  that  about  80%  of  tire  shear  wave 
velocity  data  was  measured  using  the  seismic  CPT  probe  so  there  is  very  little 
stratigraphic  bias  error. 


Discussion  of  CPT  prediction  ofVs} 


For  clays  and  silts,  the  normalized  shear  wave  velocity  (V..,)  contours  in 
Figure  8.67  are  parallel  to  normalized  sleeve  friction,  f.,K  (fs)c  contours  are  shown  in 
Figure  8.18).  This  indicates  that  overburden  stress  similarly  influences  Vsl  and  fsk. 
However,  sleeve  friction  resistance  reftec’s  a  high-strain  strength  (Douglas  and 
Olsen,  1981)  while  the  shear  wave  velocity  represents  a  modulus  at  extremely  small 
strains  (10' 5  percent).  Shear  wave  velocity  is  indexable  to  void  rtio,  confining  stress 
and  soil  type  (Riehurt,  Hall,  and  Woods,  1970,  Seed  and  Mxiss,  1970).  Sleeve 
friction  resistance  approaches  steady  state  strength  (if  there  is  no  volumetric  shear) 
(Wald,  et.ah,  1991,  Castro,  et.al.,  1989).  Sleeve  friction  resistance  should  be 
indexable  to  the  same  properties  as  steady  state  strength,  namely;  void  ratio, 
confining  stress  and  soil  type  (Castro,  1969,  Schofield  and  Wroth,  1968).  Because, 
the  Vsl  and  fSiC  contours  are  parallel  (for  clayey  soils),  it  suggests  that  they  are  alsu 
dependent  to  the  same  combination  of  geotechnical  properties. 


The  Vsl  contours  bend  within  the  sand  zones  of  the  CPT  soil  characterization 
chart  (Figure  8.67)  foi  loose,  and  medium  dense  sands.  CPT  probing  of  these  sands 
causes  volume  change  due  to  the  high  bearing  stresses  which  results  in  a  denser  soil. 
For  loose  ch  av;-  sands,  the  measured  cone  resistance  therefore  reflect  a  denser  state 
than  the  in  situ  condition  because  of  grain  rearrangement,  grain  crushing,  and 
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dcnsification.  However,  for  dense  sand,  the  Vs]  contours  appear  to  be  parallel  to  fslc. 
Dense  sands  do  not  density  during  shear.  Therefore,  for  dense  sands,  the  in  situ  void 
ratio  is  approximately  equal  to  (or  even  greater  than)  the  void  ratio  surrounding  the 
sleeve  and  cone  units  during  penetration  The  resultant  is  a  Vs)  contour  which  is 
approximately  parallel  to  a  fs)c  contour  from  clay  to  dense  sand  as  shown  in 
Figure  8.67  for  Vsl=1200  ft/second. 

The  Vsl  contours  in  Figure  8.67  can  also  be  generalized  as  shown  in  Figure  8.77 
as  a  means  to  explain  densification  effects  for  loose  to  medium  dense  sand  during 
penetration.  A  line  of  constant  normalized  sleeve  friction  resistance  is  shown  as  Line 
CF  in  Figure  8.77.  Line  CF  would  be  parallel  to  the  Vs,  contour  if  there  was  no 
volume  change  in  the  clean  sands.  However,  line  VF  is  a  typical  Vsl  contour.  The 
difference  between  Lines  CF  and  VF  could  be  the  effect  of  sand  volume  change 
during  probing  resulting  in  a  denser  soil  for  loose  to  medium  dense  sands.  Point  N 
is  along  a  constant  fsle  contour  (Line  CF)  representing  no  volume  change  for  clean 
loose  sand.  Point  S  represents  a  denser  sand  condition  due  to  probing  and  a  higher 
measured  fsle  (and  qclc)  because  the  sand  surrounding  the  sleeve  unit  is  denser. 
Therefore,  CPT  probing  of  clean  loose  to  medium  dense  sand  produces  a  denser 
condition  than  exists  in  situ  which  results  in  an  explainable  bending  of  the  Vsl 
contours. 
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Figure  8.77  Illustration  why  the  predicted  shear  wave  velocity  contours  on  the  CPT 
Soil  Characterization  chart  are  curved 
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Chapter  9 


Conclusions 


The  primary  objective  of  this  research  was  to  develop  techniques  and  correlations 
for  predicting  geotechnical  properties  using  CPT  cone  resistance  and  sleeve  friction 
resistance  measurements.  Stress  normalization  provided  the  means  for  taking 
confining  stress  dependence  into  account  when  predicting  geotechnical  properties 
using  CPT  measurements. 

A  new  concept,  the  Stress  Focus,  was  discovered  and  confirmed  which  provides 
a  basis  for  understanding  strength  behavior  (e.g.  friction  angle,  cone  resistance,  etc.,) 
as  a  function  of  confining  stress.  This  study  demonstrated  that  the  relationships  of 
sand  friction  angle  to  confining  stress  for  different  initial  densities  converge  to  a 
Stress  Focus  at  high  pressure,  where  the  strength  behavior  is  similar  to  that  of  a 
sedimentary  rock.  This  convergence  to  a  Stress  Focus  was  also  confirmed  using  the 
CPT  cone  resistance  measurement.  The  paths  of  convergence  to  the  Stress  Focus  are 
exponentially  related  to  overburden  stress  and  can  be  represented  as  straight  lines  on 
a  log-log  plot  of  strength  versus  vertical  effective  stress.  A  sand  at  a  given  relative 
density  can  be  represented  by  a  straight  line  on  this  log-log  plot  with  all  relative 
density  lines  converging  to  the  Stress  Focus.  The  slope  of  this  line,  termed  the  stress  * 
exponent,  is  inversely  proportional  to  sand  initial  relative  density. 

The  strength  of  dense  sand  is  strongly  influenced  by  dilation  effects.  Dilation 
effects  for  a  dense  sand  will  decrease  with  increased  vertical  effective  stress  until  the 
Stress  Focus  is  reached  where  its  behavior  is  similar  to  that  of  a  sedimentary  rock. 
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The  convergence  of  relative  density  lines  to  the  Stress  Focus  for  sand  is  caused  by 
consolidation  (i.e.  decreasing  void  ratio)  which  is  reflected  in  numerous  geotechnical 
properties  such  as  Mohr  envelope  curvature,  grain  crushing,  compressibility,  etc., 
Sands  of  all  relative  densities  will  have  approximately  the  same  bulk  density 
(i.e.  void  ratio)  at  the  Stress  Focus  and,  as  a  result,  approximately  the  same  strength 
behavior.  The  Stress  Focus  should  therefore  be  considered  a  fundamental 
geotechnical  property.  The  Stress  Focus  concept  replaces  the  variable  soil  type 
dependent  stress  exponent  technique  for  normalizing  cone  resistance  previously 
developed  by  the  author  (Olsen,  1984,  1988). 

The  Stress  Focus  concept  was  confirmed  using  historic  high  pressure  triaxial  test 
data  in  Chapter  2  (collectively  shown  by  Figure  2.5).  Sand  strengths  are  shown  to 
converge  to  a  Stress  Focus  at  a  vertical  effective  stress  of  approximately  100 
atmospheres.  A  simple  yet  unique  CPT  cone  resistance  normalization  formula, 
which  accounts  for  exponential  stress  dependence  of  the  tip  stress,  has  been  derived 
in  Chapter  4.  The  Stress  Focus  for  CPT  cone  resistance  is  demonstrated  using  CPT 
laboratory  chamber  data  in  Chapter  5. 

The  Stress  Focus  location  (i.e.  vertical  effective  stress  and  strength  level)  was 
shown  to  be  soil  ty  pe  dependent  using  field  CPT  data  from  uniform  soil  layers 
(Figure  6.24).  For  clay,  the  Stress  Focus  occurs  at  a  vertical  effective  stress  of 
approximately  9  atm.  For  sands  the  Stress  Focus  occurs  at  a  vertical  effective  stress 
of  approximately  70  to  200  atm. 

A  technique  to  estimate  the  stress  exponent  required  for  cone  resistance 
normalization  was  also  developed  using  CPT  data  from  uniform  soil  layers. 

Contours  of  cone  resistance  stress  exponent  were  established  on  the  CPT  soil 
characterization  chart  (Figure  6.16).  These  stress  exponent  contours  also  support  the 
validity  of  the  CPI  soil  characterization  chart  (Olsen,  1988)  to  characterize  soil  in 
terms  of  soil  type  and  state. 
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The  constant  drilling  mud  height  (i.e.  constant  pressure)  used  in  SPT  cl  amber 
tests,  for  all  confining  stress  levels,  reduces  the  confining  stresses  next  to  the  SPT 
sampler  (Chapter  7).  All  SPT  chamber  tests  in  the  past  have  had  this  problem.  This 
reduced  confining  stress  results  in  a  SPT  based  stress  exponent  which  is  too  low  if 
derived  from  the  results  of  chamber  tests.  The  SPT-N  to  N,  normalization  concept, 
developed  over  10-years  ago,  is  based  on  SPT  chamber  test  results  and  therefore  uses 
stress  exponents  that  are  too  low.  The  stress  exponent  for  SPT  borehole  applications 
were  shown  to  be  equal  to  the  CPT  cone  resistance  stress  exponents.  Therefore,  the 
CPT  determined  cone  resistance  stress  exponent  should  be  used  for  SPT 
normalization.  To  achieve  the  best  SPT  normalization,  therefore,  requires  a  nearby 
CPT  sounding  where  the  cone  resistance  stress  exponent  can  be  estimated. 

The  stress  exponent  for  shear  wave  velocity  is  shown  to  be  soil  type  dependent, 
and  it  can  be  approximated  as  45%  of  the  CPT  cone  resistance  stress  exponent  as 
shown  in  Chapter  7. 

Developing  CPT  correlations  to  geotechnical  properties  required  a  large  database 
of  CPT  and  tested  soil  sample  data  (described  in  the  Appendix).  The  largest  error 
during  con-elating  of  CPT  measurements  to  nearby  borehole  soil  samples  is 
geological  change  of  soil  type.  A  change  of  soil  type  is  considered  a  bias  condition 
because  it  will  skew  the  data  trend.  Bias  data  must  therefore  be  discounted  when 
developing  correlations.  A  quality  index  was  developed  in  an  attempt  to  account  for 
possible  soil  type  difference  between  CPT  soundings  and  boreholes — namely  the 
Academic  Quality  Index  (AQI).  The  overall  AQI  accounts  for  possible  stratigraphic 
change  and  CPT  measurement  quality.  Establishing  predictive  contour  trends  on  the 
CPT  soil  characterization  chart  was  accomplished  with  a  new  teclmique  that  uses  the 
AQI  quality  index  to  account  for  bias  error.  Correlations  were  established  for  the 
following  normalized  geotechnical  properties;  SPT  blow  count,  undrained  cohesive 
strength  from  the  unconsolidated  undrained  triaxial  test,  undrained  cohesive  strength 
from  the  field  vane  shear  test,  and  the  shear  wave  velocity. 
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The  contours  of  normalized  SPT  blow  count  on  the  CPT  soil  characterization 
chart  have  a  predictable  trend  (Figure  8. 1 7).  In  loose  to  medium  dense  sands,  the 
SPT  sampler  is  primarily  resisted  by  end  bearing  force  (while  the  sampler  side 
friction  force  is  minor).  Cn  the  other  hand,  for  dense  (and  overconsolidated)  sands 
and  clays,  SPT  sampler  side  friction  force  dominates.  This  relative  contribution  of 
SPT  sampler  end  bearing  to  side  friction  forces  for  different  relative  densities  is 
confirmation  of  work  by  Schmertmann  (1979a).  The  technique  for  predicting  the 
SPT-N  values  based  on  the  CPT  cone  resistance  and  sleeve  friction  resistance  models 
the  forces  on  the  SPT  sampler  and  demonstrates  the  potential  for  reliable  prediction 
of  SPT  N  values  from  knowledge  of  CPT  cone  and  sleeve  friction  resistances. 

C  Su 

Contours  of  normalized  undrained  cohesive  strength  (i.e.  —  ,  — r,  or  S  ,)  for 

P  ov 

unconsolidated  undrained  triaxial  tests  (Figure  8.36)  and  field  vane  shear  tests 
(Figure  oA2)  were  established  on  the  CPT  soii  characterization  chart.  Historically, 
the  Nk  (i.e.  ratio  of  net  cone  resistance  to  measure  undrained  strength)  must  be  either 
be  estimated  or  developed  as  a  site  specific  value.  An  estimation  of  the  Nk  bearing 
factor  is  not  required  with  this  new  technique  for  CPT  prediction  of  clay  undrained 
strength.  The  ratio  of  the  (Su)TxUU  versus  (Su){  v  is  the  historic  Bjerrum  p  correction 
factor  used  for  reducing  the  field  vane  shear  test  results  to  unconsolidated  undrained 
triaxial  (TxUU)  strength  levels.  The  calculated  p  (Figure  8.58),  based  on  these  new 
correlations,  has  the  same  range  published  by  NGI,  namely  0.9  to  1  for  normally 
consolidated  clays  and  decreases  to  0.5  for  overconsolidation  (or  increased 
silt  content). 

Contours  of  nonnalized  shear  wave  velocity  were  also  established  on  the  CPT 
soil  characterization  chart.  These  contours  are  parallel  to  the  normalized  sleeve 
friction  resistance  contours  when  there  is  no  volume  change  during  shear  (e.g.  clays, 
clayey  silts,  clayey  sands,  dense  sands,  etc.,).  It  appears  that  the  friction  sleeve 
resistance  and  shear  wave  velocity  are  influenced  by  the  same  geoteclinical  properties 
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and  in  the  same  proportion.  The  normalized  shear  wave  velocity  contours  were  also 

observed  to  bend  within  the  loose  to  medium  dense  clean  sand  area  of  the  CPT  soil 

characterization  chart  because  volume  change  during  shear  causes  a  denser  state. 

The  predicted  shear  wave  velocity  agreed  well  with  the  measured  values. 

Future  Research 

•  Integration  of  the  Stress  Focus  concept  into  critical  state  soil  mechanics. 

•  Continue  to  evaluate  chamber  test  results  toward  better  understanding  and 
definition  of  the  Stress  Focus. 

•  Evaluate  man-compacted  soils  (e.g.  earth  dams,  fills,  etc.,)  toward  establishing 
better  techniques  for  stress  normalization  and  prediction  of  geotechnical  properties. 

•  Investigate  the  soil  property  and  bulk  density  at  the  Stress  Focus. 

•  Seek  more  verification  of  the  CPT  predicted  geotechnical  properties. 

•  Continue  to  develop  the  AQI  as  a  subjective  index  for  excluding  bias  error. 

•  Develop  a  computer  technique  for  locating  the  Stress  Focus  using  several  relative 
density  groups  (this  is  an  iterative  procedure  which  searches  for  the  optimum 
location  for  the  Stress  Focus  by  minimizing  the  relative  density  trend  variances). 

•  Develop  improved  techniques  for  calculating  the  CPT  stress  exponents  using  field 
data. 
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Appendix 


The  Database  Contents 


Introduction 

Collecting  the  CPT  and  boring  data  represented  the  single  largest  effort  (and 
longest  time  effort)  for  the  purpose  of  developing  CPT  correlations  of  geotechnical 
properties.  There  were  three  major  sources  for  this  data,  namely,  1)  the  author's 
project  files,  2)  requests  for  data  by  letter,  and  3)  requests  for  data  after  a  lecture  was 
presented.  Overall,  the  best  CPT  and  boring  data  originated  from  the  author's  project 
files.  This  data  base  represents  the  largest  coherent  accumulation  of  CPT  cone  and 
sleeve  with  boring/laboratoiy  project/research  data  in  the  world  (1994).  This 
collection  of  data  also  contains  two  types  of  data:  field  CPT/boring  data  and 
laboratory  large  diameter  chamber  test  data.  Most  of  the  collection  emphasis  was 
directed  toward  field  CPT/boring  data. 


Requesting  Data  by  Letter 

From  1986  to  1988,  about  120  letters  were  sent  to  government,  utility  and 
consulting  firms  requesting  CPT  and  boring  data.  Only  1 5%  of  the  requests  for  data 
resulted  in  data  and  only  5%  of  the  total  database  came  from  these  blind  requests  for 
data.  Only  60%  of  the  data  received  during  this  research  program  was  useable.  Less 
than  15%  of  the  geotechnical  engineering  publications  which  describe  projects  having 
CPT  data  actually  show  a  CPT  sounding  log.  When  a  CPT  sounding  log  was 


presented  in  a  publication,  less  than  10%  of  the  total  also  show  sleeve  friction 
measurements.  Therefore,  only  a  small  fraction  of  this  database  originated  from 
publications. 


Requesting  Data  by  Presenting  Lectures 

Approximately  half  of  the  geotechnical  data  came  from  professional  contacts  and 
as  a  result  of  presentation  lectures  on  CPT  technology.  A  total  of  23  lectures  were 
presentations  by  the  author  to  consulting  engineering  companies  and  state/federal 
organizations  from  1988  to  1992.  Engineers  typically  felt  obligated  after  a  lecture 
was  presented  to  search  for  and  copy  at  least  one  geotechnical  project  for  this 
research  program. 


The  Data  in  the  Database 

This  database  contains  approximately  670  CPT  records  and  580  borings  from  90 
projects.  The  CPT  records  represent  approximately  51,000  ft  of  data  and  t  :re  are  at 
least  1,200  boring  soil  samples  with  a  total  of  approximately  8,100  laboratory  and 
field  test  values.  The  following  is  a  partial  list  of  the  types  of  data  in  the  database: 
CPT  measurements  (cone  resistance,  sleeve  friction  resistance,  and  dynamic  pore 
pressure),  SPT  blow  count,  measured  shear  wave  velocity,  water  content,  plastic 
limit,  liquid  limit,  total  density,  percent  passing  #200  sieve,  D50,  field  vane  shear 
strength,  laboratory  triaxial  strength  test  results  (and  testing  method),  Unifield  Soil 
Classification  System  designation,  word  descriptor  based  soil  classification, 
consolidation  parameters  (e.g.  Cc,  C5,  OCR,  Pmax,  etc.,),  void  ratio,  etc., 
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Test  Chamber  Data 


Data  from  several  large  diameter  test  chambers  are  also  included  in  the  database 
for  establishing  new  CPT  and  SPT  normalization  techniques.  The  CPT  chamber  data 
was  provided  by  Jamiolkowski  (1988).  SPT  chamber  data  by  Bieganousky  and 
Marcuson  (1976,  1977)  and  Gibbs  and  Holtz  (1952)  were  also  included  in  this  data 
base  collection.  A  complete  listing  and  description  of  the  CPT  and  boring  data  base 
are  shown  in  Table  A-l.  The  soil  conditions  column  provides  a  general  descriptor 
only  for  comparison  between  sites. 


Data  Sensitivity 

Most  of  the  data  used  in  this  research  program  originated  from  military, 
sensitive  security,  or  sensitive  private  projects  where  the  data  source  must  be  kept 
confidential.  Much  of  the  private  project  data  (classified  as  sensitive  private  sites) 
were  provided  by  geotechnical  engineering  consulting  firms  with  the  understanding 
that  the  data  and  source  would  be  kept  confidential.  In  all  cases,  the  data  could  be 
used  to  establish  data  correlations  and  the  data  points  on  the  correlation  plots  could 
be  published;  however,  the  provider  would  not  allow  the  site  name  to  be  published 
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Table  A-l  Summary  of  CPT  projects  contained  in  the  data  base 
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